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The specific ionization of 120 cosmic-ray tracks has been measured by counting droplets in 
photographs of delayed expansion cloud chamber tracks. For values of J/p between 10° and 
2 < 10° the specific ionization as a function of Hp is found to be in good agreement with theory. 
The predicted minimum ionization for electrons of Hp~10* (~2 Mev) is verified. These results 
are applied to the calculation of the mass of heavily ionizing cosmic-ray particles as observed 
by us and by others. The masses of all heavy cosmic-ray particles thus far reported, with one 
exception, do not seem to be in serious disagreement with a unique mass which lies in the 





limits (200+50) nto. 


: iw obtain values for the specific ionization 
due to cosmic-ray particles several hundred 
cloud chamber photographs have been taken. 
The expansion was delayed to make possible the 
counting of the individual droplets in the tracks. 
The cloud chamber was 18 cm in diameter and 
was filled with nitrogen at about 1.5 atmospheres 
pressure. It was placed in a magnetic field of 
800 gauss. The chamber was illuminated by a 
mercury arc lamp, operated with a fraction of 
an ohm series resistance and flashed on 220 volts 
d.c. for 0.1 second. The expansion was delayed 
for approximately 0.5 second by a cam mecha- 
nism set into operation by the tripping of a pair 
of Geiger counters. 

The tracks were broadened by diffusion of the 
ions in the time between the passage of the 
ionizing particle and the expansion of the 
chamber (the clearing field was shorted out by 
the tripping of the Geiger counters). In spite of 
the breadth of the tracks it was possible to 
obtain quite accurate values for the curvature 
by measuring the coordinates of the estimated 


center of the band of droplets when observed 
with a measuring microscope. 

The turbulent 
checked by measuring tracks when the magnetic 
These were found to be straight 


absence of distortions was 
field was zero. 
within the limits of experimental measurement 
except in the region less than two cm from the 
chamber walls. The droplets were easily counted 
as can be seen from the examples shown in 
Fig. 1. All to 0°C and 
760 mm of mercury pressure. 

One hundred twenty tracks that 
sharp focus, so that the number of droplets per 


cm could be easily counted, were chosen for this 


results were reduced 


were in 


preliminary survey. They were divided into 
seven //p groups. In obtaining the mean number 
of droplets per cm in each group (plotted in 
Fig. 2) the mean number per cm in each track 
was weighted proportional to the total number 
of droplets in the track. The vertical lines 
represent the probable errors. In the lowest J/p 
group the 
uncertain and the variation of ionization with 


radius of curvature is somewhat 
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Hp is rapid. In this region large deviations from 
the mean are to be expected. 


DISCUSSION 


The variation of the specific ionization in 
hydrogen is given theoretically by the formula! 


A 0s 

iJ log k+log ——-¢| (1) 
where k=1.6X10° for hydrogen. This is a 
function which varies nearly as 1/8? for B<0.9, 
then passes through a minimum at 8~0.97, 
and increases with increasing 8. For electrons 
the value of 8=0.97 corresponds to about 2 Mev. 
A smaller value of k must be used if this equation 
is to be applied to ionization in nitrogen. We 
have chosen k=2X10* as an appropriate value. 
Since k enters only in the logarithmic term the 
choice of a value for & is not critical. 

Our measurements do not give the primary ioni- 
zation but the “‘probable’’ ionization,’ i.e., the 
primary ionization plus the ionization produced 
by secondaries of energy lower than a certain 
critical energy. Secondaries of higher energies 
produce unresolvable The 
manner in which the theoretical formula must be 
modified to be applicable to probable ionization 
is discussed by Oppenheimer.’ Above the mini- 
mum point the probable ionization increases less 
rapidly than the primary ionization, but the 
changes in the shape of the curve below the 
minimum point are quite small. 

The continuous curve in Fig. 2 is Eq. (1) with 
k=2X10* and with the minimum ordinate 
arbitrarily adjusted to make the curve fit our 
observed points. The agreement of the theoretical 
curve with the experimental data is quite good. 
It is significant that our data differ from the 
classical 1/v* variation of ionization with velocity, 
as shown by the dotted curve in Fig. 2. Our 
observations show the actual increase in ioniza- 
tion with Hp, as predicted by theory. 

A point of interest in connection with all our 
measurements is the ratio of probable to primary 


clusters of ions. 


1Cf. Bethe, Handbuch der Physik, Vol. 24: 1. 
? Williams, Proc. Roy. Soc. A135, 108 (1932). 
3’ Oppenheimer, Phys. Rev. 47, 44 (1935). 
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ionization. A number or observers** have meas- 
ured the primary ionization of fast electrons by 
counting the number of distinct groups of 
droplets in sharp tracks, each group correspond- 
ing to one primary ion pair. For electrons with 
velocities comparable to those in the lower energy 
region of our observations all the observers 
report the number of primary ion pairs per cm 
(in air, nitrogen or oxygen at 0°C and 760 mm 
mercury pressure) to be about 20 or 25. By 
similar measurements on a few sharp tracks of 
energies in the minimum part of the curve we 
find between 14 and 18 primary ion pairs. Our 
total probable ionization is 50 ions or 25 ion 
pairs, which gives a ratio of probable to primary 
ionization of a little less than two. By counting 
the total number of droplets in diffuse tracks 
Wilson gives 1.8 for the same ratio. He says 
this value is low, however, because the larger 
groups of droplets were not resolved well enough 
to count. For lower energy electrons and more 
diffuse tracks he gives 3.5 for the ratio. 

From the observed range of electrons, and the 
variation of range with J/p one can calculate 
the average energy loss per cm for an electron of 
given J/p. If one takes the average energy per 
ion pair as 32 volts, this calculation indicates an 
average density of total ionization about twice 
that given by our droplet counts. 


CALCULATION OF MAss 


In cloud chamber tracks of moving charged 
particles the observable quantities are ionization, 
range and radius of curvature (if in a magnetic 
field), and the rates of change of these quantities 
with distance. These quantities depend on the 
charge, mass and velocity of the particle. It is 
possible to represent very closely the relation- 
ships between these different quantities and 
their rates of change on a nomograph of straight 
parallel lines. Such a nomograph is shown in 
Fig. 3. m/mp is the ratio of rest mass to rest 
mass of the ordinary electron, D is the ratio of 
specific ionization to specific ionization at the 
minimum point, D’ is the same ratio obtained 
by use of the 1/v? law. R is the range in cm of 


4C. T. R. Wilson, Proc. Roy. Soc. A104, 191 (1923). 


‘Williams and Terroux, Proc. Roy. Soc. A126, 289 
(1930). 
6 Skramstad and Loughridge, Phys. Rev. 50, 677 (1936). 
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Fic. 1. Sections of electron tracks broadened by diffusion 
of the ions in delayed expansions. 


air at 760 mm of mercury pressure and 15°C. 
d(l[p)/dR is the rate of change of Hp with 
distance. Any straight line drawn across the 
nomograph gives consistent values for all the 
quantities plotted. If the data include any two 
of the variables the mass of the particle in 
question is given at once. 

From an inspection of the nomograph one 
sees that //p and range afford the most accurate 
method of determining the mass, provided all 
the data are given with the same degree of 
precision. JJp and D also is a good way to 
determine the mass, for a relatively large change 
in D produces a relatively small change in m/mpo. 
D and range is a poor way to find m/mzy for the 
value of D must be accurately given in order to 
determine m/mp within relatively small limits. 
This last method has the advantage that it does 
not involve a measurement of curvature, and 
so is free of distortions which may produce 
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spurious curvatures. A significant fact is that a 
value of J/p and a minimum value for the range 
determine an upper limit on the mass. 

The ionization is theoretically independent of 
the mass, is proportional to the square of the 
charge on the moving particle, and depends on 
the particle’s velocity as given by Eq. (1). The 
fact that all particles with 8 near unity have 
approximately the same ionization indicates 
that all cosmic-ray particles carry the charge e. 
We have defined the relative specific ionization 
D as the ratio of the specific ionization at a 
given 8 to the ionization at the minimum, i.e. 
at B~0.97. These values can be calculated from 
Eq. (1) when the specific ionization at the 
minimum is known. The experimental value of 
D is found by dividing the number of droplets 
per cm in the track by our observed minimum 
value of 50. In this way the value of 8 for any 
track can be found. Then the relationship 
between m, J/p, and 6 for a charged particle 
moving in a magnetic field is given by the 
equation 

m /my = IT p(e/myc)(1 — B*)* /B, 


where #2) is the mass of the ordinary electron. 
Since the loss of energy per cm of path of a 
charged particle depends only on its velocity 
one can write R=kmf(v). R is the range in cm 
of air, m the particle’s rest mass and f(v) some 
function of velocity. This is true for relativistic 
as well as nonrelativistic energies. Then for all 


HP 


Fic. 2. Specific ionization vs. Hp. I is the number of 
droplets per cm in nitrogen at 0°C and 760 mm pressure. 
The continuous curve is the theoretical equation (1). 
The dotted curve is the inverse square law, The vertical 
lines represent the probable errors of the measurements. 
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particles of the same velocity R,/Re=m,/mo. 
When the velocities are the same we also have 
(IIp);/(Tp)2=m,/me. Then we can use the known 
ranges of protons’ to find the range R,; of a 
particle of mass m, and //p=(J/p),. Conversely, 
the range of a particle and its //p determine 
its Mass. 

It is an experimental fact that the relation 
between range and velocity may be closely 
expressed as R=kv" or R=kR(I[p)"/m""'. From 
this one gets 


R=p [n(dp dR) }. 


Thus p and dp/dR is sufficient information to 
determine the range the mass. The 
exponent ‘‘2’’ has a value varying from 3.45 to 
3.70 for 10*<J/p <10° and 80<m/my <300. 


and so 
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Fic. 3. Nomograph giving mass in terms of variables 
which can be determined from cloud chamber data. 
m,mo is the ratio of mass to the mass of the ordinary 
electron, D is the ratio of specific ionization to specific 
ionization at the minimum, D’ is the same ratio from the 
1/v? law, R is the range in cm of air at 15°C and 750 mm 
pressure and d(Hp)dR is the rate of change of Hp with 
distance. Any straight line drawn across the nomograph 
gives consistent values for all the variables plotted. 


Cf. M. S. Livingston and H. A. Bethe, Rev. Mod. 
Phys. 9, 269 (1937). 
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Fic. 4. A photograph of the heavy track previously 
reported by us. The density of ionization in this track is 
about 5.5 times that of an ordinary cosmic-ray track. The 
value of Hp=1.5 X 10° gives a mass of 250 times the mass 
of an electron. 

It is also possible to find m/mp» in terms of 
range and ionization (i.e., in terms of range and 
velocity), and thus eliminate the necessity of 
measuring curvature. But, as was pointed out 
in connection with the this is 
usually not a good way to determine 7/ mp. 


nomograph, 


APPLICATION TO HEAVY TRACKS 


We have previously reported® a track of a 
heavily ionizing particle, reproduced here in 
Fig. 4. Because of an error in our calculations 
the /7p was incorrectly given. The correct value 


8 Corson and Brode, Phys. Rev. 53, 215 (1938). 
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is JJp=1.5X10°. The number of droplets per 
cm is 265, which is about 5.5 times the minimum 
density of normal tracks. By the use of this as 
the value of D on the nomograph a value of 
m/mo= 250 is found. 

Another track of interest is one already 
published.’ This track has an I/p=5.5 X10* and 
an observable range of 18 cm (in nitrogen at 0°C 
and 760 mm of mercury pressure). These are 
sufficient data to place an upper limit of 200m» 
on the mass. The track is out of focus slightly 
which makes a count of the specific ionization 
impossible. If the ion density were 10 times 
that of normal tracks the mass would be 1257». 
The fact that undeflected tracks appear in the 
same region in the chamber makes it improbable 
that the curvature was produced by turbulence. 

The track reported by Street and Stevenson!® 
has an Hp=9.6X10' and an ion density six 
times that of normal thin tracks. Using the 1/v* 
law for ionization, D’ in Fig. 3, Street and 
Stevenson gave a mass of 130. When one 
uses, instead of D’, the column D in Fig. 3 the 
mass is found to be 160m». 

Two unusual tracks have been reported by 
Anderson and Neddermeyer.'' One of these 
tracks has an apparent J/p=5.5X10* and has a 
range of about four cm. This gives a mass of 
about 350m. The other track has an apparent 
Hp=1.4X105 and a range greater than five cm 
(the particle apparently passed out of the 
illuminated part of the chamber). This puts an 
upper limit of about 1000m, on the mass of the 
particle. 

Ruhlig and Crane” reported a track with a 
curvature of nine cm in a magnetic field of 
2850 gauss, and a rate of change of p with 
distance of 0.6+0.4. This gives a mass of 


® Brode and Starr, Phys. Rev. 53, 3 (1938), Fig. 4. 

10 Street and Stevenson, Phys. Rev. 52, 1003 (1937). 

4 Anderson and Neddermeyer, Phys. Rev. 50, 263 
(1936), Figs. 12 and 13. 

? Ruhlig and Crane, Phys. Rev. 53, 266 (1938). 
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(120+30)m. However, the track appears to 
have a range of at least five cm, which would 
place an upper limit of about 110 on the mass. 

Nishina, Takeuchi and Ichimiya™ reported a 
particle whose J/p changed from 7.4X10° to 
4.9X10° in passing through 3.5 cm of Pb. If 
one assumes that this particle is losing energy 
only by ionization then the equivalent thickness 
of air' is about 1.8X10° cm, or a value of 
d(IIp)/dR=1.4. This gives a mass of about 
200m. 


DISCUSSION 


Our measurements indicate that the variation 
of specific ionization with Hp (for JZp's between 
10° and 210°) is in good agreement with the 
theoretical equation (1). This agreement includes 
the verification of the minimum of ionization at 
IIp~10*, as predicted by the theory. In de- 
termining the velocity (and so the mass) of 
heavily ionizing particles one should use, not 
the 1/v? law, but Eq. (1). By use of Eq. (1), 
and consideration of the errors to which the 
experiments are subject, it seems that all the 
heavy tracks reported thus far, with the ex- 
ception of Ruhlig’s and Crane’s, are not in 
serious disagreement with a unique mass lying 
within the limits (200+50)m. The disagreement 
between the different values given for the mass 
of the heavy particle may be due to errors in 
determining the radius of curvature in the 
magnetic field. When the curvature is small it 
may be influenced considerably by scattering 
and turbulence. The range can be determined 
with a relatively high degree of accuracy. In 
nearly every case one observes a minimum 
range, thus placing an upper limit on the mass. 

We are indebted to Professor J. R. Oppen- 
heimer for his helpful discussions and suggestions. 


8 Nishina, Takeuchi and Ichimiya, Phys. Rev. 52, 1198 
(1937). 
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The spacing of the levels in the fine structure of alpha- and beta-ray processes and the 
existence of metastable nuclear states (isomers, isobars) are in contradiction with the existence 
of low lying levels corresponding to the rotation of the nucleus as a whole. The exchange of 
the nuclear constituents effected through rotation, together with the fact that the particles are 
not rigidly bound to equilibrium positions in the nucleus, will in some cases forbid, in other 
cases perturb the lowest levels, and cause the first state of excitation to lie considerably higher. 
Simple models illustrating these effects are discussed in II. In III an estimate is made for the 
position of the lowest excited level for heavy nuclei. It is found to vary inversely with the mass. 


the resulting nuclei have J=0, it should be 
possible to observe not only the process J=0 to 


MAY 15, 1938 PHYSICAL 
On the Rotation of the Atomic Nucleus 
E. TELLER, George Washington University, Washington, D. C. 
AND 
J. A. WHEELER, University of North Carolina, Chapel Hill, N. C. 
(Received March 23, 1938) 
I. INTRODUCTION 
CCORDING to Bohr and Kalckar' the 


lower levels of the heavier nuclei can be 
systematized if the nucleus is compared to a 
droplet. The low lying levels which are of 
importance in nuclear reactions will then be the 
vibrations of this droplet. Still lower levels are 
obtained, however, if the rotations of the droplet 
are taken into consideration. If we take as an 
example the nucleus Pb***, assume a uniform 
distribution of mass throughout the nucleus and 
take its radius to be 12X10-" cm, the moment 
of inertia 2.0010-** g¢ cm? will be obtained. 
If we furthermore assume that the lowest 
rotational state of this nucleus corresponds to 
J=0, then the pure rotational energies corre- 
sponding to J=1, 2, 3, will be 3.3, 10, 20, 
+++ kilovolts, approximately. 

According to a proposition first put forward 
by Thibaud? the rotations of the nucleus might 
be used to explain fine structure observed in the 
emission of alpha-particles. The energy difter- 
ences of rotation, however, are smaller than those 
observed in the fine structure. It is indeed 
reasonable to assume that, in general, nuclei 
will not have excited states which lie only a few 
thousand volts above the ground level. 

Such extremely low levels would give rise to 
difficulties in the following ways: 

(1) If in the process of alpha-particle emission 
the fundamental states of both the emitting and 

1N. Bohr and F. Kalckar, Kgl. Danske Videnskab. 


Selskab. 14, No. 10 (1937). 
2 J. Thibaud, Comptes rendus 191, 656 (1930). 
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J=0 but also the processes J=0 to J/=1, J=0 
to J=2, J=0 to J=3. In the latter processes 
the alpha-particle would carry away the angular 
momentum h, 2h, and 3h, respectively. According 
to Gamow’s*® formula the relative probabilities 
of these processes for a nuclear charge Z=82 
would be approximately as given in Table I. 
Instead of finding such closely spaced levels as 
given in Table I no fine structure has been 
observed in transformations in which there is 
reason to assume that the fundamental states of 
initial and final nuclei have J=0. 

(2) The explanation for the fine structure in 
the cases in which it is observed‘ is based on the 
idea that fine structure occurs if a transition to 
the fundamental state of the final nucleus would 
require the alpha-particle to carry away a very 
high angular momentum. Since such a process is 
improbable, the disintegration will lead in most 
cases to excited states which possess angular 
momenta differing from that of the original 


TABLE I. Relative probabilities of transitions calculated from 
Gamow’s formula for a nuclear charge Z=82. 





TRANSITION RELATIVE PROBABILITY ENERGY CHANGE 

J=0 J=0'| 1.00 0 

J=0 J=1| 0.84 3.3 Kev 

J=0 J=2/| 0.60 10 Kev 

J=0 J=3/ 0.36 20 Kev 
3G. Gamow, Structure of Atomic Nuclei (Clarendon 


Press, 1937), Eq. (56), p. 103. 
* Reference 3, p. 108. 
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nucleus by fewer units. These excited states 
frequently lie several hundred thousand volts, 
or even more than a million volts, above the 
fundamental state of the final nucleus. If the 
final nucleus had rotational states lying close to 
its fundamental level there would be no reason 
why the alpha-process should not lead to these 
low lying levels rather than to those observed. 

(3) The strong gamma-radiation accompany- 
ing some beta-disintegrations has been given an 
explanation similar to that of the fine structure 
of the alpha-rays.' Also in this case the hard 
gamma-rays would be difficult to understand if 
low lying energy levels were available. 

(4) Metastable states of nuclei which do not 
emit gamma-rays but rather eject electrons from 
the atomic shell (coefficient of internal con- 
version equal to unity) are found to be lying 
more than a million volts over the fundamental 
state. If low lying rotational levels exist these 
strong, discrete beta-rays should always be 
accompanied by at least equally strong gamma- 
rays of only slightly smaller energy leading from 
the metastable state to the low lying rotational 
states. 

(5) Metastable states of much longer life have 
been postulated by Weizsacker® in order to 
explain isomeric nuclei. These long lives he 
explained by assuming that all transitions from 
the metastable state to states of lower energy 
would be connected with a great change of 
angular momentum. Evidently if Weizsacker’s 
hypothesis is to be upheld, one must assume the 
existence of a special level with a high angular 
momentum lying lower than all but a few of the 
rotational states of the droplet. 

(6) The existence of the low lying levels of the 
droplet would not be consistent with the explana- 
tion given for the long periods of the potassium 
and rubidium radioactivities. According to 
Klemperer? a large change in the angular 
momentum should account for the small transi- 
tion probability. Klemperer’s suggestion has been 
confirmed for the case of rubidium by Kopfer- 
mann and Konopinski.* Low rotational levels of 
the product nucleus would make possible less 


5 Reference 3, p. 150. 

®C. F. v. Weizsicker, Naturwiss. 24, 813 (1936). 

7Q. Klemperer, Proc. Roy. Soc. A148, 638 (1935). 

§ Kopfermann, Physik. Zeits. 38, 960 (1937); Konopin- 
ski, N. Y. Meeting Am. Phys. Soc., Feb. 25, 1938. 
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strongly forbidden transitions. Still smaller tran- 
sition probabilities have been assumed to account 
for the stability of the neighboring isobars. The 
low lying levels would invalidate this explanation 
too. 

We conclude therefore that no such low lying 
levels are present in the nucleus. 

It is our purpose to show that the application 
of the exclusion principle to the constituents of 
nuclei gives the possibility of reconciling the 
droplet model with the absence of low lying 
levels such as would follow from the rotation of 
the droplet. As the simplest model we consider 
N alpha-particles on a circle and spaced at equal 
distances. The rotation of the ring into itself 
(i.e., around an axis passing through the center 
of the circle and perpendicular to the plane of 
the circle) will interchange the alpha-particles. 
Since the first exchange takes place after a 
rotation by the angle 27/N, the Bose statistics 
of the alpha-particles will exclude all rotational 
states except those having the angular mo- 
mentum J=0, J=N, J=2N, ::-. 

If we now permit rotations of a group of 
alpha-particles around any axis in space such 
simple exclusions of rotational states will be 
obtained only for nuclei for which particularly 
symmetrical configurations can be assumed. 
Thus, Wheeler® has found that for four alpha- 
particles arranged in a tetrahedron the first 
rotational excitation has J=3. If we arrange six 
alpha-particles on the vertices of an octahedron 
(this structure has been proposed by Wefelmeier'® 
for Mg**) the lowest excitation will lie at J=4. 
For twelve alpha-particles at the vertices of an 
icosahedron the lowest rotational excitation will 
be at J=6 (this icosahedron was assumed by 
Wefelmeier for Fe™* with the thirteenth alpha- 
particle in the center of the icosahedron and 
with two additional neutrons about the localiza- 
tion of which no definite statements are made). 
Apart from such symmetrical cases, J=1 or 
J=2 should be expected as the lowest allowed 
rotational state since, in most cases, there will 
be either no rotation which effects an interchange 
of the constituents of the nucleus, or else only a 
few rotations, mostly through the angle z, will 
effect such an interchange. 


9J. A. Wheeler, Phys. Rev. 52, 1089 (1937). 
10 W. Wefelmeier, Naturwiss. 25, 525 (1937). 
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We must remember, however, that the con- 
stituents of a nucleus cannot be represented as 
firmly bound to equilibrium positions. Indeed, 
the dimensions of the nucleus alone show that 
the kinetic energies of the nuclear constituents 
are larger than the nuclear binding energies, 
that is, the kinetic and potential energy in the 
nucleus are comparable and the binding energy 
is obtained as a relatively small difference 
between the two. Therefore it will be important 
to consider the following exchange between the 
constituents of the nucleus: we first perform a 
rotation which, although it does not actually 
produce an interchange of the particles, does 
bring them into positions which are close to 
positions obtainable by simple exchange; then 
we have to move the constituents of the nucleus 
only through comparatively short distances in 
order to complete the exchange. It will be 
necessary to investigate which of these two 
motions will correspond to a higher frequency: 
the rotation through a finite angle or the dis- 
placement through short distances. The first will 
be the case if the constituents are bound to 
equilibrium distances with potentials large com- 
pared with their kinetic energies and, if thus, no 
appreciable displacements can be obtained with- 
out passing over potential barriers. Then the 
nucleus may possess all rotational states which 
are allowed according to the symmetry of the 
equilibrium position. If, however, the rotation 
through a finite angle takes a longer time and 
corresponds to a smaller frequency than the 
subsequent displacements through short dis- 
tances, then the process of rotation and the 
rotational levels will be strongly disturbed and 
we may expect the nucleus to act like a body of 
higher symmetry than would correspond to any 
“equilibrium position’ of the constituent par- 
ticles. The effect may be that the predicted low 
rotational levels are absent. 

In the next section we shall consider the effect 
of exchanging the particles for simple models 
consisting of constituents satisfying Bose sta- 
tistics (alpha-particles). In the last section we 
shall discuss the angular momentum and the 
excitation energy which we may expect for the 
lowest levels in heavy nuclei according to this 
modified droplet model. 


J. A. WHEELER 


II. EXAMPLES 


In the following we select a few simple 
examples which illustrate how the possibility of 
small displacements from equilibrium increases 
the effective symmetry of a system and how this 
increase in symmetry reacts on the structure of 


the allowed rotational energy levels. 


Two particles on a circle 

Let us picture two identical particles, obeying 
Einstein-Bose statistics, free to rotate in a plane 
about the same fixed point to which both are 
coupled by inextensible rods of the same length. 
Between the two particles acts a potential which 
depends upon the angle @ between their two radii 
vectors. The potential has the same value for @ 
and 2x—@. The minimum occurs for a value of 
6= 6) close to, but definitely less than, 180°. In 
the extreme case where the potential valley is 
very deep and very narrow (rigid coupling) the 
energy levels of the system are simply those of 
a rotator with one degree of freedom and a 
moment of inertia J: 


E=£oth?m?/2]. (1) 


All values of m are allowed because the system 
has no symmetry with respect to its rotational 
axis. If, now, the rigidity of the binding is 
decreased, it must be remembered that the 
system actually possesses two equilibrium posi- 
tions, at 4 and 2r—, and that the potential 
barrier between them, while high, allows the 
possibility of resonance between the two equi- 
librium positions. Thus each of the levels con- 
sidered in the case of rigid coupling is really a 
double level with a separation 2AE proportional 
to the frequency with which penetration through 
the barrier occurs. One easily finds from the 
W-K-B approximation that 


AE= +hvejass(”) 
xX exp — f (21 V0) —Evn(ny i a6 h. (2) 


Velass(@) represents the classical frequency of 
vibration in the mth vibrational level, and the 
integral in the exponent depends upon the area 
under the potential barrier. The positive or 
negative sign in the equation has to be applied 
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Fic. 1. The diagram in the lower right corner represents 
the triangle-dumbbell model. The curves represent the 
energy levels for different cases of coupling. In case (d) 
the coupling is rigid and the levels correspond to different 
values of the rotational quantum number m. In case (c) 
the coupling is loosened and we see the effect of the dumb- 
bell penetrating the potential barrier to turn into another 
symmetrical position. In case (a) the dumbbell and tri- 
angle are allowed to rotate freely with the angular mo- 
menta bh and ah. The total angular momentum mh=ah+bh. 
In case (b) a small coupling is introduced causing repulsion 
between levels of the same m value. (See levels for m=4.) 


according as the vibrational functions in the 
two minima are superposed with the opposite 
sign or with the same sign. If rotations of the 
system are taken into account, the symmetry of 


the wave function with regard to the exchange 


of two alpha-particles requires that it should 
remain unchanged if a rotation through 180° is 
performed and the oscillator is moved from one 
minimum to the other. Since a half-revolution 
multiplies by (—1)” (m is the rotational quan- 
tum number), only the energy shift which 
corresponds to superposition of the vibrational 
functions with the same sign is allowed for 
m=(0, 2, ---, and with the opposite sign for 
m=1, 3, ---. Thus putting |AE| =Q, we obtain 
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for the total energy 
E= Eyi,(n) +h?m?/21—(—1)"Q. (3) 


As the rigidity of the binding is decreased, Q in 
general increases. When 20 becomes as great as 
(h?/21){mo?—(my—1)*}, already the odd level 
m=m,—1 is raised higher than the even level 
m=mo. This corresponds to the frequency of 
interchange 2Q/h becoming greater than the 
frequency of rotation fm/27J. Finally we shall 
find that a number of even m levels will lie 
below even the lowest odd level (m= 1). 

This approximates then the result we find for 
the case of equilibrium position of alpha- 
particles opposite each other (model for Be’), 
where only even values of J occur. 


Triangle-dumbbell model 

A second instructive example is provided by 
what we may call the triangle-dumbbell model, 
illustrated in the lower right-hand portion of 
Fig. 1. The two rigidly coupled particles in the 
dumbbell are identical, obey Einstein-Bose sta- 
tistics, and have only a single degree of freedom, 
that of rotation about their common center of 
gravity, which is also the center of rotation for 
the triangle. The latter is composed of three 
identical particles, again rigidly coupled and 
following the symmetric statistics. As extreme 
cases we may consider (a) no interaction between 
the two rotational degrees of freedom and (d) 
rigid coupling. In the first case the energy levels 
are determined by the sum of the quantum 
expressions for two free rotators of moments of 
inertia A and B: 


E=h'a?/2A +h*b?/2B. (4) 


Here the rotational quantum number } may take 
on only the values 0, +2, +4 and a the 
values 0, +3, +6 ---, because of the symmetry 
of the system. In case (d) the coupling is so 
strong that we obtain effectively the problem of 
a single rotator having no symmetry against 
rotation, and possessing a moment of inertia 
equal to the sum A+B8 of the two separate 
moments. The energy of the system is now 


E=E,+h?m?/2(A+B), (5) 


where m may take on all integral values. 
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b ! 


Fic. 2. (a) Bipyramid alpha-particle model for the 
nucleus Ne*®. (6, c), a distortion separating alpha-particles 
2,4 and bringing particles 5,1 into contact. The result is an 
exchange plus rotation through 90 degrees. 


The question as to how the transition occurs 
between the two preceding cases now presents 
itself in an interesting form, for it is clear that 
the order of the various rotational levels must 
go through some sort of mixing process. 


First, let us commence with the case of a very weak field, 
say 
V(8—a)=—W cos 6(B—a), 
which possesses a sixfold periodicity because of the 
symmetry of the two parts of the system. It is easy to 
calculate from the unperturbed wave functions 
WV = (21) exp i(aa+b~), (6) 


that the first-order perturbation energy vanishes. In the 
second order we have 


E=Wa?/2A+hb?/2B 


— (W?/36h?)(A7+B)— 4 1- yt (7) 
— 3A+3B ° 





which represents a repulsion between energy levels pos- 
sessing the same values of m=a+b. (See case (b) in Fig. 1.) 
On the other hand, proceeding from the side of strong 
interactions toward the intermediate case, we find it 
simplest to separate the wave equation 
(h?/2A)d?*/da?+ (h?/2B)d*y/dp?+{E—V(p—a)}y=0 (8) 
by writing the wave function as the product 
¥ =f(¢) exp (cmd) (9) 

of a part representing rotation of the system as a whole by 
the angle @= (Aa+B8)/(A +B) timesa part which depends 
only upon the relative orientation ¢=8—a of the triangle 
and dumbbell. The wave function of the internal torsional 
vibration, f(¢), must satisfy the equation 

(i /2)(A74 B>)d'f/de? + |E— Ve) |f=0, 


and it is also clear that it must obey a certain periodicity 


(10) 
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f(e+2x/6) =f(¢) exp —im(24A/3+22B/2)/(A+B). 
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condition. From the picture in the lower right-hand portion 
of Fig. 1, one sees that a rigid rotation of the whole system 
by an angle @=(27A/3+27B/2)/(A+B) followed by a 
rotation of the triangle relative to the dumbbell by an 
angle g=2z/6 will carry the system from its original 
position to one where the rod has been given a half-turn, 
the triangle a third of a revolution, The requirement that 
the wave function have the same value for the two equiva- 
lent configurations tells us at once (see Eq. (9)) that 
(11) 
The solutions of Eqs. (10) and (8), subject to this condition, 
may be obtained with the help of the W-K-B approxi- 
mation and give us the energy values 
E=Eyin(n)—(—1)"2hvetass(n) 
{cos 2rm(A/3+B/2)/(A+B)} 
Xexp {— S| p(¢)|dy/h} +h?m?/2(A +B), 


where m is the quantum number of the torsional vibration, 


(12) 


The figure shows the influence of the possi- 
bility of displacements of the triangle and 
dumbbell with respect to each other. In the case 
(d) of the rigid coupling the low symmetry of 
our model caused the level of the total angular 
momentum m=1 to be the first excited state. 
In the intermediate cases (c) and (0) the level 
m=1 moves up strongly whereas the energy of 
the levels m=2 and m=3 are not raised to an 
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Fic. 3. The energy levels at the left represent the low 
rotational states allowed for the bipyramid neon model 
of Fig. 2. Jh stands for total angular momentum and Kh 
is the component of this along the axis of figure of the rigid 
model. The possibility of exchange between different 
configurations of the same energy (see Fig. 2) perturbs the 
levels. The shift is shown as a function of the value of the 
exchange integral Q for which a reasonable estimate is 
0.17 mMU. Additional excited levels are found if the 
vibrations of the model are taken into consideration. 
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ROTATION OF THE 
equal extent. Eventually, as we approach case 
(a), both the energy of the lowest excitation and 
also its angular momentum (m=3) become 
considerably greater than in the case of rigid 
coupling. 


Bipyramid model 


A further example of the qualitative consider- 
ations discussed above is furnished by a problem 
which can be considered as a model for the 
nucleus Ne*®. We may picture the stable state 
of the system as that in which five alpha- 
particles are situated at the vertices of a double 
pyramid (see Fig. 2). If we assume that the 
“bond distances” are the same between all pairs 
of alpha-particles, then this configuration gives 
the maximum possible number of bonds (9). 
Diagrams (b) and (c) show that a distortion of 
the nucleus in which the 2—4 bond is stretched 
and the distance between 1 and 5 is shortened 
will carry the system over into an equivalent 
position of the same energy but different orienta- 
tion. For purposes of identification we may 
symbolize the original configuration by the 
expression 15 and the equivalent one by 42. 
By 15 we indicate that the alpha-particles 1 and 
5 are located at the two opposite poles of the 
bipyramid and the three alpha-particles at the 
base are numbered clockwise when viewed from 
1 towards 5. Thus 51 and 15 differ in the order 
of the alpha-particles in the base. Corresponding 
to the 20 ordered pairs, we have 20 different 
configurations of equal energy. Transitions simi- 
lar to the change 15 to 42 pictured in the figure 
are 15 to 23, and 15 to 34. One can see quali- 
tatively that other changes from 15, such as from 
15 to 51, will encounter higher potential barriers. 

A simple criterion for the most important 
“channels of intercommunication”’ is the value, 
r’, of the sum of the squares of the displacements 
of the individual alpha-particles which occur 
during the distortion. In terms of the distance, 
c, between neighboring alpha-particles, r* is 
easily computed to be 0.74c? for the change 
from 15 to 42, 2c? for the change to 51, 1.75c? for 
the change to 13, 2c? for the change to 31, and 
2.07c* for the change to 24. As the penetration 
factor goes down exponentially with 7’, it is 
clear that the most important type of exchange 
is the one represented in the figure. 
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Fic. 4. The lower portion is a schematic diagram for the 
variation of potential energy which occurs when the system 
is subjected to a distortion carrying it from one equilibrium 
position P to an equivalent one P’. The order of magnitude 
of the potential barrier between P and P’ is estimated to be 
1.5 mMU for the change from 15 to 24. The fictitious poten- 
tial barrier B has been introduced in such a way as to 
obtain a simple form for the zero-order wave function. The 
contour lines of this wave function are pictured above. 


Neglecting at first these exchange phenomena, 
we may divide the energy of the nucleus into 
the kinetic energy of rotation of a symmetric 
top and the zero-point energy of the 3X5—3—3 
=9 vibrational degrees of freedom. In terms of 
the moments of inertia C and A along and at 
right angles to the figure axis, and the total 
angular momentum Ji and its component Kh 
along the figure axis, we have 


E=Eot+J(J+1)h?/2A 
+(K*h?/2)(C-'—A-). (13) 


We shall not be interested in states of the 
system in which the vibrational degrees of 
freedom, 4:1, @2, -, are excited, and shall 
therefore treat the system as a rigid rotator. 

Only for certain values of J and K is it possible 
to obtain a rotational wave function whose value 
is unaltered by those rigid rotations of the model 
which effect an interchange of the alpha- 
particles, namely, 


K=0; even J: V~P (cos 0) fvin(q) ; 
K=3n0; any J: V~P;“(cos 6) 
x fexp («Kx)+(—1)/ exp (—7Kx) |fein(g). (14) 
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Here x measures rotation about the axis of 
figure of the double pyramid, and @ is the angle 
between this axis and some fixed direction in 
space. The left-hand portion of Fig. 3 shows the 
position of the low rotational levels allowed 
according to our analysis. The absolute values 
of the rotational energies have been obtained by 
using A = 4.7 X10-*5 g cm’, C=3.0X10-* ¢g cm’. 
The occurrence of exchange processes by 
penetration through potential barriers means 
now that the total function will be a 
superposition of functions of the type (14): 


wave 


20 
W=(20)-3 S fx6 (Oc, xe)fvin(Girer Gaer ***) (15) 
c=1 


Here the summation goes over all 20 stable 
configurations of the five alpha-particles. The 
Eulerian angles @, and x. and the vibrational 
coordinates g, are defined in the same way for 
the various configurations, so that WV is sym- 
metric with respect to all interchanges of the 
alpha-particles. If we write V. for the potential 
for which the cth term in VY would represent an 
accurate solution of the wave equation, and V 
for the actual potential, then by making a first- 
order perturbation calculation for the change in 
energy due to leakage through the barriers, we 
obtain 


E=Ey+E,a(J, K) 
+3 ffe*(ee, xaa)few*(Qaa)(V — Vis) 


Xf (O15, xis)fyin(Gis)dr. (16) 


In the derivation we have made use of the 
symmetry of W and the fact that from a given 
configuration, 15, the most important exchange 
processes lead only to the three equivalent 
equilibrium positions 42, 23, and 34. 

To carry out the integration over angles, we 
note from Fig. 2 that the orientation of the 
configuration 15 differs from that of 42 by a 
rotation, R, of 90°. Hence according to group 
theory, the integral of fx°* (O42, x42) fx (O15, x15) 
over all space orientations may be expressed in 
a simple manner in terms of the elements 


Dx, x(R) and D’_x, x(R) of the standard 
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rotational transformation matrix" for R=90°, 
whose values are given by 
D «x (0, 90°, 0) 

=(—1)4- KD _x, x(0, 90°, 0) =cux 


J-K (J+K)! 
Cyx=2-/ > (-1)/ muemeenms 
poe f\(J+K-—f)! 

(J—K)! 
xX ———-, (17) 
f\(J-K-—f)! 
After the integration over angles, we have left 
the vibrational part of the exchange integral, 
defined by 


Q _ J Freinlaud( ] — J is) fvin(Gis)dq. ( 1 8) 


Our result for the position of the low rotational 
levels, taking account of exchange, is expressed 
in the equation 


E=EotEkralJ, K) 


[CrK K=0; even J | 
+3 lo-| , . (19) 
| 20% K=3n+0; any J) 
The energy values are shown in Fig. 3 as a 


function of the magnitude of the vibrational 
exchange integral, Q. A reasonable estimate gives 
Q=-—0.17 mMU, so that the first excited 
rotational level, J=2, is raised up 0.75 mMU 
with respect to the ground state. However, 
already if Q be as great as —0.38 mMU, the 
effect is so pronounced that the levels J/=3 and 
J=2 become interchanged, with the level of 
greater angular momentum lying lower. The 
example shows in a striking manner the dis- 
turbance of the normal order of rotational levels 
caused by the existence of resonance between 
geometrically similar configurations of a system 
of five identical particles. 

To obtain the estimate we have just used for the value of 


the vibrational exchange integral Q, we express the poten- 
tial and kinetic energies for small oscillations from equi- 


librium in terms of suitable normal coordinates qi, go, «++, Qn? 
V = (ma/2)2wi*q;", (20) 
T = (ma/2)24;?. (21) 


NE. Wigner, Gruppentheorie und Quantenmechanik 


(Braunschweig, 1931), p. 180. 
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Here w;/2m is the classical frequency of vibration of the 

corresponding normal coordinate. The wave function for 

the normal state (zero-point vibration) has the value 


W = (ma/mh)”"'*(wiw.-+-w,)' exp |—LTmawigi/2h!, (22) 


if we suppose the vibration to occur in the neighborhood 
of one equilibrium position P. Actually the representative 
point of the system can penetrate through the potential 
barrier represented schematically in the lower portion of 
Fig. 4 and vibrate about the second minimum P’. If, as 
shown in the figure, we join together in the center of the 
barrier the unperturbed wave functions for the two 
minima, we obtain a solution of the wave equation 


(h?/2mg)2PWV/d9i?+(E— Vi)¥ =0 (23) 


for the unperturbed value of £. Here the potential V; is not 
the same as the actual potential V, however, but has been 
obtained from it by superposing a fictitious potential wall, 
very narrow and extremely high, in such a way as to bring 
about continuity in the normal derivatives of the wave 
functions on the two sides. Using gnorm to denote the 
coordinate measured parallel to the direction PP’, and 
calling the height and width of the barrier B and dq, we 
easily obtain in fact the relation 


— (h?/ma)(OWV tert /OQnorm) = Big¥ (24) 


by integrating (23). 

To obtain the correct energy value, we must remove the 
above-mentioned fictitious potential. Treating this alter- 
ation as a small change, we obtain by a first-order pertur- 
bation calculation 


AE=0= fv(V—Vi)dr=— f BaqvdS 
= (h?/ma) S Viett(OV iett/Anorm)dS, (25) 


where the integration goes over the m—1 dimensional 
surface normal to PP’. It would appear at first sight that 
the result should be divided by 2 because both parts of the 
total wave function are normalized to 1; but in the one- 
dimensional case one can readily prove that a compensating 
factor of 2 enters, just because of the fact that the pertur- 
bation is not small. It seems reasonable to expect this 
result in general, so we use (25) as it stands, without 
rigorous analysis. The integration may be carried out 
directly, from (22), and gives 


Q = —hw(maw/rh)*(r/2) exp | —(maw/h)(r/2)*}. (26) 


Here r is the distance from P to P’ and w/2zx is that 
frequency with which the system would oscillate if con- 
strained to move along the straight line PP’ with a 
potential energy which at every point P; on its path (see 
upper part of Fig. 4) is given by the value of V, not at P,, 
but at that point P2, equally distant from the surface of 
join, where the wave function has its maximum value. 
The contour representation of the wave function in the 
figure is drawn in such a way as to emphasize that in 
general no normal mode of vibration will tie in the direc- 
tion PP’. 

Following the methods used in reference 9, we find for the 
normal vibrations of neon characteristic energies ranging 
from hw, =0.8 mMU to hws = 2.2 mMU. An estimate for the 
value of the vibrational exchange integral is obtained from 
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Eq. (26) by substituting for w a harmonic mean of the 9 
characteristic frequencies, and for r? the value 0.74¢ 
1.2 mMU and r=7.8X 10748 
cm we find Q= —0.17 mMU, which is the value used above 


mentioned earlier. With hw 


in discussing the energy level perturbations shown in 
Fig. 3. It is of interest to note that the value of the vibra- 
tional exchange integral given by (26) will for no value of w 
ever exceed in magnitude the limit Q= —0.925h?/mar?. 


III. First Excirep STates OF HEAVY NUCLEI 


The actual behavior of nuclei cannot be repre- 
sented in any mathematically rigorous form; 
we must be content to estimate the order of 
magnitude of the excitation energies as they 
would follow from the simplest possible assump- 
tions. We must particular in 
what way the excitation energies will change if 


investigate in 


the number of the nuclear constituents increases. 

The reason why the first excited rotational 
levels will be strongly disturbed is that the 
phase of the wave function will be changed by 
a rotation. If a further small displacement of 
the alpha-particles suffices to complete an ex- 
change of the constituents (and restore the sign 
of the wave function) then it follows that the 
phase of the wave function must change on a 
comparatively short path. Such change of phase, 
however, requires an increase of kinetic energy 
unless a potential barrier forces the amplitude 
of the wave function to approach zero irrespec- 
tive of its phase. 


Alpha-particle model 
We shall first 
barriers are present. Let a nucleus composed of 


assume that no_ potential 
N. alpha-particles® be pictured in the 3N,- 
dimensional configuration space. A phase change 
y over a distance 7 will cause the momentum to 
increase by the amount ig/r, where r’ is the sum 
of the squares of the displacements of the 
alpha-particles. If originally the momentum and 
kinetic energy were zero, the result would be to 
increase the energy by h’?¢*/2mr’. Even if in the 
original case a considerable zero-point energy 
and corresponding momenta fp» are present, the 


2 The nucleus is assumed to be composed of alpha- 
particles rather than of protons and neutrons in order to 
simplify the argument. A detailed consideration of particles 
which possess spin and obey the Fermi statistics would 
lead in a more complicated way to results similar to those 
obtained in the text. 
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average increase of the kinetic energy is still 


{pothe/r}n?/2m— po?/2m=h?¢?/2mr?, (27) 
since the expression | Poig/r}s vanishes, because 
of the random direction of p in the unperturbed 
state. 

In order to obtain the perturbation of the 
rotational levels, the expression h?¢?/2mr? should 
be averaged over all orientations. For an order 
of magnitude estimate, it will be sufficient to 
substitute both for ¢ and for 7 their maximum 
possible values. Thus, we set 


g~T 


r~ N,(8/2)?, (28) 


where JN, again is the number of alpha-particles 
and 6 is their average distance. We thus obtain 
for the probable increase in excitation energy™ 


AE~4r*h?/2N.m.8, (29) 
It will be seen that this quantity decreases with 
the —1 power of the mass and reaches for heavy 
nuclei a value of about 100 kilovolts. Thus, 
according to this picture, we shall not expect in 
general an excitation energy of less than this 
order of magnitude for nuclei of any atomic 
number. Exceptions may occur if the nucleus 
has a spin due to the spin of protons or neutrons; 
then the lowest excitation energies might be 
determined by spin-orbit interaction rather than 
by the factors considered above. 

In the model just described, no definite 
statement can be made about the angular 
momenta of the levels. Of course, if the nucleus 
is composed of alpha-particles obeying Bose 
statistics, the wave function for the lowest level 
must not have any nodes and must thus belong 
to the quantum number J=0. The first excited 
level, then, may have any J value for which the 
rotational energy does not considerably exceed 
the AE given in (29). Therefore we have for the 
probable upper value for the rotational quantum 
number J,x. of the first excited state the relation 


RT exe(Jexe+1)/21~42°h?/2Names?. (30) 


Here J is the moment of inertia of the nucleus, 


13 This estimate of the first excitation assumes that of all 
the low lying levels that for J=0 is the only one which, 
due to the favorable phase relations, is not raised. 
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given by 


T=10-(18/22)*N 25/362. (31) 


Thus we obtain 


Jexe~2.2Ne', (32) 


which for heavy nuclei gives an angular quantum 
number J.x.=8. 

Another characteristic J value is the 
above which the order of rotational levels will 
become normal. We shall call this J value the 
quantum number of collective rotation, Joi, 
since for higher angular momenta the rotational 
energies approach those of a rigid body. We 
obtain for J..; the relation 


(h®/2T) | (Foor +1) (Teor +2) = Jeor(Feor +1) } 


~4r*h?/2N.m,6, 


one 


(33) 


which leads to 


Toi ~2.4Nal. (34) 


For heavy nuclei J..; is about 36, and the 
corresponding energy of rotation is E,..~2 Mev. 
Therefore one cannot expect to determine the 
moment of inertia of a nucleus described by our 
present model from a study of energy levels and 
their angular momenta, since in the region where 
the rotational levels assume their proper order 
there will be present a great number of vibra- 
tional excitations. 

Another interpretation for J. is obtained by 
comparing the classical frequency of rotation 
with the frequency of the small displacements 
which complete the rotational exchange of alpha- 
particles. This comparison gives 


AJ o1/2e7] ~AE/2th, (35) 


which is equivalent to Eq. (33). 

To obtain an interpretation of J.x-, we shall 
compare the time necessary to rotate an alpha- 
particle on the surface of the nucleus through 
half the distance between two alpha-particles, 
with the time #/AE required for an essential 
disturbance of the phase of the wave function 
by the displacement of the alpha-particles. We 
then have 


{r/(5/2)\hJexe/I~AE/h, (36) 


which leads to a value for J.x. differing from 
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that given in (32) only by the factor (7*/10)'. 
For J values smaller than J,x~ it is impossible 
to follow the rotation of the nucleus, even if we 
restrict ourselves to an angle 6/2r and if we use 
the best possible method of following the rota- 
tion, that is, a coordinate system fixed by a 
statistical study of all alpha-particles. 


Thomas-Fermi model 


In the above estimates, we used a model in 
which alpha-particles could exchange without 
crossing potential barriers. It is of interest to 
see how these results compare with estimates 
which one can obtain from the Thomas-Fermi 
model of the nucleus.'* We consider a nucleus 
which consists of N particles (protons and 
neutrons). The number of protons and neutrons 
shall be even and equal, so that each cell in 
phase space will be occupied by four particles. 
If the radius of the nucleus is r and the maximum 
momentum P,,, we have for the volume in phase 


space 
(37) 


For the maximum angular momentum of a 


(42/3) Pm3(4x/3)r3 = (N/4)(2xh)3. 


particle we obtain 


Pr /h=2.4(N/4)'=2.4Nol. (38) 


The close resemblance between this value and 
Jexe may be explained in the following way: if 
the neutrons and protons are supposed to move 
in spherically symmetrical fields and if their 
angular moments are composed according to the 
rules of vector addition, the lowest state will 
become an S level in consequence of the attrac- 
tion between constituent particles. If, then, one 
neutron or proton is lifted to a level of the same 
or a neighboring shell a group of excited levels 
will result the average angular momenta of 
which will be comparable with the value given 
in (38). Thus the J values for the lowest excited 
states will be similar for the a-particle picture 
and the Thomas-Fermi model. 

An estimate of the energy of the first excited 
state in the Thomas-Fermi model reveals no 
immediate agreement with the alpha-particle 
picture. If the excitation energy AE is due to 

4 Similarity of results may be expected in view of the 


comparison carried out by Euler, Zeits. f. Physik 105, 553 
(1937). 
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the transition from one shell into the next one 
we may write 


AE=P,,AP,,/mu, (39) 


where AP,,, the change of momentum of the 
excited neutron or proton, is given by 


(42/3)3Pn2APn(4x/3)r3=e(2rh)®. (40) 


g is the degree of degeneracy in the shell into 
which the neutron or proton is lifted. Dividing 
(40) by (37) we obtain 


3AP »/Pn=4e/N. (41) 


Since the density of nuclei is roughly constant 
P,, will not vary with N in a systematic way. 
However, g~2P,,7/h will be proportional (ac- 
cording to (38)) to N!. Thus we get 


8P,,*r &- 10’ ev 


AE~———_ =- 
3myNh Ni 


(42) 


This equation differs from (29) both in its 
dependence on N and in the high values of 
AE~2-10° ev which are obtained from it for 
heavy nuclei. 

The Thomas-Fermi model leads, however, to 
closer agreement with the alpha-particle picture 
if one carries out the discussion in greater detail. 
We shall discuss two such improved pictures. 

First we shall assume that we deal with an 
incomplete shell and that the first excited state 
corresponds to the transition of a neutron or 
proton from one degenerate state of this shell to 
another. Then the excitation energy will corre- 
spond to the diminished coupling of the particles 
in the incomplete shell. Its order of magnitude 
will be given 
between two particles in the same orbit. For 


by the change of interaction 


short range exchange forces this interaction will 
be proportional to N~.'® Since for the lightest 
nuclei the interaction is'® of the order of 10’ ev 
to 4-10’ ev for the heaviest nuclei excitation 
energies of 10° ev are obtained. This, as well as 
the dependence on JN is in agreement with (29). 

The above estimate was based on a model in 
which the nuclear proper function is not simply 

18 Wigner has found this dependence from the application 
of the symmetry properties of nuclear wave functions to 
the kinks in the mass defect curve. See Phys. Rev. 51, 954 


(1937). 
16 Breit, Condon and Present, Phys. Rev. 50, 825 (1936). 
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a product function of the individual particles 
but rather a linear combination of those products 
which would have the energy if each 
particle is assumed to move in a spherically 
symmetrical field. A second kind of improvement 
on the individual particle approximation would 
be to retain the product form of the proper 
function but to drop the assumption that the 
particles move in a field of spherical symmetry. 
Then in estimating the excitation energy AE the 
equations (39), (40) and (41) remain valid. 
Only g must be set equal to one since the absence 
of symmetry in the field will remove the de- 
generacy of the individual particle proper func- 
tion. Thus we obtain 


3AP,../ Po an 4 ‘N 
AE=4P,,2/3myN = (2.7-107/N) ev, 


same 


(41a) 


and (43) 


which has the same dependence on N and the 
same order of magnitude as (29). 


Comparison with experience 

Unfortunately there is no direct experimental 
procedure to determine the position of the first 
excited level of a heavy nucleus. The arguments 
adduced in the introduction against low lying 
rotational levels might be used however to test 
the orders of magnitude we have obtained in 
Eqs. (29) and (43) for the excitation energy. 
We shall discuss here only our last argument (6) 
concerning the stability of isobars and the 
radioactivity of potassium and rubidium. 

The upper limits of the beta-ray spectra of 
potassium and rubidium are 7-10° ev and 
2.5-10° ev, respectively. The explanation of the 
long periods require that there should be no 
levels in the product nuclei lying much lower 
than the beta-limit (unless beta-transitions into 
these levels are strongly forbidden). Thus we see 
that the lowest excitation energy of suitable 
angular momentum must be a few hundred 
thousand volts in both cases; the limit is more 
strict for the heavier element, rubidium. 

The number of neighboring pairs of isobars 
can be explained statistically if one postulates a 
mechanism whereby such a pair will be stable 
whenever the energy difference between initial 
and final atoms is less than 50,000 ev.'’ Now 


17 Bethe and Bacher, Rev. Mod. Phys. 8, 199 (1936). 
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with a neutrino of zero mass a beta-disintegration 
of this energy will yield too short a period of 
decay if the transition is an allowed one. A 
satisfactory explanation of the isobars can be 
obtained, however, by assuming that the angular 
momenta of the two isobars differ by three or 
more units. If the difference in angular momenta 
is great enough the energy difference might be 
10° ev or even more. This explanation will only 
hold, however, if in general there exist in the 
product nucleus no excited levels below 10° ev 
with suitable angular momenta to which transi- 
tions may occur. We see therefore that general 
experience points towards a value of ~10° ey 
for the lowest excitation energy. The calculations 
on the ‘‘thermal properties” of nuclei indicate'’ 
that the lowest excited states cannot be much 
higher than this quantity. 


The effect of potential barriers 


In conclusion we shall discuss on the basis of 
the alpha-particle model how the dependence of 
the lowest nuclear excitation on NV, will change 
if potential barriers must be crossed when the 
alpha-particles are displaced. For the lightest 
nuclei the transition from freely moving alpha- 
particles to strongly bound alpha-particles is a 
gradual one. As has been shown in the examples 
discussed in the second section the low lying 
energy levels approximate gradually the levels 
of a rigid top as the potential barriers become 
higher. 

For heavy nuclei, however, even low potential 
barriers which for light nuclei had no strong 
effect will become important. The reason is that 
in heavy nuclei many alpha-particles have to 
cross the barriers which causes an increase both 
in the effective width and height of the barrier. 
Thus for heavy nuclei low lying rotational levels 
can be avoided only if the displacement of the 
alpha-particles can be performed on such a path 
in configuration space on which the “kinetic 
energy” of the system never becomes negative. 
It is sufficient to have low potential barriers to 
effect a sudden change to the rotational spectrum 
of a rigid body. This change may be called a 
change of phase (or state) since it occurs suddenly 
for a system of many particles. 


Landau, 


18 Weisskopf, Phys. Rev. 52, 295 (1937); 
Physik. Zeits. Sowjetunion 11, 556 (1937). 
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It is interesting to note in this connection that 
the low temperature modification of liquid 
helium!’ has an unmeasurably small viscosity.”° 
Since helium does not solidify unless it is com- 
pressed it is reasonable to assume that it is kept 
in the liquid state because of the zero-point 
vibrations of the atoms*' which, for helium, are 

19 W. H. Keesom and M. Wolfke, ‘Comm. Phys. Lab. 
Leyden,’’ No. 190b. 

20P. Kapitza, Nature 141, 74 (1938); J. F. Allen and 
A. D. Misener, Nature 141, 75 (1938). 


21Simon, Nature 133, 460 and 529 (1934); F. London, 
Proc. Roy. Soc. 153, 576 (1935). 
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comparable with the heat of evaporation. The 
experiments on the viscosity of helium indicate 
that in a “quantum liquid” the particles may 
rearrange without crossing potential barriers. 

The absence of low lying rotational levels in 
heavy nuclei definitely indicates therefore that 
if nuclei are to be compared with a phase of 
matter in macroscopic experiments the correct 
analog to use is not a crystallite but a droplet 
of a “quantum liquid” such as the low tempera- 
ture modification of liquid helium. 
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New Experimental Evidence for the Existence of a Neutrino 


H. R. CRANE AND J. HALPERN 
University of Michigan, Ann Arbor, Michigan 
(Received March 24, 1938) 


A new method is used for determining the energy of recoil of the nucleus in the individual 
beta-disintegration process. A gaseous compound of radiochlorine is placed in a cloud chamber. 
The clearing field is removed long enough before the expansion to allow the ions to spread out 
so that the resulting droplets can be seen individually. A cluster of droplets appears at the 
beginning of the track, and this is believed to be produced by the recoil nucleus. From the 
number of droplets an estimate is made of the kinetic energy of the nucleus, and this is com- 
pared with that calculated from the observed curvature of the beta-ray track. It is found that 
the laws of mqmentum as well as those of energy indicate that a third particle participates in 


the disintegration. 


INTRODUCTION 


N the experiment to be described' the mo- 

menta of both the electron and the recoil 
nucleus have been simultaneously measured for 
the elementary beta-disintegration process. This 
is the first experiment which has given any 
information at all regarding the momentum 
relations in the individual disintegration event. 
Although the results are of limited accuracy, 
they strongly indicate that momentum is not 
conserved between the electron and the nucleus 
alone. Hence the laws of momentum, as well as 
those of energy, indicate that a third particle 
participates in the disintegration. 

The idea of observing the disintegration of a 
substance in the form of a gas in a cloud chamber 
has been suggested many times as a possible 
way of measuring the momentum or energy of 


1 Halpern and Crane, Abstract 5, New York Meeting 
of the American Physical Society, Feb. 25-26, 1938. 





recoil and the direction of recoil of the nucleus 
emitting the beta-ray. The difficulty with such 
a scheme is that the length of track made by 
the nucleus is far too short for observation, even 
in a cloud chamber operated at the lowest 
possible pressure. The nucleus will, however, 
produce a number of ion pairs concentrated in 
a very small region of space, and the number of 
ion pairs will be a function of the kinetic energy 
of the nucleus. It occurred to the authors that 
if these ions could be allowed to diffuse into a 
cluster several millimeters in diameter before 
the condensation were brought about, the indi- 
vidual droplets could be counted, and hence the 
kinetic energy of the nucleus could be estimated. 
It was found that this could be accomplished, 
and by applying a magnetic field to the chamber 
it was possible to know the momentum of the 
electron in each case, so that the estimated 
momentum of the nucleus could be compared 
with that of the electron. 
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APPARATUS AND EXPERIMENTAL METHOD 


An automatic Wilson cloud chamber? 15 cm 
in diameter and 4 cm deep was used, which was 
filled with air at atmospheric pressure, ethyl 
alcohol vapor, and a trace of a gaseous compound 
of the radioactive material. In order to allow 
the ions to diffuse before the alcohol vapor 
was made to condense upon them, the clearing 
field was removed a fraction of a second before 
the expansion. This was accomplished by means 
of an adjustable contact in the automatic timing 
system. Two Sept 35 mm cameras, having f :3.5 
lenses were placed 25 cm above the chamber 
and at a 30 degree angle to each other. This 
gave close-up views from two directions of only 
a part of the chamber. The tracks were illumi- 
nated by means of a carbon arc, flashed at about 
200 amperes, and were photographed at 1/15 
second exposure. A rectangular coordinate sys- 
tem of lines 3 cm apart was scratched on the 
lower surface of the glass top plate of the 
chamber. By measuring the apparent position 
of the beginning point of a track with respect 
to the lines on the plate in each of the views, an 
accurate determination of the distance of the 
point below the plate was possible. The two 
stereoscopic views were projected at the same 
time, but were not superimposed on the same 
screen. The position of a given point on a track 
relative to the coordinate system was simply 
recorded for each of the two views. Only those 
tracks which were found to begin at a point 
well within the illuminated portion of the 
chamber were included in the final data. 

The radioactive material used in the first part 
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Fic. 1. Numbers of droplets produced by the recoil 


nuclei (corrected), shown in relation to theoretical pre- 
dictions made on three different assumptions, as follows: 
curve I, no neutrino; curve II, neutrino and electron 
escape in the same direction; curve III, neutrino and 
electron escape in opposite directions. 


* Crane, Rev. Sci. Inst. 8, 440 (1937). 
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of the experiment was radiochlorine, Cl**, pro- 
duced by bombardment of NaCl by 6 Mev 
deuterons in the University of Michigan cyclo- 
tron. This emits negative electrons with a half- 
life of approximately 37 minutes, and with an 
upper energy limit of about 5 Mev.* The combi- 
nation of high energy, long half-life and fairly 
low atomic weight make this element suitable 
for our purpose. A gaseous compound of the 
chlorine suitable for introduction into the cloud 
chamber was made by the following procedure. 
The bombarded NaCl (a fraction of a gram) was 
put into a test tube containing about 5 cc of 
sulphuric acid, diluted 1 to 3, and a few grams 
of manganese dioxide. This produced free chlo- 
rine in solution. A few cc of ethylene dichloride 
were introduced into the test tube. This formed 
a layer above the sulphuric acid. Ethylene gas 
(generated by dropping ethylene dibromide upon 
zinc metal) was then bubbled through the liquid 
in the test tube. The ethylene combined with 
the chlorine in the lower layer, forming ethylene 
dichloride, which dissolved in the ethylene di- 
chloride of the upper layer. After bubbling for 
about 15 minutes the upper layer was drawn off 
and was shaken vigorously with a large quantity 
of water in order to remove traces of sodium 
(which was also radioactive) and also traces of 
free chlorine and of sulphuric acid. The amount 
of radioactive ethylene dichloride introduced 
into the cloud chamber was just that which 
could exist in the vapor state at room tempera- 
ture, which was several drops. Into the otherwise 
dry chamber just slightly more ethyl alcohol 
than necessary to saturate the volume was 
added, and the chamber was sealed up with air 
in it at atmospheric pressure. The presence of 
an excess of either of the liquids was avoided in 
order to insure that as large a fraction as possible 
of the radioactive compound would remain in 
the vapor state. The time taken for the chemical 
procedure and for clearing the cloud chamber 
was 30 to 40 minutes. Needless to say, only a 
minute fraction of the original activity of the 
NaCl appeared in the cloud chamber at the end 
of the procedure, but sufficient activity was 
always obtained if the NaCl was bombarded for 
about 5 minutes with 6 microamperes of deu- 


3Kurie, Richardson and Paxton, Phys. Rev. 49, 368 


(1936). 
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Fic. 2. Four different beta-ray tracks of intermediate energy, showing 


clusters of droplets at their beginning points. 


terons at 6 Mev. The time of bombardment 
affords a convenient control on the final activity 
in the chamber, if the chemical procedure is 
standardized. From 200 to 300 photographs were 
taken with each loading of the chamber. 

The gaseous compound of radiophosphorus*® 
(P®*), which was used as the control in the 
second part of the experiment, was prepared as 


follows. Calcium phosphide (Ca;P2) was bom- 
barded by 6 Mev deuterons in the cyclotron. 
This was decomposed by dropping water onto 
the solid Ca3P2 in an atmosphere of nitrogen. 
The resulting mixture of PH; and nitrogen was 
passed into the cloud chamber, which had 
already been flushed with nitrogen. The neces- 
sary amount of absolute ethyl alcohol was added, 
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Fic. 3. An electron of high energy, in which the mo- 
mentum of the electron is sufficient to account for the 
observed recoil of the nucleus. A stereoscopic pair of 
photographs is shown. 


and additional nitrogen flushed through the 
chamber to reduce the activity to the optimum 
value. It was necessary to take precautions 
against having traces of oxygen present, because 
of its action on PH3. 

In determining the number of droplets to be 
assigned to the recoil nucleus, we made the 
simplifying assumption that the ions originating 
at a given point on a track diffuse in such a 
way as to fill with uniform density a sphere 
whose center is at the source and whose radius 
depends upon the time and speed of diffusion. 
On this assumption, the ions for which the recoil 
nucleus is responsible will fill a sphere of diameter 
equal to the width of the electron track. In the 
two-dimensional photograph a circle of this 
diameter is drawn around the nuclear cluster, 
the center of which is supposedly at the point of 
disintegration. The number of droplets inside 
the circle which are not due to the nucleus but 
to the electron can easily be calculated on the 
above assumptions to be nearly pR, where p is 
the linear density of droplets along the elec- 
tron track, and R is the radius of the circle. 
Accordingly we have subtracted pR from the 
observed number of droplets in each cluster. p 
was determined by counting the droplets along 
several centimeters of the particular track. It is 
obvious that this correction, as it is here de- 
termined, is not dependent upon the track lying 
in a plane perpendicular to the direction of the 
camera: it is correct for any projection, provided 
p is determined in the same projection. 
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Fic. 4. Two separate tracks of low energy which are 
clear cases in which the nuclear recoil is larger than that 
expected from the electron momentum alone. 


RESULTS 


The results of the droplet counts on the 17 
cases obtained with radiochlorine are shown as 
the dots in Fig. 1. No cases of electrons having 
than 1 Mev energy are included, for a 
reason which will be given later in the discussion. 
The number of droplets produced by the electron 
within the region over which the droplets due 
to the nucleus extend, is subject to statistical 
fluctuation. To give an idea as to the amount 
of uncertainty due to this cause, a_ vertical 
line of length equal to (pR)! is drawn through 
each point. It is more difficult to estimate the 
uncertainty in the actual counting of the drop- 
lets, because this depends largely upon the 
quality of a given track. We estimate that this 
error ranges from 2 to 5 droplets. 

In Figs. 2 to 4 are shown a number of the 
tracks obtained. Although all tracks were ana- 
lyzed stereoscopically, both views are reproduced 
in only two of the examples. The ruled lines on 
the chamber top are visible, and are 3 cm apart. 

The theoretical expectations for the number 
of ions produced by the nucleus were derived on 
three different assumptions. 1. That the mo- 
mentum of the nucleus is equal to that of the 
electron (equivalent to assuming that there is no 
neutrino). 2. That a neutrino of zero or very 
small rest mass escapes in the same direction as 
the direction of flight of the electron, and that 
its momentum is (W—£)/c, where W is the 
upper limit energy of the beta-ray spectrum of 
Cl*8 (taken as 5 Mev) and where £ is the energy 
of the electron, calculated from its curvature in 
the magnetic field. 3. The same as 2 except that 
the electron and neutrino escape in opposite 
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directions. In all the above calculations the 
kinetic energy of the nucleus was calculated 
from the momentum (its mass is known) and it 
was assumed that one ion pair (or two droplets) 
would be produced for each 30 electron volts of 
kinetic energy. The three possibilities are plotted 
in Fig. 1. It is clear that if all possible directions 
of emission of the neutrino with respect to that 
of the electron are considered, curves 2 and 3 
will represent upper and lower limits for the 
number of droplets, and all the area between 
them will be permitted. 

It is well known that a moving ion is able to 
produce other ions when its kinetic energy (in 
the moving coordinates of the collision) is only 
slightly above the ionization potential of the gas 
through which it is passing. But the question 
of the total number of ion pairs a heavy ion of 
a few hundred electron volts energy will produce 
before it comes to rest is difficult to answer. 
It would not be expected to produce more than 
the usual one ion pair per 30 ev, and could not 
under any circumstance produce more than one 
per 16.7 ev, which is the average ionization 
potential of air. On the other hand it could lose 
rather large fractions of its energy by thermal 
collisions without ionization. The numbers of 
droplets indicated by the points in Fig. 1 may 
therefore in some cases give too low values for 
the recoil momentum of the nucleus, but are 
less likely to give values which are too high. 

Possibilities other than ionization by the recoil 
nucleus exist for producing ions in the neighbor- 
hood of the disintegrating atom. The escaping 
beta-ray may disturb the outer electronic struc- 
ture of the atom, which may in its readjustment 
emit an electron of a few hundred volts energy, 
or may emit a soft x-ray quantum which will be 
absorbed in the immediate vicinity. To test this 
we carried out a control experiment, with 
radiophosphorus, P*®, instead of Cl**. P® is a 
good control because it has about the same 
atomic weight as Cl**, emits negative electrons, 
but has a much lower upper energy limit. The 
greatest possible recoil energy for the P* is 85 
electrons volts, which should produce less than 
6 droplets. Out of 25P* disintegrations in the 
gas obtained, only two had detectable clusters 
of droplets at their beginning points. This 
negative result makes us feel certain that the 


clusters of droplets observed in the experiment 
with Cl** were actually caused by the recoil of 
the nucleus. 


DISCUSSION 


Meaning of the results obtained 


It can be seen from the data presented that 
in the cases in which beta-rays of small mo- 
mentum are emitted from radiochlorine the 
momentum received by the nucleus is usually 
much greater than that which it could have 
received from the beta-ray alone. This statement 
is not subject to any assumption about the 
efficiency of ionization by the nucleus: it is true 
even if we assume that the entire kinetic energy 
of the nucleus is utilized in ionization. This 
result has a bearing upon certain of the hy- 
potheses which can be made as to the mechanism 
of beta-disintegration. 

1. Suppose that the nucleus retains the excess 
energy temporarily, emitting it later by any 
means we may care to postulate. Momentum 
would then be conserved between the beta-ray 
and the nucleus, and the data would have to 
fall on curve 1, Fig. 1. 

2. Suppose that the beta-ray carries off all the 
energy of disintegration, partly as kinetic energy 
and partly as extra mass, as has been suggested 
by the group in this laboratory,’ by Jauncey,® 
and undoubtedly by others. Since the curvature 
of the beta-ray track in the magnetic field gives 
the momentum directly (assuming constant 
charge) the data in this case also would have to 
fall on curve 1, Fig. 1. 

3. If the beta-ray spectrum owes its continu- 
ous distribution to a failure of the law of the 
conservation of energy, there is no reason to 
believe that the law of the conservation of 
momentum does not fail also. The present data 
are not able to eliminate this possibility. How- 
ever, the experiment is inherently capable of 
throwing light upon the question if carried far 
enough, because if the above were true we should 
expect to find some cases in which the recoil 
momentum of the nucleus would be less than 


4 Breit, Rev. Sci. Inst. 8, 141 (1937); Zahn and Spees, 
Phys. Rev. 53, 524 (1937). 

5 Jauncey, Phys. Rev. 52, 1256 (1937); Phys. Rev. 53, 
106 (1938); Phys. Rev. 53, 197 (1938); Phys. Rev. 53, 
319 (1938). 
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the minimum which could be given to it by any 
combination of electron and neutrino energies 
or angles. Such cases would fall in the area 
below curve 3 on the right-hand side of the 
diagram in Fig. 1. 

Of all the possibilities mentioned, the hy- 
pothesis that a neutrino is emitted simultane- 
ously with the electron seems to fit the data best. 


Lower energy component in the Cl** beta-ray 
spectrum 

Kurie, Richardson and Paxton® investigated 
the beta-ray spectrum of chlorine after bombard- 
ment with deuterons and found that it could be 
resolved into two components of about equal 
intensity and of K. U. upper limits 6.1 and 1.5 
Mev. The actual observed upper limits, which 
we prefer to use in connection with our experi- 
ment, are about 5 and 1.2 Mev. They also found 
that gamma-rays of 2.4 and 1.9 Mev were 
emitted. The most reasonable assumption seems 
to be that the 1.2 Mev beta-ray spectrum and 
the gamma-rays are associated, since the sum of 
the energies, 1.2+2.4+1.9=5.5 Mev, is about 
the same as the energy of the other beta-ray 
component. It is therefore reasonable to say 
that no gamma-rays are associated with the 
5 Mev beta-ray component. 

In our experiment we must be careful to 
include only those disintegrations which belong 
to the 5 Mev component, if we are to compare 
the results with the predicted recoil momenta. 
An electron which belongs to the 1.2 Mev 
spectrum, emitted practically simultaneously 
with the two gamma-ray quanta might give 
considerable recoil to the nucleus, but would 
tell nothing about the neutrino. We have there- 
fore included in our data only those cases in 
which the electron has more than 1 Mev energy; 
it is fairly certain that these belong to the 5 Mev 


spectrum. 


Use of other gases and radioactive substances 


The use of hydrogen in the cloud chamber 
instead of air would decrease the linear density 
of ionization of the electron track, which would 
be an advantage, in that it would decrease the 
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number of droplets to be subtracted from the 
nuclear cluster. On the other hand the droplets 
would fall faster, make sharp photography more 
difficult, and, most important of all, the number 
of ions produced in hydrogen by the nucleus 
would probably be much smaller than in air. 
In the collision of the recoil nucleus (A** after 
the disintegration) with He, the kinetic energy 
of the A** in the moving coordinates of the 
center of gravity of the system would be very 
small indeed, and ionization might not result. 
It is therefore better to use a gas of atomic 
weight comparable to that of the Cl**. 

The use of He® as the radioactive material 
would give a much greater recoil energy, because 
the upper limit of its spectrum is about the 
same as that of Cl** and its mass is much smaller. 
He® has a half-life of 1 second, which means 
that it would have to be introduced into the 
chamber before each expansion. It would be 
difficult to maintain ideal conditions in the 
chamber under these circumstances, but perhaps 
not impossible. 


Direction of emission of neutrino 

The present experiment is inherently capable 
of giving the distribution in angle of emission 
of the neutrino with respect to the direction of 
emission of the electron. An accurate knowledge 
of the recoil momentum and the electron mo- 
mentum would make it possible to solve for the 
angle in each case. Before the results can be 
interpreted in this way we must obtain exact 
information as to the efficiency of ionization by 
slow ions. The higher recoil energy of He® would 
lessen this difficulty. 
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Cross sections for the combined processes of absorption+ inelastic scattering of fast neutrons 


(E 24.5 Mev and E27 Mev) have been measured for a number of elements by the use of 


absorbing cylinders and aluminum or iron detectors. A smooth increase with increasing atomic 


number is found. For neutrons of energies 27 Mev the ratio of these cross sections to the cross 


sections for absorption+total scattering (Dunning) is nearly constant at a value of about 0.4. 


The existence of a scattering of fast neutrons in which little or no energy is lost has been 


demonstrated experimentally and approximate cross sections for the process have been ob- 


tained. Cross sections for an “‘effective absorption”’ of medium fast neutrons have also been 


determined. These cross sections, which are smaller than those obtained for the faster neutrons 


mentioned above, show a maximum at about Z=50 and a marked decrease for the heaviest 


elements. The cross sections obtained by Aoki and Lea from experiments on the yield of gamma- 


rays from the action of fast neutrons on matter are in agreement with the present results. 


INTRODUCTION 


HE absorption scattering cross sections of 

atomic nuclei for fast neutrons have been 
measured for a considerable number of elements 
by Dunning.' These cross sections are the sum 
of the cross sections for three types of interaction 
of fast neutrons with nuclei: (a) absorption 
(transmutations due to simple capture, proton 
emission, double neutron emission, or alpha- 
emission), (b) elastic scattering, and (c) inelastic 
scattering. Experiment has not yet shown de- 
cisively what fraction of the total cross section is 
due to each of the three named processes. In the 
present investigation two series of experiments 
have been performed which yield cross sections 
of atomic nuclei for the combined processes of 
absorption and inelastic scattering for fast 
neutrons having energies greater than 4 and 7 
Mev, respectively. As Weisskopf has pointed 
out,? the results of experiments on the inelastic 
scattering of such fast neutrons are especially 
suitable for theoretical investigation because the 
number of excited states in which the nucleus 
may be left when the neutron is reemitted will 
be large enough to make statistical considerations 
valid. Effective absorption cross sections for 
medium fast neutrons (so called in order to 
distinguish them from the faster neutrons used 
in the first-mentioned experiments) were re- 


1 Dunning, Phys. Rev. 45, 586 (1934). 
2 Weisskopf, Phys. Rev. 52, 295 (1937). 


“I 


ported in a previous paper.* These measurements 
have been repeated and extended to include a 
greater number of elements. Interesting com- 
parisons can now be made between the results 
obtained with neutrons of different energies, and 
the conclusions thus reached are of assistance in 
interpreting the measurements of Aoki‘ and Lea‘ 
on the intensities of the gamma-rays excited by 
fast neutrons in various substances. 


SCATTERING OF FAST NEUTRONS 


The considerations of Bohr®: 7 have shown that 
the phenomena which occur upon the impact of 
a neutron (or other heavy particle) with an 
atomic nucleus must be treated as two inde- 
pendent processes, namely the formation and 
eventual breaking up of an intermediate meta- 
stable compound nucleus. Frenkel made _ the 
suggestion that the compound system formed by 
the impact of a fast neutron with a nucleus may 
reemit a neutron with an energy corresponding 
to the ‘‘temperature’’ of the residual nucleus, 
in a manner analogous to an ordinary evapora- 
tion process.* This process of inelastic scattering 
was treated in a quantitative manner by Weiss- 
kopf? who applied statistical methods to the case 


3 Seaborg, Gibson, Grahame, Phys. Rev. 52, 408 (1937). 

# Aoki, Proc. Phys.-Math. Soc. Japan 19, 369 (1935). 

5 Lea, Proc. Roy. Soc. (London) A150, 637 (1935). 

6 Bohr, Science 86, 161 (1937). 

7 Bohr and Kalckar, Vidensk. Selsk. Math.-fys. Medd. 
14, 10 (1937). 

’ Frenkel, Physik. Zeits. Sowjetunion 9, 533 (1936) 
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of collisions of fast neutrons (E>3 Mev) with 
heavy nuclei (A>100). He found that the 
reemitted neutrons should have, on the average, 
energies much less than the energies of the 
incident neutrons. For example, neutrons with 
an energy of 10 Mev will be reemitted with a 
Maxwellian distribution having a mean energy 
of about 2 Mev. The nucleus is left in an excited 
state as a result of this energy loss by the 
neutrons and subsequently returns to the ground 
state by the emission of one or more gamma-rays. 

A study of the cross sections for absorption and 
inelastic scattering of this kind for fast neutrons 
may be made by using the radioactivity induced 
in aluminum or iron. These substances become 
radioactive under the action of fast neutrons 
through the reactions :°: !° 


AP??+n—-Mg"? +p Meg??—-Al*?+e- (10.3 min.) 
Fe+n—-Mn*+p Mn *®—->Fe*®+e- (2.6 hr.) 


In each case the final stable nucleus is identical 
with the initial nucleus. From the upper energy 
limits of the emitted electrons!'!—“ and the heights 
of the potential barriers for protons of the 
nuclei formed,!®> one may estimate the least 
energy required for the neutron to produce an 
observable radioactivity. From such considera- 
tions one estimates that neutron energies of 
about 4.5 and 7.0 Mev are required to activate 
aluminum (10.3-min. period) and iron, respec- 
tively. 

The effective absorption of such fast neutrons 
was measured by observing the decrease in 
activity of an aluminum or iron detector caused 
by surrounding a source of neutrons (200 mg 
Ra+Be) with a cylinder of any one of a number 
of elements. If the inelastically scattered neu- 
trons have a Maxwellian distribution of energies, 
the fraction having energies greater than 7.0 Mev 
probably will not exceed a few percent. Hence 


® Fermi, Amaldi, d’Agostino, Rasetti, Segré, Proc. Roy. 
Soc. (London) A146, 483 (1934). 

10The small amounts of Na*™ and Al?* produced in 
aluminum under the experimental conditions adopted 
have been taken into consideration. Iron detectors have 
the advantage that Mn* is the only radioactive element 
known to be formed by the action of neutrons on iron. 

11 Henderson, Phys. Rev. 48, 855 (1935). 

2 Brown and Mitchell, Phys. Rev. 50, 593 (1936). 

13 Gaerttner, Turin, Crane, Phys. Rev. 49, 793 (1936). 

4 Alichanow, Alichanian, Dzelepow, Nature 136, 257 
(1935). 

1’ Bethe, Rev. Mod. Phys. 9, 172 (1937). See Table 
XXXIII. 
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we may state that fast neutrons which have 
been inelastically scattered by the Bohr process 
will not, in general, retain enough energy to 
induce radioactivity in an iron detector. That 
this is also approximately true for an aluminum 
detector may be inferred from the fact that the 
results obtained with iron and aluminum de- 
tectors were found to be nearly the same. The 
success of the method probably also depends, in 
part, upon the fact that even those neutrons 
which retain sufficient energy to induce radio- 
activity in the detector after an inelastic scatter- 
ing process must be expected to have a reduced 
probability for doing so as a result of their 
diminished energy. 

The quantity measured in these experiments is 
the total cross section for all processes by which 
the neutron loses a large fraction of its energy or 
is absorbed. The values called absorption +in- 
elastic scattering cross sections in the text are 
therefore more strictly cross sections for absorp- 
tion+inelastic scattering with large energy 
losses. If a form of inelastic scattering occurs in 
which only a small fraction of the energy is lost, 
the experiments will not detect the process to 
any marked extent. 

The cylindrical absorbers used in these experi- 
ments were 17.8 cm long and approximately 10 
cm in diameter. The source of neutrons (13 cm 
long and 1.78 cm in diameter) was placed 
centrally in a 1.78 cm hole through the axis of 
the cylinder. Flexible targets of iron or aluminum 
(made of small strips held together with Cello- 
phane tape) fitted around the curved surface of 
the cylinders. The radioactivity induced in these 
targets was measured with a thin-walled steel 
Geiger-Miiller counter. In the measurements 
made without the absorbers, the geometrical 
conditions were accurately preserved by the use 
of identical blanks. When aluminum targets were 
used, measurements were made alternately with 
and without the absorbing cylinders in order to 
eliminate the effects of any changes in the 
sensitivity of the counter. When the longer 
(2.6-hr.) period induced in iron was used, the 
sensitivity of the counter was checked regularly 
during the measurements. 

The results of these measurements are pre- 
sented in Table I and plotted in Fig. 1 together 
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SCATTERING 


with Dunning’s cross sections! for the total 
scattering -+absorption of fast neutrons. 

The cross sections ¢ in Table I have been 
calculated from the fraction transmitted 7 by 
use of the relation 7=e-*%*', where WN is the 
number of nuclei per cm® and ¢ is the average 
path length of the neutrons in the absorbers. 
The probable errors given refer to the errors in 
the measurements only and are not intended to 
include systematic errors inherent in the method. 
Under the experimental conditions here used the 
average length of path of a neutron originating 
on the axis of the cylinder and passing in a 
straight line through the absorber and detector 
is 24 percent greater than the wall thickness 
stated in column 2 of Table I. (This 24 percent 
is a weighted average which takes account of 
the obliquity of the path of the neutron in both 
the absorber and the detector.) 

Two corrections to the cross sections in Table I 
must be considered. Neutrons which are scattered 
in the absorber with negligible energy losses will 
have, on the average, longer paths through the 
absorber than the calculations take into account. 
An approximate calculation of this correction 
may be made by assuming isotropic scattering 
and neglecting scattering processes after the 
first. Dunning’s scattering cross sections have 
been used for Al, Zn, Sn and Pb, and inter- 
polated values have been taken for Sb and Bi. 
It is assumed from Table II that 40 percent of 
the collisions are inelastic with large energy 
losses. (The calculation is not very sensitive to 


TABLE I. Measurements of absorption+inelastic scattering 
cross sections for fast neutrons. 














| 
| | FRACTION TRANSMITTED 
| 








WALL 
TuHIckK- | Fe*(n,p)Mn* AF?(n,p) Mg?’ Cross 
NESS DETECTION DETECTION SECTION 
ELEMENT cm E2>7 Mev E>4.5 Mev x 10% cm? 
Al | 4.2 | 0.7540.03 | | 9341.3 
Zn 4.2 0.59+0.03 | | 15.4+1.5 
Sb |} 4,2 0.73 +0.03 | 17.942.3 
Pb 4.2 | 0.69+0.03 | 21.342.4 
Bi | 4.2 0.71+0.02 | | 23.2+2.0 
Cc | 4.3 0.8340.02 | 4.2+0.5 
Al | 4.2 | | 0.784+0.02 | 7.8+0.8 
Zn |} 4.2 | 0.58+0.02 | 15.941.0 
Sn | 4.1 0.72+0.02 17.1+41.5 
Sb | 4.2 | 0.74+0.02 | 17.6+1.5 
Ph =| 4.2 | | 0.7140.02 | 19.7+1.6 
Bi | 4.2 | | 0.7440.01 | 19.8+0.9 
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Fic. 1. Upper curve, total scattering+ absorption cross 
sections for fast neutrons (Dunning). Lower curves, 
inelastic scattering+absorption cross sections for fast 
neutrons as measured with iron detectors @ or with 
aluminum detectors O. 


this quantity.) The result of such a calculation 
indicates that the cross sections in Table | 
should be reduced by about 13 percent (15 per- 
cent for Zn) to account for the effect of single 
scattering in the absorber. For scattering proc- 
esses after the first the calculation becomes very 
complex, but a careful consideration of the 
problem suggests that a complete correction for 
single and multiple scattering would reduce the 
calculated cross sections by about 20 percent 
(23 percent for Zn). 

The second correction to be considered is that 
which arises from the increased obliquity in the 
detector of those neutrons which have been 
scattered in the absorber with negligible energy 
losses. The magnitude of this correction may be 
estimated from a calculation similar to the first 
by which it is found that the average path 
through the detector is 6 percent greater than it 
would be in the absence of the absorber. On 
account of the exponential nature of the calcu- 
lation of the cross sections, this correction will 
amount to 16 percent of the total cross section 
(13 percent for Zn) and acts to increase the cross 
sections over the values given in Table I. Since 
these corrections act in opposite directions, the 
net result of both corrections is a correction 
smaller than the probable error of the measure- 
ments themselves, except possibly in the case of 
zinc.!>* 

16 A corrected value for the cross section of zinc with 
iron detection gives ¢=14.1X10-%* cm*. This value 


corresponds to a sticking probability of 0.43 for Table II 
instead of the value 0.47 there reported. 
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Ra+Be 
Fic. 2. Experimental arrangement for demonstrating 
scattering of fast neutrons with small energy losses 


(elastic or inelastic). 


Collie and Griffiths'®’ made measurements of 
the effective absorption of fast neutrons in 
spherical absorbers of several elements, and 
their results are in fair agreement with com- 
parable values in Table I. Ehrenberg!’ measured 
the effective absorption of fast neutrons in a 
silver cylinder using aluminum and silicon de- 
tectors. His results yielded a cross section for 
silver somewhat larger than that expected from 
the curves in Fig. 1. However, Ehrenberg used 
the wall thickness of his silver cylinder for the 
path length of the neutrons in the silver. Since 
the average path length of the neutrons was 
presumably greater than the wall thickness, it is 
reasonable that his result should have been too 
large. The magnitude of the correction involved 
cannot be estimated from his published data. 

Figure 1 shows a gradual increase in the cross 
sections for absorption + inelastic scattering from 
the lightest to the heaviest elements. No large 
differences were observed between the results 
obtained with aluminum and with iron detectors, 
although the results do suggest slightly larger 
cross sections for the neutrons detected by iron 
(E=7 Mev) than for the neutrons detected by 
aluminum (E=4.5 Mev), at least in the case of 

_ the heavier elements. The cross section given for 
carbon is doubtless too large because here 
elastic scattering may cause a neutron to lose 
nearly one-third of its energy in a single collision. 

The cross sections listed in Table I give the 
sum of the cross sections for absorption and for 
inelastic scattering with large energy losses. 
Theoretical considerations lead one to expect 
that for any but the lightest elements the true 
absorption cross sections for fast neutrons will be 
small compared to the cross sections for inelastic 


16 Collie and Griffiths, Proc. Roy. Soc. (London) A155, 


434 (1936). 
17 Ehrenberg, Nature 136, 870 (1935). 
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scattering.'* Thus one may say to a good approxi- 
mation that the cross sections given in Table | 
refer to the single process of inelastic scattering 
with large energy losses. 

One infers that the differences between 
cross sections for inelastic scattering +absorption 
as given in Table I and Dunning’s cross sections 
for total scattering+absorption give at once the 
cross sections for elastic scattering or for in- 
elastic scattering with small energy losses. This 
procedure assumes that the absolute magnitudes 
of both sets of cross sections are correct. Dun- 
ning’s cross sections were determined under good 
geometrical conditions, and none of the dis- 
turbing factors would tend to make his cross 
sections too large in any event. The values in 
Table I cannot be sufficiently in error to alter 
the conclusion that a large difference exists 
between the two sets of values. Thus it appears 
that the difference between these results and 
those of Dunning is real and attributable to 
elastic scattering or to inelastic scattering with 


the 


small energy losses. 

It is to be noted that 
scattering of fast neutrons may occur in which 
there is a less complete amalgamation of the 
incident neutron with the nucleus and which 
might be designated by the term surface scatter- 
ing. Such a process might often give rise to small 


a form of inelastic 


energy losses. 

Since no direct evidence for scattering with 
small energy losses (either elastic or inelastic) has 
been reported, experiments were carried out to 
give qualitative proof that such scattering does 
occur. A view of the apparatus used is given in 
Fig. 2. Two iron targets (composed of thin 
strips stacked together) were activated simul- 
taneously in the positions shown, and the in- 
tensity of their activities compared. If a neutron 
may be scattered and still retain sufficient energy 
to induce radioactivity in iron, then target 
No. 1 (Fig. 2) will be less active than target 
No. 2 because in the latter case the mass of iron 
not directly between the source and the target 
may scatter into the target neutrons that would 
otherwise be lost. The activity observed in 
target No. 2 was 25+5 percent greater than 
that observed in target No. 1. This result is in 
qualitative agreement with the supposition that 
~ 18 Bethe, Rev. Mod. Phys. 9, §65 (1937). 
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SCATTERING 


the difference between Dunning’s curve and the 
lower curves in Fig. 1 gives the cross section for 
scattering of fast neutrons with small energy 
losses (elastic or inelastic). 

Additional evidence for such scattering was 
obtained as follows: An aluminum detector in 
the form of a cylinder surrounded the source of 
neutrons, and measurements were made with and 
without a large iron block (back-scatterer) 
entirely surrounding this assembly. When the 
iron back-scatterer was in place, the 10.3-min. 
activity in aluminum was increased 4.5+1.0 
percent by the action of neutrons scattered from 
the iron block into the detector. 

Table II shows, in column 2, cross sections for 
scattering of fast neutrons with small energy 
losses (elastic or inelastic) obtained by difference 
from Dunning’s values for absorption-+total 
scattering and the values in Table I for absorp- 
tion+inelastic scattering as measured with iron 
targets. Column 3 gives the ratio of the ab- 
sorption+inelastic scattering cross section to 
Dunning’s cross section for the same elements. 
This ratio, which is nearly constant at a value of 
about 0.4, represents the fraction of the neutron- 
nucleus collisions which are inelastic with large 
energy losses. This ratio is the sticking proba- 
bility, £, as defined by Bethe'® in the equation 
t=0/7R*, where o is the absorption+inelastic 
scattering cross section and where, in this case, 
Dunning’s cross sections are used for 7R?°. 


SCATTERING OF MEpDIUM FAsT NEUTRONS 


In a previous paper*® cross sections for the 
combined processes of absorption and inelastic 
scattering were reported for medium fast neu- 
trons; these measurements have now been re- 
peated under improved experimental conditions 
and have been extended to include a greater 
number of elements. In these experiments a 
TABLE II. Cross sections for scattering with small energy 


losses and sticking probabilities for fast neutrons 
(E=7 Mev). 


STICKING 
Cross SECTION PROBABILITY 
ELEMENT < 10% cm? (SEE TEXT) 
Al 15 0.39 
Zn 18 0.47 
Sb 26 0.41 
Pb 37 0.37 
Bi 39 0.39 
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hic. 3. Experimental arrangement for measuring effective 
absorption of medium fast neutrons. 


source of fast neutrons was enclosed within a 
lead block of such dimensions that the neutrons 
to be studied passed through 30 centimeters of 
lead before reaching the detector. The detector 
used was a Geiger-Miiller counter surrounded by 
copper. Neutrons which impinge upon the copper 
nuclei may undergo inelastic scattering, and the 
gamma-rays resulting from this process are de- 
tected by the counter. In order to be detected in 
this manner, the neutron must have an energy 
not less than the energy of the first excited state 
of the copper nucleus. Since this method of 
detection certainly detects neutrons with energies 
less than a million volts, the average energy of 
the neutrons to which it responds will be much 
less than that of the aluminum and iron detectors 
used in the previous section. Such neutrons will 
be called medium fast neutrons in order to dis- 
tinguish them from the faster neutrons employed 
in the experiments described above. 

The exact experimental arrangement used was 
a modification of that described in the previous 
paper and is shown in Fig. 3. The source of fast 
neutrons (200 mg Ra+Be) was mounted at the 
center of the cavity A within the lead block, 
which was large enough to absorb nearly all of 
the gamma-rays emitted by the radium source in 
the direction of the counter. The counter used 
was made from a steel tube (12.5 cm X3.8 cm, 
0.13 mm wall thickness) and was filled with a 
mixture of argon and ethyl alcohol vapor (10 cm 
A+0.3 cm C,H;OH).!* Boron carbide and Pyrex 
glass surrounded the copper detector, as shown, 
to prevent slow neutrons from the source, or 
neutrons slowed in the floor and walls of the 
room from acting on the detector system. The 
entire assembly was mounted on a steel table well 


1’ Trost, Zeits. f. Physik 105, 399 (1937). 
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removed from neighboring walls or other hydro- 
genous material. 

When the radium-beryllium source was intro- 
duced at the cavity A (Fig. 3), the counting 
rate increased from the natural background of 
35 per minute to 160 per minute. This increase 
in counting rate was due chiefly to the gamma- 
rays produced by fast neutrons in the detector.* 
A small fraction of the increase caused by the 
radium-beryllium source will, of course, be due 
to the incomplete absorption in the large lead 
block of the primary gamma-rays from radium 
and its decay products. The magnitude of this 
effect was determined by methods previously 
described,* but another and more accurate deter- 
mination has been made by replacing the 
radium-beryllium source by 150 millicuries of 
pure radon.”° With the arrangement shown in 
Fig. 3, the effects due to primary gamma-rays 
were found to be very small; however, in the 
work which follows, the appropriate corrections 
have been made in every case. 


Measurements of effective absorption cross 
sections for medium fast neutrons. 


TABLE III. 


























FRACTION TRANSMITTED 
WALL Bee — 
THICK Cross 
| NESS Cu(n, ny)Cu | Pb(n, ny)Pb}| SecTION 

ELEMEN1 cm | DETECTION | DETECTION < 10% cm? 
C(p=1.7) | 4.3 | 0.959+0.015| 1.1+0.4 
Al | 4.2 | 0.933-40.015| 2740.6 
S(p=2.0) | 4.0 | 0.918+0.015) 5.71.1 
CI(CCh) 4.3 | 0.944+0.015| 5.141.5 
Fe 4.1 | 0.817+0.015| 5.7+0.5 
Cu 3.5 | 0.823+0.015) 6.5+0.6 
Zn | 4.2 | 0.83340.015| 6.7+0.6 
Cd(p=7.7) | 4.1 | 0.855+0.015| 9.2+1.0 
Sn 4.1 | 0.886+0.015) 7841.1 
Sb | 4.2 | 0.864+0.015| 10.5+1.2 
I(p =4.5) 4.1 | 0.922+0.015) 9.24+1.7 
Hg 3.9 | 0.851+40.015| 10.0+1.1 
Pb 4.2 | 0.881+0.010) 9.1+0.8 
Bi 4.2 | 0.927+0.010 6.3+0.9 
C(p =1.3) 4.1 | 0.9840.02 0.6+0.6 
. 4.2 Rawr 3.4+0.8 
S(p = 2.0) 4.0 0.93+40.02} 5.0+1.4 
Cl(CCh) 4.0 | 0.92+0.02) 8.342.1 
Fe 4.0 | 0.80+0.02) 6.2+0.7 
Cu 3.5 | 0.7740.02| 8.9+0.9 
Cd(p=7.3) | 4.3 | 0.82+0.02| 12.14+1.4 
Hg 3.9 | 0.86+0.02} 9.441.5 
Pb 4.2 | 0.90+0.02) 7.341.6 











20 The radon was enclosed in soft glass. Such a gamma- 
ray source gives a small intensity of neutrons from the 
action of alpha-particles on the constituents of glass. 
Their presence was easily demonstrated by the use of a 
BF;-filled proportional counter and their effect was 
allowed for in this determination. 


AnD? GG. F. 


SEABORG 


The effective absorption cross sections for 
medium fast neutrons were measured by sur- 
rounding the source at A _ with cylindrical 
absorbers and observing the decrease in counting 
rate of the Geiger-Miiller counter. Those neu- 
trons which have been inelastically scattered in 
the absorber with large energy losses will have a 
greatly diminished probability of exciting gamma- 
rays in the detector by a second inelastic 
collision. Experimentally, this will manifest itself 
as an effective absorption of the neutrons. The 
results of such measurements are presented in 
Table III together with results obtained in a 
similar manner by using lead instead of copper 
as the detector. Figs. 4a and 4b show both sets 
of data graphically. The cross sections given 
have been calculated from the mean free paths 
determined by the transmissions observed. The 
suitability of this method of calculation has been 
checked by showing that the absorption is 
approximately exponential with thickness of 
absorber. The values given are somewhat smaller 
than those reported in the previous paper be- 
cause of a more accurate determination of the 
intensity of primary gamma-rays reaching the 
counter from the radium source. 

The curves in Figs. 4a and 4b show a gradual 
rise with atomic number up to about Z=50 and 
a small decrease in cross section for the heaviest 
elements studied. In particular, the cross section 
for bismuth (Fig. 4a) is certainly much smaller 
than those for lead and mercury. A marked 
decrease for lead and bismuth has been observed 
by Aoki in experiments on the intensity of the 
gamma-rays excited in various elements by 
2.4 Mev (D+D) neutrons.‘ The values obtained 
by Aoki give a direct measure of the relative 
cross sections for inelastic scattering, provided 
that the number of gamma-rays emitted by the 
participating nuclei and the efficiency of their 
detection remain the same from element to 
element. Aoki’s results for the elements here 
studied are presented for comparison in Fig. 4c. 
The general agreement between the three sets of 
values shows that the small cross sections ob- 
served in the heavy elements must be real. 
These low cross sections for heavy nuclei indicate 
that medium fast neutrons either lack the 
energy necessary to undergo inelastic scattering 
with a high probability or, more likely, undergo 
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SCATTERING 


inelastic scattering with much smaller energy 
losses. The experiments already described with 
more energetic neutrons (Table 1) do not show 
this decrease in the cross sections for the heavy 
elements. 

The absolute values of the cross sections in 
Table III are very much smaller than those 
obtained in Table I for more energetic neutrons. 
This result is undoubtedly due in large part to 
the fact that neutrons which are inelastically 
scattered in the absorbers may still retain enough 
energy to be detected by the copper or lead 
detector. It may also indicate that medium fast 
neutrons which are inelastically scattered lose on 
the average a smaller fraction of their energy 
than the faster neutrons employed in the first 
experiments, in agreement with the calculations 
of Weisskopf who showed that the energy of the 
inelastically scattered neutrons should be ap- 
proximately proportional to the square root of 
the energy of the incident neutrons.? In addition 
it is not improbable that the sticking probability 
of fast neutrons increases with increasing energy 
of the neutrons as a result of the more closely 
spaced levels associated with higher energies of 
excitation in the compound nucleus. 

Measurements on the intensity of gamma-rays 
produced by the action of fast neutrons on 
matter have been made by Lea® as well as by 
Aoki. The results of the two workers differed 
markedly for the heavier elements. Lea found a 
continual increase in the cross sections with 
atomic number, even up to lead and bismuth, 
whereas Aoki found a falling off for these 
elements as already shown (Fig. 4c). Moreover 
the yield of gamma-rays seemed to be larger in 
the experiments of Lea than in those by Aoki. 
These differences are evidently attributable to 
the higher energies of the neutrons used by Lea. 
He used fast neutrons from the interaction of 
polonium alpha-particles with beryllium (with- 
out the interposition of lead or other substances 
which would have slowed down the neutrons by 
inelastic scattering), whereas Aoki used the less 
energetic (2.4 Mev) neutrons from the action of 
deuterons on deuterium. 

Several investigations have dealt with the 
intensity of gamma-rays produced by radium- 
beryllium neutrons acting on lead.*: *!~*3 In such 
experiments it has always been necessary to pass 
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Fic. 4. a, Effective absorption cross sections for medium 
fast neutrons as measured with a Cu detector; }, same, 
measured with a Pb detector; c, relative cross sections for 
the production of gamma-rays by 2.4 Mev neutrons 


(Aoki). 


the neutrons through large thicknesses of some 
dense substance in order to absorb the primary 
gamma-rays, and hence the neutrons emerge 
with reduced energies, as already pointed out. 
This fact is adequate to explain fully the fact 
that Lea, with his more energetic neutrons, 
observed large yields of gamma-rays from lead 
whereas Kikuchi and co-workers failed in their 
earlier experiments to find any gamma-rays 
whatsoever from the process. 

We take this opportunity to express our 
gratitude to Professor G. E. Gibson and Dr. W. 
F. Libby of the department of chemistry and to 
Professor J. R. Oppenheimer of the department 
of physics for valuable advice and suggestions. 
Thanks are due also to Dr. L. R. Taussig of the 
University of California Medical School for a 
gift of radon used in the gamma-ray control 
experiments. 


*t Kikuchi, Aoki, Husimi, Proc. Phys.-Math. Soc. Japan 
18, 297 (1936). 

* Gibson, Seaborg, Grahame, Phys. Rev. 51, 370 (1937). 

*3 Takeda, Proc. Phys.-Math. Soc. Japan 19, 835 (1937) 
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Honl used his dispersion theory to develop a formula for 
the x-ray scattering power of an atom for wave-lengths 
near the absorption edge. He expressed his formula as the 
difference Af between his value and f calculated by some 
method such as that of Hartree. In this paper it is shown 
that the F values obtained in the powder crystal experi- 
ments of Brindley, Wollan and Harvey, and Miller may be 
used to verify this formula when they are corrected by the 
proper temperature factor given by recent experiments in 
diffuse scattering. The diffuse scattering experiment also 
gives f values which may be used to verify the Hénl theory. 


1. INTRODUCTION 


HE study of the scattering of x-rays by 

matter has been considerably delayed by 
the difficulty of separating the contributing 
factors in the intensity measurements of the 
scattered x-rays. In the ideal case of the Thomson 
scattering from a group of free electrons, the 
scattered intensity at an angle ¢ with the primary 
beam is independent of the variable (sin }@)/X. 
The grouping of these same electrons into the 
tightly bound formation of the atom introduces 
a function of (sin 3¢)/ by virtue of the fact that 
there are definite space relations among the elec- 
trons. However, the complications due to these 
relations are not nearly so troublesome as the 
difficulties caused by the natural vibrations of 
the electrons, for the phase angle of a scattered 
wave from an electron is a function of the dif- 
ference between the frequency »v of the primary 
radiation and the natural frequency v, of the 
electron in the atom. On the classical picture, if 
vis much greater than p,, the scattered radiation 
is practically in phase with the primary radi- 
ation. If v is slightly less than »,, the scattered 
ray is approximately 180° out of phase with the 
primary ray. Between these limits any value of 
the phase angle is possible. Thus all scattering 
experiments which will be free from the influence 
of this phase change must be made with a 
primary radiation whose frequency is much 
larger than any of the natural frequencies of the 


It is shown that all existing data on the f values of zinc are 
in good agreement when they have been corrected by the 
Jauncey and Bruce temperature factor and the Hénl Af. 
A working formula for Hénl’s Af for any element and for 
any wave-length is given together with a table of theoret- 
ical values of Af for zinc. From a consideration of all experi- 
mental results for zinc a set of f values has been con- 
structed and these are given in a table. This table shows 
that there is electron distortion and reasons are advanced 
for believing that the distortion of the electron atmosphere 
extends deeper into the atom than the valence electrons. 


electrons. Practically, this means y>vx, where 
vx is the natural frequency of the K electrons. 

Hartree! and other wave mechanists have cal- 
culated the scattering factor of an atom for 
v>vx. They have obtained values which give 
better agreement with experiment than the 
simple classical theories. Hénl*? has developed a 
very complete quantum-mechanical treatment of 
the phase angle problem for v near vx. He has 
made use of his theory to predict the structure 
factors of Fe, Mo, and W. H6nl shows that his 
numerical values agree fairly well with certain 
experimental results. 

In addition to the above-mentioned scattering 
problems from free atoms, we encounter other 
factors when we go to the crystalline state. As 
the atoms forming crystals are crowded together. 
the form of the potential barrier* of the atom is 
changed. The outermost conduction electrons are 
free to move in a periodic field. The effect of the 
conduction electrons has been discussed by 
Zener* and Scharwiachter,*® but both theory and 
experiment indicate that any effect on the f 
values is negligible at values of (sin }@)/\ greater 
than about 0.2. The motions of the next layer of 
electrons are less profoundly affected by the 
interatomic spacing, but they do move in orbits 
which are probably affected by the anisotropy of 


1D. R. Hartree, Proc. Camb. Phil. Soc. 24, 89 (1928). 
2H. Hoénl, Ann. d. Physik 18, 625 (1933). 

3]. C. Slater, Rev. Mod. Phys. 6, 209 (1934). 

*C. Zener, Phys. Rev. 48, 573 (1935). 

5 W. Scharwichter, Physik. Zeits. 38, 165 (1937). 
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SCATTERING 


the crystal. We now have a picture of the crystal 
composed of a relatively few conduction elec- 
trons wandering to and fro in a rather free 
manner and of atoms whose scattering power has 
been made anisotropic. 

We must now modify the whole picture by 
considering the thermal vibrations of the atoms. 
The effect of temperature on the scattering 
power of the atoms in a crystal has had 
several theoretical treatments® for the case of 
isotropic crystals. Zener’ has made a theoretical 
study of the thermal vibrations in anisotropic 
crystals. For the cases of both isotropic and 
anisotropic crystals, it is assumed that each atom 
vibrates as a whole in the crystal lattice. A truer 
picture may be that the amplitude of vibration is 
different for the inner electrons than for the 
outer electrons—some distortion of the outer 
electron atmosphere being perhaps due to the 
thermal vibrations. However. in this paper we 
shall separate the thermal vibrations from the 
distortion of the outer electron atmosphere. 


2. Tue ATOMIC STRUCTURE FACTOR 


Hoénl’s method of attacking the problem of 
the scattering power of an atom in the region of 
anamolous dispersion was first suggested by 
Waller,* who used a formula involving induced 
dipoles, quadrupoles, and octopoles in an atom 
which was under the influence of an electromag- 
netic wave. Hénl succeeded in developing the 
solution into a rapidly converging series by the 
use of spherical harmonics. The fact that Hénl 
did not include in his calculations the effect of 
spin or relativity is of no consequence in this 
paper since hvx for zinc is less than mc? by a 
factor of 50. The structure factor obtained by 
Honl differs from the Hartree f by 


1 ma 
Af= —§,® —- —— —+quadrupole terms, (1) 
2 f+E 
where, according to Hénl’s notation, 
£, = (27e-*/9)(402’—Q;’), (2) 


x af ae 1 
a (ee ae 
(1—6)*7, \¢e—-1 ¢+41 


6 See C. Zener and G. E. M. Jauncey, Phys. Rev. 49, 17 
(1936). 

7C. Zener, Phys. Rev. 49, 122 (1936). 

51. Waller, Zeits. f. Physik 51, 213 (1928). 
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Fic. 1. The solid curve is the approximate Af for zinc. 
It is assumed that other terms in (1) than the first are 
insignificant. The dotted curve is the classical theory. 


s—s)*—hyx/Rhe 
pal de = (4) 


x=X/Xk, (5) 


and s is the screening constant for the K shell. 

A table of values of 6 is given elsewhere.’ These 

values indicate a slowly varying function of the 

atomic number with the value of 6 for zinc being 

0.205. 
On integration, (3) gives 

x 1 | 

QO.’ = ———— log |1 -—-- (3a) 
(1—6)? x? | 


|1—x| 
and Q;'= (2+ log |— )/ 1—6)*. (3b) 
\1+x 


If only the £,“ term is considered, (1) becomes 


4x? | 1 | 2x? 
Af=0.26; ———- log 
(1—6)? 


ee ee 


| x*| (1—6)' 


x |1—2| 
a log —it. (la) 
(1—6)® = [14x] 


This formula can be used for calculating Af for 
any atom and any wave-length and should be of 
use to other physicists when they are comparing 
experimental results with theory. 

The solid curve of Fig. 1 is the graph of & 
vs. x for zinc. In this simplification of Af the 
wave mechanical picture differs from the clas- 


*H. Honl, Zeits. f. Physik 84, 1 (1933). 
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sical picture principally because of the fact that 
wave mechanics considers the scattering power 
of the K electrons to be a function of the atomic 
number. For instance, in zinc the effective 
strength of the K shell is 1.31° rather than the 
classical value of 2. The broken curve shows the 
relation between £,“ and x for zinc on the 
Kramers-Kallmann-Mark theory.!° 
The second and third terms in (1) 
functions of (sin }¢)/X. In the second term 


27e-4 4x? x3 
a 5 ae, 
9 (1—6)*? (1-—6)8 


and in reality depends upon the absorption coef- 
ficient of the scattering material. This means that 
while it may have an appreciable value on the 
short wave-length side of the absorption edge, 
its magnitude on the long wave-length side is 
insignificant. This explains why the second term 
has been omitted in comparison with experi- 
mental data. In the experiments on zinc the 
value of the second term is two percent of the 
first term. 

The third term is a quadrupole term which 
would normally be significant in the Z absorption 
edge problem. Its effect is very small in the case 
of the K region and for the element zinc. 

Table I gives the values of £; and of 9; for 
various values of x in the case of zinc. 

There seems to be little data to test the angle- 
dependence part of Hénl’s theory. Rusterholtz! 
gives an empirical formula which he obtained by 
using primary rays of Zn Ka and a scatterer of 
copper. This made x =\/Ax = 1.038. His formula 
came out to be 


involve 


(6) 


Af =3.68—0.03(sin }¢)/A. 


The H6nl formula for these conditions is 


BRUCE 








0.025 2 cos* @ 
Af =3.50— +0.016— 
f-—3.50 1+cos* ¢ 
Numerically these two formulas agree pretty 


well over the experimental range. It is noted 
that the quadrupole term is included. In this 
case the quadrupole term is of some importance 
because of the closeness of x to unity. 

FACTORS 


3. POWDERED CRYSTAL STRUCTURE 


Structure factors for powdered zinc have been 
obtained for Cu Ka-rays (A= 1.54A) by Brindley® 
and by Wollan and Harvey" and for Mo Ka-rays 
(A=0.71A) by Miller.“* These experiments give F 
values rather than f values. The true atomic 
structure factors, f, are obtained from the experi- 
mental F values by means of 


f=Fe™. (7) 
In each experiment the front face method was 
used and the powder was obtained by sub- 
limation.: 5 

If we call fo.7; the structure factor of zinc for 
A=0.71A and fi.54 the structure factor for 


\=1.54A and if we define Af’ by 


Af’ = fo.ri— fis, 
Af’ =Afo.21—Afi.s4, (9) 
where Af, is obtained from (1) or (la) for 
x=/dx. For a given reflection e” is the same 
regardless of the wave-length and hence by (7) 
and (9) 


(8) 


then 


Af’ =(Fo.a— Fi.ssae™. (10) 
Values of F are listed in Table II. The values 
shown in column 4 are the averages of Brindley’s 
and of Wollan and Harvey’s values for \= 1.54. 


Column 5 shows Miller’s values for \=0.71. 




















2G. W. Brindley, Phil. Mag. 21, 790 (1936). 
18 E. O. Wollan and G. G. Harvey, Phys. Rev. 51, 1054 
10See A. H. Compton and S. K. Allison, X-Rays in (1937). 
Theory and Experiment, p. 292. 4R. D. Miller, Phys. Rev. 51, 959 (1937). 
1 A. A. Rusterholtz, Zeits. f.,Physik 82, 538 (1933). 1 W. A. Bruce, Rev. Sci. Inst. 8, 451 (1937). 
TABLE I. Values Fs £, and m for various values of x in the case of 3 Zinc, 
x 0.2 0.4 | 0.6 | 0.8 | 09 | 44 | 12 | 14 | 16 | 1.8 2.0 
aes Se Se a aS cumin és Sa ea 
H | 0.164! 0.301} 0.125 | —0.692 | —1.72 — 2.66 2.17 | —1.78 | 1.61 | -1.53 | —1.46 
m® | —0.206 | —0.726 | —1.52 | -2.49 | —3.01 | | 
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1 2 3 4 | 5 | 
| | | 

LINE ¥ eM Fi.ss Fo.n 
0002 | O | 1.10 19.7 | 223 | 
1010 | 90 1.03 22.7 | 24.2 ~« | 
1011 65.0 1.06 19.3 | 220 | 
1012 47.0 1.14 16.7 | 18.7 | 
1013 35.6) » | ‘s 
1030 90.0, 1.20 13.7 | 15.9 | 
1122 61.7 1.21 120 | 13.9 
2021 76.9 1.17 124 | 139 | 
2023 55.0 1.42 8.1 | 10.2 | 
1015 23.2 . | 
ian | 428 1.69 | 67 | 82 
2131 80 1.29 77 | 97 | 
0006 0 \ Ds | 
2132 70.6) 1.49 6.3 | 7.8 | 











4. TEMPERATURE FACTOR AND STRUCTURE 
FACTOR FROM DIFFUSE SCATTERING 
From the theory of diffuse scattering from 
single crystals of zinc'®: '7 we have 
S= Secon + Since/(1 +a vers ¢)', (11) 
where Seon = f?(1 —e-?™*) /Z (12) 
and . a = 1 —YE;?/Z. (13) 
In diffuse scattering experiments the inten- 
sities involved are so small that it has been found 


advantageous to use the continuous spectrum, 
and so the equation!’ 


Sp p 
( ) = J (Son-+Sine)=Ian (14) 
ML exp K 


has been developed. In this, J(A) is the distribu- 
tion of the energy in the continuous spectrum. 
The integral involving Sine can be calculated and 
(14) then is reducible to the form 





GO = [Gu/N10)a, (15) 
where ¢=(sin }¢), and G(t) and J(A) are experi- 


mental quantities. By a transformation of vari- 
ables, this takes the form 


nt) = f als), s)ds, (16) 


where A(t) is a known function of ¢ (an experi- 


16 G. E. M. Jauncey and W. A. Bruce, Phys. Rev. 50, 408 


7G. E. M. Jauncey and W. A. Bruce, Phys. Rev. 51, 
1068 (1937). 
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TABLE II. Values of F. 


6 7 8 9 10 
Af’ Af 
\ 1.54 X 1.54 

fo.44 Fy.sseM | (y 0.44) Ly 071) iH 
3.7 21.7 2.0 ? 86 2.43 
5.2 | 23.4 1.8 1.55 2.44 
3.9 20.4 3.5 2.86 2.44 
12 | 19.0 2.2 2.28 2.44 
8.7 | 16.4 2.3 2.64 2.45 
17.1 | 14.5 2.6 2.30 2.45 
16.7 | 145 2.2 2.75 2.45 
14.7 | 11.5 3.2 2.98 2.45 
13.7 | 11.3 2.4 | 2.54 2.46 
13.2 | 9,94 3.2 2.58 2.47 
12.3 | 9.4 2.9 2.24 2.47 





mental curve) and J(A) is a known function of A 
(also an experimental curve). Eq. (16) is an 
integral equation of the first kind.'* This equa- 
tion possesses a unique solution for our experi- 
mental conditions but there is no standard pro- 
cedure for obtaining the solution. However, a 
graphic solution has been obtained and is ex- 
plained in detail elsewhere.'’ By means of this 
graphic solution values of f and M have been 
obtained. These values of f are shown in column 
6 of Table II. These f values are assigned to an 
average wave-length of 0.44A. 

At a given temperature the temperature factor 
depends only upon the angle of orientation y and 
the scattering angle ¢, and the formula for M 
at 7 =295°K is given by the empirical relation 


M = (2.34 cos? y+0.68 sin? ¥)(sin? 3¢)/d*, (17) 


in which the constants of the equation are those 
found by Jauncey and Bruce.’ The trigono- 
metrical functions of the orientation angle are 
those given by Zener’s theory.’ The temperature 
factors obtained by use of (7) and the orientation 
angles in column 2 are listed in column 3. 

There are several advantages in the diffuse 
scattering method, despite the low intensity and 
the tedious mathematics. In a soft metal like 
zinc where it is possible for small particles to be 
strained by van der Waals forces, it is very good 
to be able to use large single crystals that can 
be annealed, after they are grown and polished, 
and thus made entirely strain free. 


18 See R. Courant and D. Hilbert, Methoden der Mathe- 
matischen Physik, p. 97. 
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5. RESULTS AND COMPARISON WITH THEORY 


The f values for \=1.54A are shown in column 
7 of Table II. These are obtained from the F 
values by multiplying the corresponding values 
shown in columns 3 and 4. Column 8 shows 
experimental values of fo.4s—/f1.54. Experimental 
values of fo.71—fi.s4 are shown in column 9. The 
theoretical values of fo.71—/fi.54 are shown in 
column 10. For these it is necessary to include 
the 7,“ term and there is a slight dependence 
on angle. However, the variation is smaller than 
the experimental error and so the average of 
the values shown in column 9 is compared with 
that of column 10 thus 


Af’ exp = 2.40 ; 
0.71A to 1.54A. 
Af’ theo = 2.45. 
The theoretical value of fo.4;— 1.54 is 2.45 and is 
practically constant with the angle. A com- 
parison with the experiment is shown thus 
Af'exp=2.57 
Af’ theo = 2.49 
This is a clear verification of the validity of 
Hoénl’s theory. It also gives one considerable 
faith in the f values of zinc as given by these 


Joasa to 1.54A. 


AND 


s. DBD. CORNELL 
TABLE III. Structure factors for zinc. 
(sin 4¢)/d 0.1 | 0.2 | 0.3 | 04/05 | 06! 08 | 1.0 
ffor ¥=0° | 27.0} 23.5) 20.2) 17.0) 14.4| 
} | | | 12s} 9.5 | 7.3 
f for ¥y=90° | 28.4) 25.9) 21.5) 18.2) 14.9 
| | 





four experiments and in the MW values as given 
by Jauncey and Bruce. Table III lists the aver- 
age of the experimental f values of the four 
experiments corrected for the change in structure 
factor due to the nearnessof the K absorption edge. 

It is interesting to note that the f values for 
(sin $¢)/A=0.2 for ~y=0° and 90° differ by 2.4 
electrons which is greater than 2. By reference 
to James and Brindley’s paper'® on the calcula- 
tion of quantum-mechanical f values, it is seen 
that at this value of (sin }¢)/A the scattering for 
the two valence electrons is negligibly small, and 
so this variation of 2.4 electrons cannot be due 
to the valence electrons as stated in a previous 
paper.'? The electron asymmetry thus extends 
deeper into the atom than the valence electrons. 

The author wishes to thank Professor G. E. M. 
Jauncey for his encouragement and guidance in 
this work. 


19 James and Brindley, Phil. Mag. 7, 81 (1931). 
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The Band Spectra of SrCl and SrH 


KENNETH R. More* Anpd S. D. CORNELL 
Sloane Physics Laboratory, Yale University, New Haven, Connecticut 
(Received March 30, 1938) 


The spectra of SrCl and SrH are observed in absorption. The SrCl analysis given by Parker 
is confirmed. Measurement of the relative intensities of corresponding bands of the electronic 
doublet *II—* indicates that 2 is the ground state. Rotational analyses of several bands of 
the C and D systems of SrH lead to values of B, of 4.008 and 1.925 for the C and D states, 
respectively. Vibrational analyses of the two systems give w,.=1206.2, x.w,.=17.0 for the 
X *Z state; w,.=1347, x.w,=23.5 for the C 23 state and w,=1014.1 and x.w,=15.4 for the D 2 
state. The perturbations in the C and D states are discussed in terms of the interactions of 
the C, D, and E states. 


co-workers. All but one of these systems have a 
strong (0,0) band, together with a weak (0,1) 


2W. R. Fredrickson, M. E. Hogan, Jr. and W. W. 
Watson, Phys. Rev. 48, 602 (1935). 

3W. W. Watson, W. R. Fredrickson and M. E. Hogan, 
Jr., Phys. Rev. 49, 150 (1936). 

*R. F. Humphreys and W. R. Fredrickson, Phys. Rev. 
50, 542 (1936). 


INTRODUCTION 
TOTAL of five band systems of the SrH 
molecule have been reported in a series of 
papers by Watson and Fredrickson!“ and their 


* Sterling Fellow. 
1W. W. Watson and W. R. Fredrickson, Phys. Rev. 39, 


765 (1932). 
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band, and very weak (1,1) and (2,2) bands. 
The exception to this is the D system, D? SX 22. 
Watson, Fredrickson and Hogan* report two 
bands of this system, corresponding to transi- 
tions from a common state of unknown v’ value 
to the ground state vibrational levels for which 
v’’ equals 2 and 3. No other bands of this system 
have been reported. The D state differs from 
the other states of the SrH molecule in that the 
B constant has a value of about 1.9, while the B 
constants for the other states vary from 3.3 to 
3.9. This same difference is found in the case of 
the corresponding states of CaH as reported by 
Grundstrém.® It is interesting to note that in 
the D system of CaH the bands of the v’=0 and 
v’’=0 progressions are more prominent than are 
the corresponding bands of the other systems, 
and that the (0,0) band is absent. It is also 
observed that the vibrational frequency for the 
D state is appreciably less than that in the 
ground state. 

The present investigation was undertaken 
primarily for the purpose of extending the 
analysis of the D system of SrH. Since the work 
on SrH just discussed was done in emission it 
was felt that it would be profitable to attempt 
the problem in absorption. During the course 
of the experimental work it became evident that 
some interesting information was being obtained 
about the C system of SrH, and about the band 
systems of the SrCl molecule. These latter bands 
were observed because of the presence of some 
chloride as an impurity in the material used for 
the core of the furnace in the earlier part of the 
work. Indeed the red bands were extremely 
strong and were well developed at relatively low 
furnace temperatures. 

A high temperature vacuum furnace which was 
constructed for such spectroscopic investigations 
was used for the study of the absorption spectra 
to be considered. The furnace will be described 
in another article. 

In order to obtain the spectrograms, light 
from a tungsten filament lamp operated con- 
siderably over its rated voltage was passed 
through the core tube of the furnace and focused 
on the slit of the spectrograph by means of a 
suitable optical system. For the study of the SrH 
and SrCl bands, strontium metal was distributed 


5B. Grundstrém, Zeits. f. Physik 75, 302 (1932). 


along the core tube and the furnace filled with 
hydrogen to a pressure of about 70 cm of mer- 
cury. In the case of the SrCl bands, zirconium 
oxide cores were used and the furnace heated to 
a temperature of 800° to 1000°C. The chlorine 
came mainly from the core because the chloride 
bands were much weaker when graphite cores 
were used. Temperatures of the order of 1600° to 
2000°C were found necessary for the observation 
of the D system of SrH. Graphite cores were 
used in this case. 

The spectra were photographed in the second 
order of a 21-foot grating in a stigmatic mount- 
ing. The dispersion was about 2.5A per mm. 
Eastman II-O and I-F plates were used in the 
regions below and above 5000A, respectively. 
Exposure times of the order of 2 to 5 minutes 
were ample. It was thus possible to obtain an 
exposure before the strontium metal all diffused 
out of the core. 


STRONTIUM CHLORIDE 


Two band systems of strontium chloride have 
been observed in emission by Hedfeld® and by 
Parker,’ and in absorption by Walters and 
Barratt.* The latter apparently used rather low 
dispersion. According to the analyses of Hedfeld 
and of Parker the red system is due to a transi- 
tion of the type *II—*2. It was not possible to 
tell from the emission data which of the two 
states involved was the ground state. 

Both band systems previously reported are 
present on our spectrograms. Our data confirm 
the analyses already published, and extend some 
of the sequences to higher members. It is of 
interest that sequences extending as far as the 
(12,12) band are observed in absorption at 
temperatures of about 1500°C. Since the vibra- 
tional frequency for the ground state of SrCl is 
300 cm, an appreciable number of molecules 
are excited to levels up to about 3500 cm™ 
above the level v’’=0. 

If the ground state of SrCl is *II, with the 
electronic interval of 607 cm~ which is observed, 
the relative intensities of corresponding bands 
of the two members of the doublet observed in 


°K. Hedfeld, Zeits. f. Physik 68, 610 (1931) 

7A. E. Parker, Phys. Rev. 47, 349 (1935). 

§Q. H. Walters and S. Barratt, Proc. Roy. Soc. A118, 
120 (1928). 
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absorption can be calculated from the Boltz- 
mann factor for the temperature used. For a 
temperature of 800°C the ratio of intensities is 
about 0.5 to 1. A spectrogram taken at a tem- 
perature which was certainly not greater than 
800°C has been photometered in order to de- 
termine the relative intensities of the two 
members of the electronic doublet. The Q heads 
of the Av=0 sequences were compared. Within 
the limits of accuracy the intensities were equal. 
This indicates that the ground state is *2. It 
should be pointed out that even though the plate 
used for this comparison was not calibrated the 
result should be reasonably accurate. The galva- 
nometer deflections produced by the background 
and by the absorption band heads, when com- 
pared with the clear plate and total blackness 
deflections, indicated that the blackenings in- 
volved fell on the straight line portion of the 
plate calibration curve. At any rate conditions 
were such that a true intensity ratio of 0.5 to 1 
could hardly give equal blackenings. The result 
that the ground state is *2 is in agreement with 
the work of Jenkins and Harvey® who report the 
ground states of CaF and SrF and BaF to be 22. 


STRONTIUM HYDRIDE 


A. The C system 

The (0,0) and (0,1) bands of the C *2ssX 22 
band system which have been reported by 
Fredrickson, Hogan and Watson? and _ by 
Humphreys and Fredrickson‘ are present on the 
absorption spectrograms taken with the furnace 
at a temperature of about 1500°C. These bands 
do not show the predissociation cut off at 
K'=19 which was reported for the emission 
bands excited at low hydrogen pressures. This 
is in agreement with the work of Humphreys 
and Fredrickson,‘ who extended the R branch of 
the (0,0) band to members as high as K’’=49 
by operating a strontium arc in hydrogen at 
pressures above 400 mm of mercury. It is of 
interest that the band lines obtained in absorp- 
tion are relatively sharp from the origin to 
K” =18, beyond which point they are somewhat 
diffuse. 

In addition to the (0,0) and (0,1) bands, the 
(1,0), (1,1) and (2,2) bands are present on the 
absorption spectrograms. The (0,0) band is much 
®F. A. Jenkins and A. Harvey, Phys. Rev. 39, 922 (1931). 
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more intense than the others. The (1,0) band 
was not observed in emission by the authors 
mentioned. It, of course, lies far beyond the pre- 
dissociation limit at K’=19 in the (0,0) band. 
Like the high K value lines of the (0,0) band, 
the lines of the (1,0) band are diffuse. The 
assignment of frequencies in the (1,0) band is 
given in Table I. This analysis was made with 
the help of the lower state combination differ- 
ences given by Watson and Fredrickson.! The 
weakness of this band and the presence of per- 
turbations make the somewhat 
doubtful in some cases. The (1,1) and (2,2) 
bands are very weak but the positions of the 
heads are reasonably certain. 

The rotational and vibrational constants of the 
C state are given in Tables V and VI, respec- 
tively. The B,’ value was obtained from an 
analysis of the (1,0) band, and the By’ value is 
that given by Humphreys and Fredrickson.‘ 
The vibrational constants are obtained from the 
vibrational array for the C system which is 
given in Table II. This array was constructed 
from head data. The lower state frequencies 
found from this array are in good agreement 
with those found from the D system which is 
discussed in the following. 


assignments 


TABLE I. Assignment of frequencies in the C *2<X * (1,0) 
band of SrH. 





=" R ad K”’ R P 
0 27619.5 12 45.0 49.9 
1 27.2 + 27604.3 13 58.0 49.9 
2 35.4 597.7 | 14 70.8 49.9 
3 43.7 91.0 | 15 85.1 49.9 
4 54.1 84.7 | 16 96.8 49.9 

99.4 
5 60.6 78.5 801.5 
6 67.2 73.3 17 19.9 49.9 
7 74.8 66.1 18 31.1 49.9 
8 85.5 59.5 19 44.2 57.8 
9 96.1 49.9 20 59.4 
708.5 61.5 
10 08.5 49.9 21 
20.6 
22 91.2 57.8 
11 28.0 49.9 59.5 
32.7 57.8 
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band TABLE II. Wave numbers of.the band-heads of the PaB_e III. Assignment of frequencies in the D %3<-N*d 
nen C?2<—X 2 system of SrH. band system of SrH. 

' pre- b Ny" 0 1 2 0,0) Banpb 3,0) Bann (4,0) Bann 
and. . K" R p ike R p ike R p 
jand, ol : | | 3 0 0 23614.6 0 245048 
rhe 0 | 26251 25080 2 20754.0 2 t 0 een ; 007 “— 
nd is | . | fw lf a all we 
with 1 | M880 26380 ; lt me Blt ee 
iffer- 26AOS 7 685.0 30.3 | 7 40.3 485.0 | 7 27.2 373.1 
<= Secs erste s ge Gals Bk Btls me 
ite : - 31.9 87.0 17.1 - ie _ 
had B. The D system . hn ps W357) ied | ant 22.2 
(2,2) Bands attributed to the D system were found “703 794 | 11 i. 1 sh 
f the throughout the region from 4000 to 5600A, ¥” 19 ile: tee done : 65.1 74.3 
with the strongest ones between 4080 and 4480A. +: wo. aenl S40 1 | 13 ms 7 
f the No bands were observed between 4480 and 496.6 86.9 | 13 oo 48 | ve site 
spec- 4720A because of the extreme breadth of the 69.7 os) on bend 76.1 71.6 
n an strontium resonance line at 4607A. The band - ax aaal 08.2 = 90.6 | 15 4 =e 
ue is lines from 4720 to 5000A were all weak. The 27.8 es os 6. pu iia 
son.4 band (v’,2) reported by Watson, Fredrickson and .< sa. oe 8 57 
1 the Hogan*® was present. The region of the (v’,3) 87.5 | 17 062.4 on 
h is band was not photographed. 7 03.0 36.2 | ” as iia 
icted With the aid of the known lower state com- 18 ISLA 22088.9 | 52.1 57.1 
ncies bination differences, and the combination differ- ~ a wonem oe “¥ 
ment ences for the level of unknown v’, six bands have 574 | on ati 
ch is been analyzed. | 20 906.2 | eanent Yat" 
The bands at 19580 and 20752 have upper 21 B30 | 
state combination differences equal to those _ ae s See 
(1,0) reported for the (v’,2) and (v’,3) bands, and (2,0) BaNp (0,1) Bann | (1,1) Bap 
a, lower state combination differences which corre- xk” K P | Ke k Pp | K K P 
, spond to values of 1 and O respectively for v’. 0 °° &| o | @ 20565.6 
peoeame The band at 20561 cm~ gives a value of 1 for v’’ 2 ees | 2 19581.7 2 628 | a4 
9.9 ; . ° ; ° 3 80.0 22653.6 | 3 75.3 195511 | 3 56.9 30.2 
and upper state combination differences slightly 4 - 36.3] 4 662 317] 4 474 13.3 
9.9 lower than those of the two already discussed. 6 387 5809/6 380 4880/6 Ite ops 
9.9 The bands at 22686, 23611, and 24500 cm, ‘ trp | 7 ge 8] FB ee 
- which incidentally are the strongest observed, $545. 302) 8 497. ; 3 - 
, all indicate that v’’=0, and that the upper ‘ mae 71.6 | 9 54.2 380.9 
9.9 state combination differences decrease uniformly ’ 660 OL 7 | _ ee 
in going from one band to the next. The wave ww 388 531 | - 7 243 123 
numbers 20752, 22686, 23611, 24500 and the ~~ =~ @& = | 1 3980 or 
- change in B’ indicate that these bands for which ™ 03.0 150 | ” nan 
9.9 v'’=0 belong to the v’’=0 progression with one 12 4675 ~—— 373.6 | a 638 623 
78 member missing. The missing member falls in _ sed 13 34.2 23.1 
the region of the strontium resonance line. It is ™ 32.3 334 | | 0 - 
also evident that the bands 20752, 19580, 18442) us oss | = 2 
(v’,2), and 17338 (v’,3) form a v’ progression. = 15 oes 31.8 
It is probable that this is the v’=0 progression. » ~ 7 
7.8 The rotational quantum assignments of the six io 337 BLS 
9.5 bands discussed are given in Table III. The i 06-2 
— vibrational array for the D system, using origin —— | 7 
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TABLE IV. Wave numbers of the band-origins of the 
D *Z<—X *E system of SrH. 








\ v” 0 1 2 3 
\ e — 
v’ 
0 20751.4 19580.0 18441.8 17337.5 
1 20561.8 
2 22686.4 
3 23611.0 
4 24500.5 








data, is shown in Table IV. The rotational and 
vibrational constants of the D *X and X °Y states 
are given in Tables V and VI. 

The origins of the bands of the D system are 
given by the following formula: 


= {1206.2(v’’ +3) —17.0(v’’+4)?} 


A comparison of the intensity distributions in 
the case of emission and absorption is of interest. 
The strongest bands observed in emission by 
Watson, Fredrickson and Hogan* were those 
analyzed and now called the (0,2) and (0,3) 
bands. The strongest bands observed in absorp- 
tion are the (3,0) and (4,0), followed in the order 
of decreasing intensity by (2,0), (1,1), (0,1), 
(0,0) and (0,2), the latter three being about 
equal. This indicates that the distribution may 
be represented by a broad Condon parabola. 
In the case of the absorption bands the (0,2) and 
(0,3) bands are weakened by the Boltzmann 
factor, since the vibrational energies of states 
for v’’=2 and 3 are 2310 and 3404 cm~ above 
the v’’=0 state. This effect is not present in 
emission. This intensity distribution is similar to 
that found by Grundstrém® for the D system 
of CaH. 

The Mecke rule, which states that B,/w, is 
constant for all states of a molecule, is violated 
by the D state of SrH, as it is for the D state of 
CaH. The ratio B,/w, equals 1.9X10-* for the 
D state, 3.1X10-* for the ground state and 
2.9 10-* for the C state. 

Most of the bands of the C and D systems 
show perturbations after the first few members 
of a band. The (1,0) band of the C system 
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shows perturbations similar to those reported for 
the (0,0) band by Humphreys and Fredrickson.‘ 
It was not possible to trace these perturbations 
in detail because the band was not extended to 
high enough values of K. The (0,0), (1,1), (2,0), 
(3,0), and (4,0) bands of the D system all show 
perturbations of the type illustrated by Watson, 
Fredrickson and Hogan’ in the case of the (0,3) 
band (called (v’,3) in their paper). At a K’ 
value of 10 to 15 the observed frequencies 
suddenly deviate by a marked amount from the 
frequencies calculated from a fourth power 
equation which fits the observed values up to 
the K’ in question. The deviations in some cases 
reach values of 50 to 100 cm. The values of K’ 
at which the perturbations begin are different 
in the different bands, but all fall in the range 
10 to 15. The states which interact with the C 
and D states to cause these perturbations, and 
the predissociation in the C state will be dis- 
cussed after the dissociation products and the 
heats of dissociation are considered. 


DISCUSSION 


In order to obtain some idea of the probable 
products of dissociation to be associated with 
each molecular state, approximate values of the 
heats of dissociation for the X, C, and D states 
were determined by the method of Rydberg. 
The values of the vibrational constants found 
for these states were used in this determination. 
Since the vibrational constants for the A ?II, 
B*>, and E?Il states of SrH are not known it 
was not possible to calculate the heats of 


TABLE V. Rotational constants of the C*Z and DE states 





of SrH. 
C % D*z Dx 

By=3.942 Bo= 914 B,=1.816 
B,=3.810 B,=1.890 a=0.024 
a =0.132 | Bs=1.865 B.=1.925 
B,=4.008 B;=1.841 D=-—0.3X10" 
D=—1.4x10 | 

| 





TABLE VI. Vibrational constants of the X 2Z, C2, and D *> 
states of SrH. 





X 2y C2 Dz 





(FROM ORIGIN DATA) (FROM HEAD DATA) (FROM ORIGIN DATA) 
we = 1206.2 we = 1347 we = 1014.1 
Xétdte= 17.0 Xéide™= 23.5 Séve= 15.4 
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dissociation for those states. However, some idea 
of the values was obtained from the probable 
dissociation products, once the results for the 
X, C, and D states were used to make a correla- 
tion with the known atomic states. The heats of 
dissociation are given in Table VII. 

The potential energy curves for the X, C, and 
D states, obtained by using the Morse function, 
are shown in Fig. 1. The states of the Sr atom 
into which the molecules dissociate from the 
various molecular states are shown at the right 
of the figure. Potential curves for the A, B, and 
E states are also shown. They were obtained from 
the known positions of the minima and the 
probable dissociation limits. 

It is seen that the D state is the only known 
state which can cause the predissociation in the 
C state. This determines the heat of dissociation 
of the D state, because the dissociation limit 
must be lower than the predissociation cut off in 
the C state. The value of the dissociation energy 
thus found for the D state is lower than that 
found from the vibrational constants by the 
Rydberg method. It thus seems reasonable to 
assume that the values given by the Rydberg 
method are high rather than low for the X and C 
states. By taking a value of the dissociation 
energy of the ground state somewhat lower than 
that calculated, it is possible to correlate the 
dissociation limits with the atomic terms as 
shown in Fig. 1. The probable configurations and 
dissociation products are given in Table VII. 
It should be pointed out in.this connection that 
if the dissociation energy of the ground state is 
higher rather than lower than the value calcu- 
lated from the molecular constants, the *P 


TABLE VII. Probable configurations and dissociation 
products of the SrH molecule. 


PROBABLI DISSOCIATION DISSOCIATION 





STATE CONFIGURATION PRODUCTS ENERGY ,(cm™!) 
X 5s0*4doa Sr, 1S+H, 28 10,000 

A II 5so*S pr Sr, *P+H, 2S 11,000 
By 5sa*5 po Sr, 3P+H, 2S 10,000 

E ?II 5s0*4dx? Sr, 3D+H, 2S 9,500? 
Dz 5sc5po4da Sr, 3D+H, 28 7,500 


C 2 5s0°6s0 


Sr, 3S+H, 2S 12,500 
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Fic. 1. Potential energy curves for SrH. 


state of the atom could be the dissociation limit 
of the D state. In this case the A, B, D, and 
perhaps the £& states would all approach the *P 
term as a dissociation limit. It seems reasonable 
however to take the correlation given because the 
dissociation energy in the D state is considerably 
lower than that calculated by the Rydberg 
method. 

From the relative positions of the potential 
energy curves shown it is seen that the levels 
of the D state are perturbed probably by levels 
of the E state. The curves of the C and D states 
do not actually cross, because levels of the same 
A value cannot do so. In view of the indicated 
shapes of the C and D curves, they merge in 
the manner shown. This merging of the curves 
occurs in the region of predissociation of the C 
state. The predissociation is to be accounted for 
by this interaction of the C and D states, to- 
gether with the relative positions of the point of 
merging and the dissociation limit of the D 
state. 

In conclusion the authors wish to express their 
appreciation to Professor W. W. Watson for 
many stimulating discussions of this work. 
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On the Saturation of Nuclear Forces 


C. CRITCHFIELD AND E. TELLER 
George Washington University, Washington, D. C. 
(Received March 12, 1938) 


The forces between heavy elementary particles, i.e., neutrons and protons, are investigated 
on the hypothesis that they emit electron-positron pairs. It is also assumed that the interaction 
of the heavy particles with the field of light particles is large compared with the kinetic energy 
of the light particles. It is shown that potentials result which are of the same order of mag- 
nitude as the kinetic energy of the light particles. When many heavy particles interact the 
total potential energy is found to be proportional to the number of heavy particles. 


INTRODUCTION 


HE mass defect of most atomic nuclei is 

roughly proportional to the number of 
heavy elementary particles contained in the 
nucleus. This shows that the binding energy per 
elementary particle is, in first approximation, 
independent of the number of particles consti- 
tuting the nucleus. 

Since the heavy particles obey the Fermi-Dirac 
statistics the average kinetic energy of a heavy 
particle will vary with the two-thirds power of 
the density whenever the density is sufficiently 
high and the velocities of the heavy particles are 
small compared with the velocity of light. Thus, 
under these conditions, the kinetic energy in a 
given volume element of a heavy nucleus will be 
proportional to the five-thirds power of the 
density. If, on the other hand, the total potential 
energy is a sum of potential energies acting 
between pairs of heavy particles and if the range 
of the potential energy is small compared with 
the nuclear radius the total potential energy in a 
given volume within the nucleus will be pro- 
portional to the square of the density. If attrac- 
tion is assumed between all heavy particles, the 
binding energy would then increase with in- 
creasing density and stability could not be 
reached before the nuclear radius becomes as 
small as the range of the attractive forces. In this 
case, however, the mass defect for heavy nuclei 
would increase more rapidly than the number of 
heavy particles contained in the nucleus.! 

In order to escape this difficulty Heisenberg? 
and Majorana have assumed an _ interaction 


1 These arguments were first put forward by Heisenberg 
at the Solvay Congress, 1933. 
2 W. Heisenberg, ‘“‘Rapports de Congres Solvay, 1933.” 


81 


between heavy particles which is of the exchange 
type. For such interactions it follows that a 
neutron or proton can only interact with a 
limited number of other neutrons or protons. 
The following mechanism for this exchange 
interaction has been proposed :* a neutron emits 
an electron and a neutrino and turns into a 
proton; a proton absorbs the emitted pair and 
becomes a neutron. Thus an exchange of charge 
between a neutron and proton takes place. 

This picture, however, cannot explain in a 
simple way why the forces between two protons 
(between which no such simple exchange of 
charge can take place) is of the same order of 
magnitude as the force between a neutron and 
proton.*: 5 Moreover, the probability of electron- 
neutrino pair emission can be estimated by 
extrapolating the probabilities of beta-decay 
processes. Performing this extrapolation either in 
accordance with the theory proposed by Fermi® 
or with that proposed by Uhlenbeck and 
Konopinski’ one obtains a too small value for the 
interaction between heavy particles. 

Gamow and Teller§ and Wentzel® have there- 
fore suggested that an emission and reabsorption 
of electron-positron pairs should be responsible 
for the nuclear forces. The equality of proton- 
proton and proton-neutron forces is then more 
easily understood and also sufficient freedom is 
obtained to explain the magnitude of the nuclear 

3D. Iwanenko and I. Tamm, Nature 133, 981 (1934). 

4M. A. Tuve, L. R. Hafstad and N. P. Heydenburg, 
Phys. Rev. 50, 806 (1936). 


5G. Breit, E. Condon and R. Present, Phys. Rev. 50, 
825 (1936). 

6 E. Fermi, Zeits. f. Physik 88, 161 (1934). 

7K. Konopinski and G. Uhlenbeck, Phys. Rev. 48, 7 
(1935). 

8 G. Gamow and E. Teller, Phys. Rev. 51, 289 (1937). 

®G. Wentzel, Helv. Phys. Acta 10, 107 (1936). 
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forces. Since, however, according to this hy- 
pothesis the charge is not exchanged and since a 
study of the neutron-proton interaction forces 
has eliminated the possibility that the major 


‘part of the neutron-proton interaction'’ is due 


to spin exchange we must conclude that if 
nuclear forces are due to electron-positron pair 
emission the main part of the forces between 
heavy elementary particles is not of the exchange 
type.'"' The purpose of the present paper is to 
show that the interaction between heavy parti- 
cles resulting from the pair-emission hypothesis 
is compatible with the proportionality between 
mass defect and number of heavy particles. 

In order to do this it will be necessary to drop 
the assumption that the total potential energy 
within the nucleus can be represented as a sum of 
interactions between pairs of heavy particles and 
that the potential energy is proportional to the 
square of their number. Indeed, this is correct 
only if the virtual states into which the light 
particles may be emitted are occupied only witha 
very small probability. If, however, the term in 
the Hamiltonian, 7/7’, which gives rise to the 
creation and annihilation of the light particles is 
large compared with the kinetic energy, Exin, of 
the light particles in the virtual state a saturation 
phenomenon may result and the potential energy 
may increase less strongly than with the square 
of the number of particles. It then follows from 
the argument given in the beginning of the 
introduction that if the potential energy increase 
less strongly than the five-thirds power of the 
density the density of heavy nuclei will not 
depend strongly on the number of particles 
within the nucleus, and the mass defect divided 
by the number of particles in the nucleus will 
remain roughly constant. 

We shall assume therefore, in the following 








10 Tf the neutron-proton interaction be due to pure spin 
exchange there could not be an attraction between a 
neutron and a proton for both parallel and anti-parallel 
spins. 

“If the probabilities of electron-positron, electron- 
neutrino and positron-neutrino emissions are of the same 
order of magnitude, an exchange of charge would remain 
possible and there would be no necessity for abandoning 
exchange forces. This suggestion is implied in Feenberg’s 
discussion (Phys. Rev. 51, 777 (1937)) of heavy particle 
interactions. The order of magnitude of nuclear forces 
would then still require an explanation (see Kemmer, Phys. 
Rev. 52, 906 (1937)) which can, perhaps, be given in terms 
of heavy electrons and neutrinos (see Bhabha, Nature 141, 
117 (1938) and Kemmer, Nature 141, 116 (1938)). 
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discussion, that J/'>F yin. This assumption is 
just the opposite of the one usually made, 
namely /]’<E,in, which served to justify the 
treatment of the emission and absorption of 
light particles by a perturbation method. 

In order to simplify the treatment we assume 
that not only the charge but also the spin of a 
heavy particle remains unchanged if light parti- 
cles are emitted. Furthermore, the Coulomb 
interactions will be neglected. Thus, finer features 
such as spin dependence of nuclear potentials and 
the difference between proton and _ neutron 
mass will not be obtained from the present 
development. 


Part I. ONE HEAVY PARTICLE 


We consider at first one heavy particle and 
discuss the effect of the operator //’ giving rise to 
emission and absorption of electron-positron 
pairs. The kinetic energy of the light particles, 
Exin, will be neglected in accordance with the 
hypothesis Eyin<H’. The effect of Exin will be 
taken into consideration as a perturbation in 
Part ITI. 

It is convenient to adopt the “hole theory” of 
positrons in developing the formalism for the 
emission and absorption of pairs. Thus, emission 
of a pair will correspond to the transition from a 
state of negative to a state of positive energy, and 
an annihilation of a pair, a transition from a 
positive to a negative state. 

Since we assume that a heavy particle after 
emission is in the same state in which it was 
before emission of a pair there is a possibility that 
a second pair be emitted before the first is re- 
absorbed. In general, a single heavy particle 
might emit successively in this way any number 
of pairs without reabsorbing them. If, however, 
pairs can be emitted only into a limited number 
of states, the exclusion principle will limit the 
number of pairs which can be emitted. 

In particular we shall suppose the interaction 
IT’ to be such that emission (or absorption) can 
take place only into four states Aft, AJ, Bt, BJ. 
At and AJ are two spherically symmetrical 
electron states with positive kinetic energy, 
differing only in spin direction. B¢ and BJ on the 
other hand, are states of negative kinetic energy. 
The absence of an electron from B¢ is equivalent 
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to the existence of a positron in a corresponding 
state with its spin pointing downward. We shall 
assume that the dependence on space coordinates 
of the A’s and B’s is the same. In general all 
negative states orthogonal to the B’s shall be 
assumed as filled. Bt or BJ, however, shall be 
filled only if the factor Bt or BY explicitly appears 
in the proper function. 

The formalism used in quantizing the wave 
equation will be adopted. Associated with the 
electron wave function Af is the operator at 
which acts on the functional y, i.e., on the wave 
function of all electrons. aty is zero if the state 
At is represented in y as filled and aty=(y, Af) 
if Af is empty. (y, Af) is the properly antisym- 
metrized product wave function of y and Af. 
This means that af gives rise to the appearance 
of an electron in the state A? if At be empty and 
is equivalent to zero if Af be full. In accordance 
with the formalism, at* acting on (y, Af) gives y 
and acting on y gives zero, so that af* causes the 
disappearance of an electron from Af if Aft be 
full and is equivalent to zero if At be empty. In 
the same way aj, aj* are connected with AJ; 
bt, b¢* with Bt; and dj, b)* with By. For nor- 
malized functions the commutation laws hold: 


xx*-+-x*x=1, xy+yx =0, 


xy*+y*x=0 (xF#y), x*¥y*+y*x*=0, (1) 
x, y=at, ay, bt, df. 
The term giving rise to pair emission and 
absorption in the Hamiltonian will then be 
IT’=a0, (2) 


where @ is a positive constant and the operator O 

may be written” 
O=apht*+btat*+aypoy*+byay*. 

The first and third terms represent the produc- 


tion, the second and fourth terms the absorption 
of an electron-positron pair. 

If we write the pair emission operator, O, as the 
sum of two commuting operators, Of and Oy: 
Ot =atht*+btat*, O=O01t+0), 
Oj=afpby*+bjay*, Of0),—OjO0t=zero, 


iteration of Of gives, because of (1), 


(3) 


(4) 


12 A somewhat more general form for H’ may be obtained 
by adding A(afaft*+ala}*+bfot*+5)b)*) where X is a 
new constant. However, as long as ¢—A>>Exin Our main 
argument remains unchanged. 
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O04 =atbht*btat*+btat*atot*, 
(Ot)? = Of and similarly 
(O4)*=O4. 
The proper values of Of and OJ are therefore +1 
and zero. 

It is readily seen that the proper functionals of 
Ot belonging to the proper value —1 will be 
those for which the state (A¢— Bt) /v2 is full and 
the state (A¢+B)/v2 is empty. For the func- 
tionals belonging to the proper value +1 the 
state (At—Bt)/v2 is empty and (At+B?)/v2 is 
full; for the proper value zero both states 
(At—Bt)/v2 and (At+Bt)/v2 have to be full 
(i.e. there is an electron in Af but there is no 
corresponding positron present) or both have to 
be empty (i.e., we have a positron but no 
electron). All these proper functionals are de- 
generate since the presence of electrons in any 
state orthogonal to Aft and Bf will have no effect 
on the proper value. Similar statements are 
obtained for OJ if AJ and By are used instead of 
At and Bt. 

Since Of and Of commute the proper values of 
O=Ot+0O] will be the sums of the proper values 
of Ot and OJ. Thus we obtain +2, +1, and zero. 
The corresponding proper functionals will be 
products of the proper functionals of Of and Oj. 
The proper functionals for the lowest proper 
value, —2, will be characterized by the states’ 
(At—Bt)/v2 and (AJ—By)/v2 being filled and 
the states’ (A¢+B)/v2 and (AJ+By)/v2 being 
empty. We again have degeneracy since any state 
orthogonal to Aft, AJ, Bt, BY may be filled or 
empty. This degeneracy means merely that the 
state of the nucleus will not depend upon 
extranuclear electrons or positrons. The lowest 
proper value of //’ is therefore —2¢. Since a is 
assumed to be large all other states lie high and 
therefore will not be excited. We shall be con- 
cerned only with the lowest state." 

We propose to consider the energy —2e as part 
of the proper energy (or mass) of the proton or 
neutron. Because of its negative value it cannot 
represent, of course, the complete proper energy. 


13 The exclusion of the proper value zero means that the 


wave functions of atomic electrons must be orthogonal to 
[Af] and [AJ]. This amounts to a boundary condition 
which all atomic wave functions must satisfy near the 
nucleus. If the spatial extension of [Af] and [AJ] is very 
small compared with that of an atom the influence of this 
boundary condition will be small. It might account for an 
essential part of the isotope effects in atomic spectra. 
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SATURATION OF 
This energy, — 2c, can be, however, incorporated 
into the proper energy only if it can be shown 
that the corresponding part of the energy for N 
particles (neutrons or protons) will be —2Noe 
independently of their configuration. This proof 
is given in the next section. 


Part II. SEVERAL HEAVY PARTICLES 


The operator /I'(N) corresponding to the 
emission and absorption of pairs by N heavy 
particles will be the sum of operators //;’ 


N N 
H'(N) => HM, =cDO,, (5) 


k=1 k=1 


where //,’ corresponds to the emission and 
absorption of pairs by the &th particle. We may 
thus write: 


O; = a,tb;.t* +b,.ta;t* +aygdi4* + b.ary*, (3k) 


where a, and a,* are operators associated with 
the electron state A,f which is the same spherical 
function as Af with its center being located at the 
kth particle. Similar statements hold for a,J, a,4*, 
bf, O.t* and by), b.4*. The proper functions of (5) 
are not, in general, simple because A; is not 
orthogonal to A,,f, and so on for AJ and the B’s. 

However, two particular cases can be in- 
vestigated at once: the first is that in which the NV 
particles are far enough apart so that all A’s and 
all B’s can be considered orthogonal. Then the 
lowest proper value of //’(N) is the sum of the 
lowest values of //;’ or just —2No. The proper 
functionals belonging to this value are those 
functionals in which all (A,t—B,4)/v2 and 
(Ai —Bit)/V2 are filled and all (Axyt+B,4)/v2 
and (A;J+B,4)/v2 are empty. 

The second simple case is obtained when the 
positions of all heavy particles coincide" and all 
O;.’s become identical with O;. Then (5) becomes 


H'(N)=No0, 


and the lowest proper value is again — 2 Nc. In the 
corresponding proper functionals (A if—Bit)/v2 
and (A,J—B,))/v2 are filled and (A ,¢+B,t)/v2 
and (Ai4+B,)) are empty. 


4 The fact that the Pauli principle would prohibit more 
than four heavy particles’ occupying the same point is 
important in the determination of the behavior of particles 
in nuclei but not of essential interest in this consideration 
of the potential in which they move. 
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To show that the lowest proper value of //’(.V) 
is always —2No we define the quantities 


tina f Aut? ntdv= f Aust mydv 
= { Bit*Batdo= [ Bu*Bade (6) 
with eu, = 1. 


The commutation relations may be written: 
LX + Xm Xk = Emks 
Xk,» Ve=Art, Anh, Of, Or Deh, (7) 
XV * +Vm*x, = zero, xF¥y etc. 
It is possible to choose a set of 4N orthogonal 


states Ayf, AaJ, Bt and Bf related to A,t, Ax, 
Bt and B,J respectively by the linear relations: 


X= DLTaXy, X=At, Al, Bt, or Bh, 
d 


Xy,= DT 1Y,. (8) 
k 


Operators x, and x,* associated with Y, are 


defined by: 
x)= > Th Ix, . x= } T, WX, 
k d 
x=at, aj, bt, or bj, (9) 


at= (Ty) *x*, mt=Ta*n*. 
k " 


It will be convenient to take for 7, an 
orthogonal transformation,'® so that 


Ta* =Ty". (10) 


In fact, under such a transformation, J/’(N) 
retains the simple form that it has in (5). 


H'(N) =e ST wTi,*(antd,t* 


k dy, Pp 
+bh,ta,t*+a,]b,4* +),4a,4*) 
=o) 5,°(a,tb,t* (11) 
A. Pp 


+h ta,t* +a4b,4* +h4a,4*) 
=> a>-O,, 
r 


where O,=(a,th,t*+),fa,t* +a,gay*+ajay*). 


15 Tt will be seen that 7, is the transformation which 
brings the Hermitian matrix {¢,,«} to principal axes. 
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On the other hand, if the transformation T is 
orthogonal it follows that the orthogonal set of 
functions X, cannot be normalized. This is so 
because, if we assume the set X, both orthogonal 
and normalized, then any orthogonal transfor- 
mation (8) would yield a set X; which would be 
both orthogonal and normalized which is in 
contradiction to the above statements. We there- 
fore define the scalars c, by the relations'® 


f XXdo= 59, X=At, Aj, Bt,or By, (12) 


leading to the commutation relations 


tit,” +x,*x) = > Tin* a = €6\4, 
k,m 


(13) 
OV*+y,*x,=zero, x=at, aj, bt, or by 


(x #y) etc. 


The operator O, may be again considered as the 
sum of the commuting operators Oy\f and OyJ: 


Ost =aytht* +d,Fat*, O,=Ost+On, 


Ong =A, ray* +o tary’, 
OnstOnd — OntOnd = zero. 


(14) 


The proper values of Oy will then be the sums of 
those for O,f and O,{J and the proper functionals 
will be product functions of those for O,f and 
Oy). From (13), iteration of O,f shows that 


(OyF)? = 6,70 ie 
at) on" at, (15) 
(Oy$)? =cn7O>DK, 


so that the lowest proper value of each is —c, and 
the lowest proper value of O, is —2c,; this value 
belongs to those functionals in which the states 
(A,t—Byt)/(2e,)! and (Ay¥—Byal)/(2q)! are 
filled and the states (A,t+By,4)/(2c,)! and 
(Ay¥+By4)/(2c,)! are empty. . 

According to the commutation relations (13) 
the operators O, commute and 


O,O, — O,O, = zero 


and therefore the proper values of [/’(N) will be 

given by o times the sums of the proper values of 

the O,. In particular, the lowest proper value will 

be —2¢>-cq and the corresponding functionals 
N 


will be those in which all (A,#—By4)/(2c,)? and 


16 The numbers ¢ are the proper values of the matrix 
{emx}. 
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(A,}—ByJ)/(2c)! are filled and all (A,¢+ By) / 
(2e,)! and (A,¥+By,4)/(2q)! are empty. From 
(13) 


Da=D, DTT mr*€im = 2 eee = N. (16) 
d k 


A kim 
The lowest proper value of //’(N) is then —2No 
and independent of the relative position of the 
heavy particles. 


Part III. KInEtIC ENERGY or LIGHT 
PARTICLE FIELD 


We now introduce the kinetic energy of the 
light particles, K, as a perturbation in the 
calculation of the proper values of 


IT=H'+K. (17) 


We are interested in the lowest value of /7. The 
lowest level of /7’ has been seen to be an infinitely 
degenerate one of energy —2c¢. The introduction 
of K as a perturbation splits these levels and the 
problem of interest is to find the lowest of the 
resulting states. 

We shall adopt the usual convention of putting 
the kinetic energy of the state in which no 
electrons or positrons are present equal to zero. 
The average value of K in the state y is therefore 
calculated by finding /y*Kydv and subtracting 
from this value the average value of K in all 
negative energy levels. The average value for the 
lowest proper functions of //’ is then 


bf (4t—Bay*K(At— Bade 
+3 { (A}-BY*K(A—By)do +i 


— [ BYKBtde— [ BYKB Udi, (18) 


where xk is the average value of K in all states 
orthogonal to A’s and B’s and kp» is the average 
value of K in all negative states orthogonal to 
the B’s. 

According to the perturbation theory the 
lowest energy value into which a degenerate state 
will split is given by the minimum average value 
of the perturbation for the set of degenerate 
proper functions under consideration. It is evi- 
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dent that the matrix of K, calculated in ac- 
cordance with the preceding paragraph, is posi- 
tive definite. Also the matrix of «— xo is positive 
definite with the minimum value zero for that 
state in which no positive levels are occupied and 
all negative levels are full. Therefore the lowest 
proper value of H will be, in first approximation, 
the lowest value of I’, i.e., —2¢, plus (18) with 
k= Ko. Since the space dependence of Af, AJ, Bt, 
and BJ is taken to be the same 


farxa tde= f 4 KA Jdv 
~—_ f Bt*KBtdv = — f BY*KBido 
and the integrals in (18) may be expressed by 
Exin = 2 f arKa fdv. (19) 


and we obtain 
Emin = —2o+ Lixin. (20) 


In a representation including the A’s and B’s, 
K also has nondiagonal elements. The order of 
magnitude of these elements is the same as that of 
Exin. For example, the element in K coupling the 
state in which (At—B?)/v2, (A}—BJ)/v2 are 
full and (At+Bt)/v2, (A¥+BJ)/v2 as well as all 
other positive energy states are empty with the 
state in which (At—B?)/v2, (A}J+BJ)/v2 are 
full and (A¢+Bt)/v2, (A}—BJ)/v2 as well as all 
other positive states are empty is readily calcu- 
lated to be }Exin. However, the diagonal elements 
for the two states differ by 20. According to 
second order perturbation formulae the correc- 
tion to the calculation of the minimum energy 
associated with this transition will be F*,;,/8c. 
Similar transitions to electron states orthogonal 
to the A’s and B’s involve an energy difference 
approximately equal to c. We may then say that 
expression (20) is correct if contributions of the 
order of (E*,in/o) and higher orders may be 
neglected. 

The order of magnitude of the additional 
energy, Exin, can be estimated if the spatial 
extent of At is known. If the characteristic linear 
dimension of A? (its effective radius for example) 
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be “‘d’’ then for sufficiently small values of d: 
Exin~cp~ch/d. The significance of this energy is 
seen in considering many heavy particles. 

For N heavy particles the total Hamiltonian is 


H(N)=C[H/]+2XK;, (21) 


where K; is the operator A for the ith heavy 
particle alone. If contributions of the order 
F*yin/o are neglected the lowest proper value of 
(21) is, by the same reasoning which led to (20), 


Emin = 2No+ fwxK aude 


— [ eXKade, (22) 


where Wo is that functional belonging to —2.Ne in 
which all negative levels orthogonal to the B’s are 
filled and all positive levels orthogonal to the A’s 
are empty. x is the functional of negative states 
including the B’s. The average value of K 
indicated in (22) cannot be calculated in general 
without a specific knowledge of the dependence 
of Af, etc., on space coordinates. Since the form 
of this dependence is not known and since we 
make no hypothesis regarding it, we are limited 
to the evaluation of (22) for those special cases 
which are independent of the shape of the 
electron functions. 

If the N heavy particles are sufficiently widely 
separated the electron functions (A;, etc.) 
pertaining to each heavy particle may be con- 
sidered orthogonal. Then the average value of K 
will be found lowest in that state in which all 
positive levels orthogonal to the A,’s as well as all 
(A.t+B.4)/V2, (Arb +Bih) are empty and all 
negative levels orthogonal to the B,’s as well as 
all (A.f—B,f)/v2, (Arb —Bid)/V2 are filled. The 
average value of }-K,; will be NE, in: 


Emin( ©) = —2No+NExin. (23) 


The other simple case occurs when the N heavy 
particles may be considered to occupy the same 
point in space. The proper function of mini- 
mum energy in this instance is that in which 
(A,t—B,)/v2, (Ai) —B.))/v2 and all negative 
levels orthogonal to B,;t and B,J are filled and 


(Ait+Bit)/v2, (Ai +Bi4)/v2 and all positive 
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levels orthogonal to Ait and A,J are empty. As 
before these A’s and B’s are chosen arbitrarily 
from the JN identical possibilities. Then 


Emin(O) = —2No+Exin. (24) 


CONCLUSION 


The difference between E,,j,(%) in (23) and 
Emin(O) in (24) is the binding energy obtained if 
the N particles are brought from infinite dis- 
tances to the same point. Even if the positions of 
all heavy particles do not coincide the value of 
Emin(O) will remain essentially unchanged as 
long as the mutual distances remain small com- 
pared with the linear dimension, d, of the wave 
functions A and B. 

The interaction which we have obtained differs 
essentially from the ordinary interaction between 
pairs; if the interaction is the sum of interactions 
between pairs, then by bringing the heavy 
particles close together the potential energy will 
be proportional to the number of pairs, that is, to 
}N(N—1); for the type of nonexchange forces 
which we have discussed, however, the interac- 
tion at close contact will be Emin( ©) —Emin(O) 
=(N—1)Exin. Since, for finite separation, no 
lower potential is to be expected than for 
coinciding heavy particles, it follows from the 
discussion given in the introduction that the 
densities of heavy nuclei will remain essentially 
constant and that their binding energies will be 
roughly proportional to the number of heavy 
particles which they contain. The maximum 
potential energy per heavy particle will be Exin. 

It is evident that the range of the nuclear forces 
is of the same order of magnitude as the linear 
dimension, d, of the wave function A or B. From 
a study of nuclear radii it follows that this range 
is of the same order of magnitude as the electron 
radius e?/mc?. This means that the wave function 
A or B will have dimensions comparable to the 
electron radius and that in the pair emission 
theory of nuclear binding electron and positron 
wave-lengths which are small compared with the 
electron radius do not play an important role. 
Putting d=e*/mc? we obtain Exin~137mc*. This 
energy is of the same order of magnitude as the 
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minimum of the potential energy per particle 
usually assumed in nuclei. 

The formalism as given in this paper does not 
explain the spin dependence of the nuclear forces 
and does not give the magnetic moments of 
proton and neutron. An explanation of these 
phenomena might be forthcoming if it is assumed 
that the spin of the heavy particle may change 
during the emission of a pair. A dependence of 
the spin of the emitted light particles on the spin 
of the heavy particle has to be introduced. If, 
moreover, a heavy particle can emit electrons 
into several orthogonal states simultaneously, 
sufficient freedom may be gained to describe our 
empirical knowledge of nuclei. 

The form of the wave functions into which 
emissions can take place remains, in the present 
state, arbitrary and thus the detailed dependence 
of the potential on the configuration of heavy 
particles cannot be calculated. It seems, however, 
that under rather general assumptions about the 
functions A and B, the first interaction between 
two heavy particles approaching each other will 
be a repulsion. At small distances it has been 
shown that attraction is obtained. The low 
binding energies obtained for nuclei like Li® or 
B® in which, owing to the Pauli principle, the 
average distance between particles will be greater 
than in Be* or C” may be due to this repulsive 
potential at greater distances.” 

It is felt that the most serious objection to the 
theory proposed in this paper is the introduction 
of the finite distance operator O in Eq. (2). Such 
an operator may cause creation and annihilation 
of electrons or positrons at points different from 
the position of the heavy particle and no way is 
known for reconciling such processes with the 
principle of relativity. This means that the 
present theory will meet with difficulties if one 
tries to extend it to relativistic velocities of the 
heavy particles. 

In conclusion the authors wish to express their 
thanks to Professors Bethe, Fermi, Gamow and 
Oppenheimer for helpful discussions. 


17 The questions mentioned in the two previous para- 
graphs will be discussed in detail elsewhere. 
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low free oscillations in the limiting case of eccentricity oscillations, but as well of the oscillations forced 
i’ or unity, and in obtaining formulas which give by a uniform oscillating electric field, such as 
, the rough approximations to the frequencies and that of an exciting electromagnetic wave of 
eater decrements for eccentricities a little less than length long compared with the shorter axis of 
ilsive unity. Macdonald? also has obtained an expres- the spheroid. Although the methods developed 
sion for the fundamental frequency of the free are applicable to all the harmonics, we shall 
o the oscillations of a straight antenna of negligible limit the analysis in this paper to the funda- 
ction diameter by treating it as the limiting case of a‘ mental or first harmonic and, to the approxi- 
Such cone of very small vertical angle. Obviously mation to which they are needed for the solution 
ation this is a very poor physical approximation, and of the forced case, to the third and fifth har- 
frogs therefore it is not surprising that Macdonald’s monics. The spheroid will be supposed to be 
ow te value of the frequency differs from Abraham’s perfectly conducting, the only resistance taken 
‘the by over twenty percent, and from the experi- into account being that of radiation. Heaviside- 
the mental value by nearly as much. Recently Lorentz symmetrical units are used throughout 
en Stratton* has investigated the solution of the the analysis. 
f the scalar wave equation in spheroidal coordinates, 
but his differential equation is not the one to THE DIFFERENTIAL EQUATION 
their which the problem here considered leads. Ap- Let x, y, z be a right-handed set of rectangular 
- and parently no one has made much Progress 1 coordinates with the x axis along the long axis 
solving the problem of the forced oscillations of a of the prolate spheroid, and put p?=y?+2%. We 
prolate spheroid, which are, of course, of far introduce the prolate spheroidal coordinates 
mn more Practical importance than the free oscilla- £, n, @ by the transformation 
tions in the application to the antenna. 
‘Abraham, Ann. d. Physik 66, 435 (1898) ; 2, 32 (1900). “a we . (1) 
wee Electric Waves (Cambridge University p=f}(1—#)(n*—1)}?, (2) 
*Stratton, Proc. Nat. Acad. Sci. 21, 51 and 316 (1935). tan ¢=2/y, (3) 
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Fic. 1. Spheroidal coordinates. 


where —1<£ &<1, 15 n< @, O< ¢<27z. The co- 
ordinate surfaces £=const. constitute a set of 
confocal hyperboloids of revolution about the 
x axis, and the coordinate surfaces »=const. 
are a set of confocal prolate spheroids of eccen- 
tricity 1/n. The constant f represents the semi- 
interfocal distance, which is approached by the 
semi-major axis a of the prolate spheroids as 7 
approaches unity. The unit vectors &), mi, ¢1, in 
the directions of increasing £, n, ¢, respectively, 
constitute a right-handed set in the order 
named. The radius vector 7, whose square is 
r? = f?(n’?+—1), becomes effectively equal to fn 
at great distances from the origin O, and the 
cosine of the angle @ which the radius vector 
makes with the x axis approaches £ when r 
becomes very great. A section through the x axis 
of the coordinate surfaces is shown in Fig. 1. 
We are interested in fields whose time de- 
pendence is specified by the factor exp (—iw/). 
Evidently the only nonvanishing field compo- 
nents are E:, E,, Hy. If we put A=pH,, it has 
been shown by Abraham! that the field equations 





give 
u\1 1 0A 
= i(*) = neceapaeneseanrass “et, (4) 
KJ ef ((1—&)(n?—-&)}? On 
u\}1 1 OA 
B= -i(“)—— —__—_———_—, (5) 
KJ ef {(n?—1)(n?—&)}! 0€ 
1 1 
H,=- ——_—_——-A, (6) 
f(Qd-#)(W’-1)}! 
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I. ADAMS, JR. 


where A satisfies the differential equation 
(1—£)0°A /0#+ (mn? —1)0°A /dr? 

+2(7—-2)A=0. (7) 
Here «x and wu are respectively the permittivity 


and permeability of the medium surrounding 


the spheroid, and 
€=(w/v)f=2zf/d, (8) 


where v=c/(xu)! is the normal wave velocity in 
the medium, c represents the velocity of light 
in vacuum, and X is the wave-length. 

The variables are separable in (7), and if we 
put A(é, 7) =X(£) Y(m) we obtain the two ordi- 
nary differential equations 


(1-—#)@X/d?+(a—€F)X=0, (9) 


(n?—1)d?Y/dy?+(—at+en*?)Y=0, (10) 
where a is the constant of separation. 
SOLUTION OF THE EQUATION IN & 
First we solve (9). Putting 
X (€) =(1—#) !u(é), (11) 
we find 
d du u 
_{A-8)—| ~—tau= eu, (12) 
dé dt} 1-2 
from which we infer that the characteristic 


solutions 41, Ue, -** u,, +++ form an orthogonal 


set of functions. We express the characteristic 
values and the characteristic functions as power 


series in é’, viz. 
a=antane+areé+are'+aues+:--, 
and 
ui(E) = Uo(E) Han(E) e+ 2(E) A+ uis(E)e+---. 
Evidently 
ay =l(l+1), J=1, 2,3---, 
and uj» is the associated Legendrian polynomial 
uo = Pu(&) = (1 —&) (d/dé) P(E), 


where /=1 represents the fundamental or first 
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harmonic, /=2 the second harmonic, /=3 the 
third harmonic, and so forth. 

By the usual perturbation methods we find 
for the fundamental 





1 4 
a, =2+-é— é 
5 58-7 
8 124 
+——° —- ————e*+ > (13) 
3-55-7 5°-73-11 
1 
u,= Py(&) — ‘ P3i(é)ée 
3-5? 
2 1 | 
+} — P(t) +——— Pa (8) ¢ € 
32. 54 32. 52-72 
31 + 
saat magn nol ey nn ae) 
13.58.72 3-54. 72-13 
. | 
- —P7(£) ¢e6+ (14) 
3¢-5-72-11-13 
For the third harmonic we get 
7 152 
as = 12-+—& + ———¢' - - -, (15) 
3°5 34-53-11 
6 2 | 
n= Pull) +] —Pul)~——Pal 0 
52-7 33-7 
+ 4 
-|= Paul) — Plt) 
54-7 3°-5-7-1. 
_ ———— P7,(£) }e#+ (16) 
2.7-112-13 





and for the fifth harmonic 


19 
as = 30+——€+ ---, (17) 
3-13 
| 5 
us = P3\(&) +) —Pail€) 
| 33-11 
15 
————P;,(t)te&+---. (18) 
11-132 


SOLUTION OF THE EQUATION IN 7 FOR LARGI 
VALUES OF THE VARIABLE 


Since we are interested in a diverging wave 
system the boundary condition at infinity de- 
mands a solution of (10) of the form 


Y(n) =v(n)e*, (19) 


where v() is a series in descending powers of ». 
If we put s=7/en the differential equation for v 
becomes 

d*v dv a-é 

2°— +-2(2+1)——-———_1=0. (20) 

dz? dz 1+é&2:* 
We expand the denominator of the last term by 
the binomial theorem and look for a solution 
of the form 


v=1+a\2+a92"+a32°+::-, 
in which we use the characteristic value a; 
already obtained for the particular harmonic 
under consideration. 


After some labor we obtain for the funda- 
mental the two way infinite series 


if 2 2 4 62 \! 
+ — ef + ——— ¢P-— ——_+ - 
el 5 587) 3-58-7586 79-11 I 
if 69 62 3943 \1 
+-|1- —e— — ¢4-+-———_—_——_ e+ 
st 52-7) 3-547 3-58-78-11 | 9? 
94 2014 34,588 )! 
yo [1 a 5 Ce 5 omneninen Wes Se 
Sel 52-7 32-58 3-55-73-11 | 9? 
289 75,914 | 1 
+—|1- saansteilddlices a Fr ies 
5-71 38-52 38-54-72. 11 |! 


(Formula continued on next page.) 
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3 ¢ 46 761,914 ) 1 
+- {1 —¢°+- —_———¢! -- 
5.7 « 3.52 33.54.72.4] } 3 
1 371 \1 
iid | 
3-7 3-52-11 J n° 
1 i 182 1 
+- om {1 - e+:-: = 
3-7 «€ 52-11 | a 
1 |! 
+ {1+ - 
3-11 }n 
1 72 |! 
—— tite, too] (21) 
3-11 « } 


The calculation of this series gives us an 
absolute check on the correctness of the previ- 
ously computed characteristic value (13), for 
the successive terms in the series for a are of 
just the magnitude necessary to cause each 
alternate coefficient of z to start with a power of 
« two greater than the previous coefficient. 
Shortly we shall obtain an absolute check on the 
correctness of every coefficient in (21). 

Similarly, we obtain for the third harmonic 


Y3(n) =e 





1 2 38 |! 
+6-|1- —F4+-——_¢ - 
el 32-5 38-53-11 be 
1 22 398 | 1 
-15}1-- ftp ——' + 
é 32.52 38.58 Jy? 
1 8 3242 | 1 
-15|1- #4 
el 52 39-58-11 | 93 
1 1667 1 
-10—]1-—— -¢ 
é 2-35-5*-11 n* 
1 2521 1 
-10-41- —é 
Pa 33 . 5°. 1 1 | n° 
75 1 | 1 
meh | 
11 € J n° 
75 1 1 
-— —|1+ | + (22) 
ll é n' 





and for the fifth harmonic 


1 1 1 1 
Vata) =e +15 (14---}-—105—{14---} 
€ n € n° 
1 1 1 1 
—420—}1+ }_4+945—{1+.---! 
Pa n e? n* 
1 1 
+945—{1+---}—+ | (23) 
Pau n° 


to more than the needed degree of approximation. 

While (21) serves for the calculaton of the 
frequency and decrement of the fundamental 
free oscillation for eccentricities less than 0.8, 
the series converge altogether too slowly to be 
useful for the important cases of eccentricities 
nearly equal to unity. To treat the latter we 
must obtain solutions of (10) in the neighborhood 
of the origin (n=1). 


SOLUTION OF THE EQUATION IN 7 FOR SMALL 
VALUES OF THE VARIABLE 


If we put t=7°—1 the differential equation 
(10) becomes 
41(1+2)d*?YV/d?P +2 Y/dt 
+f&i+i—a}Y=0. (24) 
As this equation is of the second order the 
complete solution is a linear combination of two 
independent primitives. We look, therefore, for 
series solutions in powers of é the coefficients 
of which are polynomials in positive powers of f. 
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Denoting the two independent primitives by U 
and I’, we find for the fundamental (a=a;), 


1 1 1 
U,=t—-—f'é-—{ —-__f- rye 
2-5 2-53-7 =. 28-5-7 


(1+7)?—1 4 1 
+(140 | 2+( — ie 
(1+72)?+1 5 5§ 


284 151 1 
+(— 14+ es 
53-7 58-7 22-5+7 


8632 115,646 667 
+( 3 


),= U; log 


t+——_F 


3-597 34-59-7 — 34-58-7 


1 
inh Sininaiiiads ret _ | (26) 
23. 33-5+7 


The most general solution of (10) is, therefore, 
Y,=A,U,+B,!; (27) 


in the case of the fundamental, where the 
constants A; and B,; may be functions of the 
parameter ¢e. Now, in order to satisfy the 
boundary condition at infinity, (27) must be 
identical with (21). Therefore we expand the 
logarithm in a power series and replace ¢ by 
n°?—1 in the former, and expand the exponential 
in a power series in the latter, and compare 
coefficients. In this way we find that 





A\,= — fea, B,=3(it/e)hi, (28) 
where 
1 187 26,021 
a,=1———-e+—_-¢ --- meen + © 2 
2-52 23. 54-7? 24. 34-5°-7? 


= 1—0.020,000e +0.000, 763 
—0.000,026e°+---, (29) 
19 2609 32,593 


by = 1 ——_¢ —-——_—__¢+ ——¢6-+ 
: 24. 34.55.72 





= 1 —0.380,000e —0.010,649¢! 
+0.000,164e°+-+--. (30) 
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In carrying through this process each term in 
A, and B, was calculated by comparing at least 
two different groups of terms in the two solutions, 
and every term in each solution was used at 
least once. Therefore the calculation of A, and 
B, constituted an absolute check on the accuracy 
of every numerical coefficient in (21) on the 
one hand and in (25) and (26) on the other. 

Similarly in the case of the third harmonic 


(a=az3) we find 


5 I 5 
U;=t+ e-(. e+-— ye. , (31) 
y 3-5 23 3* 


1 
+4044 f 
3 2 


vi 


(1+/)'—1 
V3= U; log - - 
(1+2)?+1 


2 11 5 
+( - t—— eer] (32) 
38-5 2-38-52? 3? 


and the solution satisfying the boundary condi- 
tion at infinity is 


Y;=A3U3;+B3V3, (33) 
where 
1575 2 
A 2s—— -€'3, B;= one bs, (34) 
3-5?-7 2 @& 
and ad3=1+---, (35) 
37 
b;=1- . e+-:-: (36) 
2-33.52 


In the case of the fifth harmonic (a=a;5) 


7 «621 
Us=t+-f4+—! +--+, (37) 
2 ? 
(1+2)'—1 
V;= U; log —————— 
(1+¢)'+1 


; te 
+1404. +-t+- ete (38) 
a$$ 2 F 
and the solution satisfying the boundary con- 
dition at infinity is 


Y;=A;U;+B;V;, (39) 
where 
3,274,425 7 
A;=negligible, B;=—-—— b;, (40) 
25 @ 
and b;=1+->>:, (41) 
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FUNDAMENTAL FREE OSCILLATIONS 


The boundary condition at the surface of the 
conductor for free oscillation is that E:=0 for 
all ¢’s. Therefore it follows from (4) that 


dY/dn=0, n= No, (42) 


where 1, mo is the eccentricity of the conducting 
surface. 

For large values of m we use (21) for the 
fundamental. Performing the differentiation and 
solving for iw v by successive approximations 


we find 


w lf 53 1 3764 1 
j-=—i1-- sa merinee 
v 2a| 875 not 65,625 y' 
8,242,450 1 
176,859,375 ms 
(3)5) 41 971 1 
+i—)14+——+— — 
2a 15 no* 7875 not 
40,412 1 203,351,122 1 
we ead ditllpennndenans et as oa 
590,625 no 4,775,203,125 m8 
1 1 1 
=—! 1 —0.0606— — 0.0573 
a | No" No 
1 | (3)? f 1 
—0.0466——+---$+7 !1+0.2667 
No 2a No 
| 1 
+(0.1233—+0.0684 
no" no" 


1 
+().0426—+-:-- | . (43) 
No 


where a is the semi-major axis of the spheroid. 


TABLE I. Wave-lengths and logarithmic decrements for 
fundamental free oscillation of spheroids of 
eccentricities between O and 0.8. 





1/no \ ‘4a é 

0.0 1.814 3.628 
0.1 1.809 3.618 
0.2 1.794 3.588 
0.3 1.770 3.537 
0.4 1.734 3.461 
0.5 1.686 3.356 
0.6 1.625 3.214 
0.7 1.549 3.024 
0.8 1.455 2.772 
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(b) 


Fic. 2. (a) Wave-length \ and (b) logarithmic decrement 
6 of fundamental free oscillation of spheroid of eccentricity 
No- 


If we write 

lw v= x+y, (44) 
the ratio of the half-wave-length \ 2 to the 
major axis 2a of the spheroid is 


\ 4a=r 2ay (45) 
and the logarithmic decrement 64 is 
6=2rx/¥. (46) 


Unfortunately the series (43) converge too 
slowly to permit accurate calculations for eccen- 
tricities much greater than 0.7. In Table I the 
quantities \ 4a and 6 are given for eccentricities 
from 0 to 0.8. The calculated values are plotted 
on the graphs of Fig. 2, the portions of the 
curves corresponding to values of the eccentricity 
near to unity being plotted from the data in 
Table II, which will be computed later. We 
note that both the wave-length and the decre- 
ment decrease with increasing eccentricity. 

Since (21) does not converge rapidly enough 
to describe the fundamental free oscillations for 
eccentricities close to unity, we must resort to 
the solution (27). We shall consider only eccen- 
tricities so close to unity that f)=m?—1 and its 
positive powers are negligible. If then we put 


no+1 . 
j= tos — -2| (47) 
no— | 








ij-~— 
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7 
4a 
1.020- 
1LO10F 
6. To” 10° ig? 
1.000 T — + 1 — a 
0 0.050 0.100 | 
(a) 
- 
0.50 
025- ‘ 
fo" 16° ww? io? & 
0 = 4 : + * 
0 0.050 0.100 | 
(b) 


hic. 3. (a) Wave-length \ and (b) logarithmic decrement 
6 of fundamental free oscillation of spheroid of major axis 
a and minor axis }. 


and write 


1 9 886 
my=1- —fo— ——- 


5 53-7 3455-7 


e]4---- 





= 1 —0.200,000e/ +0.010,286¢4/ 
—0.000,500e7+---, (48) 
the boundary condition (42) gives 
b, = (4/9)téla,/ my. (49) 


We shall compute the wave-length and decre- 
ment for six values of / corresponding to values 
of ratios of the minor axis 26 to the major axis 
2a of the spheroid extending from 0 to about 
1/74. 

First consider the limiting case (/=0) of 
eccentricity unity, the ratio of b to a being zero. 
The boundary condition (49) reduces to },=0. 
Therefore ¢ is entirely real and the decrement 
is zero. This does not mean that no energy is 
radiated from the spheroid, but rather that the 
energy stored in the electromagnetic field is 
infinite compared with the energy radiated 
during a period of oscillation. Solving for e 
we find 

e= 1.5708 (50) 


correct to five significant figures. This is exactly 
7 2. Hence, as a=f in this case, the half-wave- 
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length is exactly equal to the major axis of the 
spheroid. These values of both wave-length and 
decrement agree with those obtained by Abra- 
ham for the limiting case. 

To discuss the remaining cases we multiply 
(49) by m,/a;. Combining the series represented 
by a,, 6; and m,, we get the formula 


1— |0.360,000 +0.200,0007! & 
— }0.018,612 —0.082,286/} e 
+ }0.000,093 —0.000,4817! €° 
=(0.444,4447/6. (51) 


We know that ¢ cannot differ greatly from 
7/2 over the range considered. Hence we put 


e=(27/2)(1—g-—17h). (52) 


A few trials indicate that an accuracy of one- 
tenth of one percent or better in the wave-length 
and better than one percent in the decrement is 
obtained if we retain terms in g, h, gh, h®? and 
neglect all higher powers. Then (51) leads to 
the simultaneous equations 


i) 


.2214—0.97351)g+ (1.5472 —2.4039/) kh? 


( ) 
=5.1677](h—2gh), (53) 


(2.2214—0.9735/)h — (3.0944 —4.8078/) gh 
= 1.72261(1—3g—3h?), (54) 


which we may solve for g and h/ for any given 
value of / in the range considered. In this range 
it is unnecessary to distinguish between a and f. 

In Table II are entered the results of the 
computation. The second column contains the 
ratio of b to a for the assumed value of /. These 
results are plotted in Fig. 3. 


FORCED OSCILLATIONS 


We now turn our attention to oscillations 
forced by an electric field E= Eye~*“* parallel to 
the x axis. As the assumption of a uniform field 
TABLE II. Wave-lengths and logarithmic decrements for 


fundamental free oscillation of spheroids of 
eccentricities near to unity. 





l b/a d/4a a 
0.000 0.00 1.000 | 0.000 
0.020 1.02 (10)-™ 1.001 0.098 
0.050 | 3.34 (10)-* | 1.004 | 0.247 
0.080 1.42 (10)-? | 1.009 0.396 
0.100 | 4.96 (10)- | 1.015 0.494 
0.125 | 1.35 (10)? 1.023 0.613 
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over the region occupied by the conducting 
spheroid is a valid approximation only when the 
wave-length is long compared with the shorter 
axis of the spheroid, we are limited in our 
discussion of oscillations in the neighborhood of 
resonance to eccentricities close to unity. This, 
however, is the range of most interest. 

If y is the angle which the unit vector &, 
(Fig. 1) makes with the x axis, we find 


1-#\3 
cos y=nf - :). (55) 
ar? 


Therefore the boundary condition at the surface 
of the conductor is 


1-£ 
E:+E n( — - :) =(. (56) 


no" — & 
With the use of (4) for E: this condition becomes 
(0A /An)o=i(xk/u)'ea(1—#)Epe**', = (57) 


where we have replaced fn by its equal a. 

Evidently A must be a linear combination of 
the solutions representing the odd harmonics, 
that is, 


A = { Ciu,(&) Yi(n) + C3us(€) Ys(n) 


+ Csus(&) Yan) +> +} (1B) 'e-, 
where 1;(£), u3(¢), u5(£) are given by (14), (16), 
(18), respectively, and Y,(n), Ys(n), Ys(n) by 
(27), (33), (39). 

Determining the coefficients Ci, Cs, 
the boundary condition (57) we have 


a=i(* “) ak, 1-——¢ 
54-7 


(58) 


C; from 











dy, 
Fg dE), 0 
358-7 
K\? 1 
c= i(- ak, € 
m 3°57 
2 dY;3 
= ete} /(—), (60) 
3?- 54 dn /o 
x\} 1 dY; 
c.=i(*) ua, ef+--- I/(—) (61) 
m 5-7 dn /o 
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To save writing we shall put 
k=2nc(x/p) ako. (62) 


Then the current J) at the center (€=0) of the 


antenna is 


Iyp=2nc(A) <0, 





n=n0, 
2 
= itd | 1- é 
54.7 
8 Y;(no) 
+— e°+ (0) ——— 
3-58. (d Y;/dn)o 
1 2 Y3(no) 
+—"|1- ets tuo) —__—_— 
3-5? 3- 5° (dY3/dn)o 
1 Y;(no) 
 seroramnenen a Le he -+1u,(0)——— 
34-5 » 7 (d Y;, dn) o 
ions ‘[ iwt (63) 
where 
1 89 733 
u,(0) = 1+——e ——_——-e ++ e+ : : -," 
2-5° 23-54-7? 24. 34- 55-7? (64) 
3] 13 | 
u3(0) = —-; 1 -———€'+ :: +}, (65) 
2l 2-33.52 
15 
u;(0) =—{1-+4----}, (66) 


from (14), (16) and (18). We note that (0) 
the reciprocal of the series a; given by (29). 
From (27) we find 


37 
Vilno) => on (67) 


(5)--3 
to a sufficient degree of 


eccentricities so near to unity as those under 
consideration, where 


3ibym, 2 
-- a vi (68) 
3 


approximation for 


1 89 733 
e— e*+ —¢+ 
2.52 93.54.72 94. 34.55.72 
= u,(0) = 1+0.020,000€ 
—0.000,3636+0.000,004e8+---, (69) 








¢=1+ 











(63) 


(67) 


(68) 
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and from (33) 


525 12 
Y3(no) Te cea a Cx. (70) 
23 ras 
dY3 1575 7 bss s 
(—) =—— — ——_+4+—_—-¢4 (71) 
dn 0 2? é€ l 3°5*-7 
where 
23 
¢3;=1+————e+ ::-,"7 (72) 
2 . 33.52 
2 1 
m;=1—-l+——el+:--. (73) 
3. 2-3-5 


As the term in the fifth harmonic is quite 
negligible we shall not trouble to write it down. 

The second term in (71) is negligible compared 
with the first, and we need retain only the first 
term in the series for 63; and c3 and the first two 
terms in the series for m;. Let us put 


54-7 3-58-7 


= 1—0.000,4576+0.000,024e8+---, (74) 


In the series ss; we need retain only the first 
term. Thus we obtain for the current J) at the 
center of the antenna 


I (les\c;*) (4/e8a,/9) 
a. (4/€°a,/9)?+ (b,m,)* 


(les \¢;") (bym,) lé | 
-i{_ len a I ist, (76) 
(4/<a, 9)?+ (b\m,)? 300 


The third harmonic contributes only a very 
small term to the imaginary part of the current. 
The series a;, b;, m, and s; are specified by 
(29), (30), (48) and (74), respectively. The 
series for c,’, obtained by squaring (69), is 


c= 1+0.040,000e —0.000,327€4 
—0.000,007€°+---. (77) 


The square of the current amplitude Joo at 
the center is 
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(lesic;* )? 
fae] " 
(4/e8a ,, 9)? + (b,m)° 


€ 
| 1-. b(t +4) | (78) 
150 


5 
to a sufficient degree of approximation. To find 
the frequency of resonance we must equate to 
zero the derivative of this with respect to e, or, 
more conveniently, with respect to &. This 
leads to the formula 


1 —0.00560/ 


b, =0.14750/ ° (79) 


(1 —0.43837/)? 


where, in the very small right-hand member, 
we have replaced ¢ by 7/2 in the numerator and 
denominator of the fractional factor. The justi- 
fication for this lies in the fact that the values 
of « computed from (79) are found to differ very 
little from 2/2 in the range under consideration. 

We have computed the wave-length A and 
the current J) at resonance, and also the current 
amplitude Jo and the angle @ by which the 
current leads the electromotive force, for six 
values of 1. These are contained in Table III 
and are plotted in Fig. 4. In the second column 
of the table is given the ratio of the minor to 
the major axis of the spheroid corresponding to 
the assumed value of /. 

It will be noted that the real part of the 
current is substantially the same for all values 
of / considered, and that the current leads the 
electromotive force at resonance for all values 
of 1 greater than zero by an angle @ which 
increases with increasing 7. For /=0 the half- 
wave-length is exactly equal to the major axis of 
the spheroid as in the case of free oscillation, 
but for values of / greater than zero the resonant 
wave-length is slightly shorter than the wave- 
length of free oscillation. Experiments on an- 


TABLE III. Wave-lengths and currents at resonance. 





l | b/a | A/4a 1o/ke~ *# Too/k ¢ 
0.000 (0.00 1.0000 | 1.045 1.045 0°.0 
0.020 1.02 (10)-™) 1.0004 | 1.045—0.0287 | 1.045} 1°.5 
0.050 (3.34 (10)~*| 1.0027 | 1.045—0.0712 | 1.047 | 3°.9 
0.080 1.42 (10)-*} 1.007 1.045 —0.1132 | 1.050 | 6°.2 
0.100 /4.96 (10)-*} 1.011 1.045—0.1412 | 1.055 | 7°.7 
0.125 |1.35 (10)-*} 1.017 1.046—0.175z | 1.061 | 9°.5 
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Fic. 4. (a) Wave-length A, (b) current amplitude Je and 
(c) lead @ of current at resonance. 


tennas in the form of right circular cylinders 
indicate that the half-wave-length at resonance 
is from five to six percent longer than the 
antenna.‘ This is in satisfactory qualitative 
agreement with our calculated ratios of \/2 to 
2a for the prolate spheroid, for the diameter of 
a prolate spheroid falls to half of its maximum 
value at a distance x=0.866a from the center 
whatever the eccentricity may be, and conse- 
quently the prolate spheroid which best fits a 
right circular cylinder is one which is slightly 
wider at the center and somewhat longer,® as 
illustrated in Fig. 5. We may call this prolate 
spheroid the equivalent spheroid for the given 
cylinder. Nevertheless it is very desirable that 
experimental measurements of the ratio \/4a 
should be made for actual spheroids in order to 
corroborate the theory. 

We have computed resonance curves for/=0.02, 
0.05 and 0.10. The data for these are given in 
Tables IV to VI, where ¢, represents the value of 
¢ at resonance. The current amplitude Jo at the 
center of the antenna and the phase ¢ by which 
the current leads the electromotive force are 


4C. R. Englund, Bell System Tech. J. 7, 404 (1928). 
5 Rayleigh, Phil. Mag. 8, 105 (1904). 
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depicted in Fig. 6. It will be noted that the 
resonance is sharper and that the phase shift 
+7/2 
more rapid the smaller /. 


The distribution of current along the length 


toward as we depart from resonance is 


of a thin antenna at resonance has been meas- 


ured. Since £=x/fn=x/a, the quantity J;/J, 
represents the ratio of the current, at a point 
whose distance from the center of the antenna 
is ga, to the current at the center. The experi- 


mental measurements indicate that, at resonance, 
I:/Ip=cos $ré (80) 


within the experimental error. 

In developing an expression for this ratio from 
our theory, we obtain sufficient accuracy if we 
neglect the third and higher harmonics. Then 
we get the formula 


é! 8 e® 
+——4 1 - ——¢ - ————_ f+: :- }. (81) 
23-5-7 3?-5? 24. 338.5-7 


Let p be the ratio of the frequency for which 
it is desired to evaluate (81) to the frequency 


inaction oma” 


Fic. 5. Prolate spheroid equivalent to cylinder. 








for resonance with /=0. Then 


p=e, 1.5708, (82) 
and (81) becomes 
T:/Ig=(1—£) [1 —0.24674p* | 1 —0.05640p° 
+0.00371 pt} 2 +0.02174p*| 1 —0.08773p"} = 
—0.00099p%E®+---]. (83) 
TABLE IV. Current for 1=0.02, b/a=1.02 (10)-™™ 
ep = 1.5701. 
cler Io/ke~*#t Took ° 
0.875 0.011 —0.120: 0.120 S4°.8 
0.950 0.081 —0.292: 0.303 74°.6 
0.975 0.275 —0.474: 0.548 59°.9 
0.990 0.726—0.495; 0.879 34°.3 
1.000 1.045 —0.028: 1.045 pe 
1.010 0.757 +0.453: 0.882 — 30°.9 
1.025 0.313+0.464: 0.500 — 56°.0 
1.050 0.107 +-0.302: 0.320 —70°.5 
1.229 0.022+0.1377 0.139 —80°.7 
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When p=1 this reduces to 
Ie/Ip= (1 —&) [1 — 0.233782 +0.0198¢! 
—0.0010¢°+---], (84) 


which is exactly cos ré/2 to the fourth decimal 
place. When /=0 it should be remembered that 
the term in the third harmonic disappears from 
(76) and therefore, in that case, we make no 
error in neglecting it. 

In Table VII are listed the values of J:/Jy 
over the range of & for various values of p near 
to resonance. 

Next we shall calculate the rate of absorption 
of energy and the radiation resistance of the 
antenna. Since the element of distance ds at the 
surface of the spheroid in the direction of 


increasing & is 


‘ ne —& . 
ds =f( ——— } dz, (85) 
i-¢# 


the electromotive force d& corresponding to the 
increment dé is 


dé = Eye ‘*' cos yds = Eye ‘*'adé (86) 


from (55). 

At any & the ratio of the part J,’ of the 
current in phase with the electromotive force to 
the part Jo’ in phase with the electromotive force 
at the center of the antenna is 


TI’ (1—-&))uy(&) (1-2) my (8) 
= — (87) 
Py u,(O) Cy 


from (58). 
Hence we obtain from (76) for the mean rate 
of absorption of energy. 


TABLE V. Current for 1=0.05, b/a = 3.34 (10)~, 


€-= 1.5606. 
€/€r To ke~*#t | Too/k | Ce) 
0.850 | 0.042 —0.2387 0.242 80°.0 
0.900 | 0.104 —0.347: 0.363 73°.3 
0.950 0.352 —0.529% 0.635 56°.4 
0.975 | 0.707 —0.525i 0.880 36°.6 
1.000 | 1.045 —0.071i 1.047 | 3°.9 
1.025 0.783 +0.4187 0.888 —28°.1 
1.050 | 0.450+0.482: 0.660 —46°.9 
1.100 | 0.182+0.360i 0.403 —63°.2 
1.150 | 0.100+0.270 0.288 —69°.7 
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Fic. 6. (a) Current amplitude Jg@ and (b) lead @ of 
current plotted against ratio ¢e/e, of frequency to frequency 
at resonance. 


dU (les ,c;") (4l ea 9) a 
dts wn (4le8a, /9)?+(bym,)? cy 
at 
x | (1—#)!ui(E)dé. (88) 
e l 
Now 
41 aa 
(1—#)'n()de= | Pi,(t)uy(t)dt=4/3 (89) 
. 1 7} 


from (14). Consequently, from the relation 
a,¢,;=1 and from (62), 


dU 64 K i a*f? a;! 
dt/y 27 7 > Sy 


(les\c;")? 


x - - (90) 
(4]e%a, /9)?+ (b\m,)? 
TABLE VI. Current for 1=0.10, bja=4.96 (10)-4, 
€p-= 1.5539. 
€/€r 1o/ke-**! Too/k Prey 
0.850 0.144—0.428: 0.451 71°.8 
0.900 0.306 —0.546: 0.625 60°.8 
0.950 0.669 —0.575: 0.882 40°.7 
1.000 1.045—0.1417 1.055 pe | 
1.050 | 0.822 +-0.3601 0.897 — 23°.6 
1.100 | 0.504 +0.450: 0.676 —41°.8 
1.150 | 0.326+0.4102 0.524 —51°.5 








830 L. PAGE 


in terms of the amplitude Ey of the exciting 


field, or 
‘© 16 (“) Too? f? ay! 
x). 27 cS, 
xf 1+ b,(1+3 | (91) 


in terms of the amplitude Jo of the current at 
the center of the antenna. 

The radiation resistance R, defined as the 
quotient of the mean rate of absorption of 


energy by 3J?, is 


32 1 f? 
R=—*(~) — “14+ b(1+ | (92) 
27 \k cn? Sy 150 


In practical units, the radiation resistance R, 


in ohms is 


4a\? a;4 
R,=87. 670(“ y( ) - 
K nN Sj 
e 
x] 1+" (1430 | (93) 
150 


provided we replace f by a. In Table VIII 
given the radiation resistance at resonance for 
the six values of / considered before. 

In Table LX we give the quantity («/)'R, for 
various fractions ¢/e, of the frequency at reso- 
nance, both for /=0.050 and for /=0.100. The 
results are plotted in Fig. 7. 

Finally we turn our attention to the electro- 
magnetic field of the oscillating antenna. The 
three nonvanishing field components £;, E, and 
H, are given in terms of A and its derivatives 
by (4), (5) and (6), respectively. Since the 
coefficient C; in (58) may be expressed in terms 


TABLE VII. Current distribution I;/Io. 





0.95 | 1.00 | 1.05 | 1.10 


+0.1 | 0.988 | 0.988 | 0.9877| 0.987 | 0.987 
+0.2 0.953 | 0.952 | 0.9511 | 0,950 | 0.949 
+0.3 0.894 | 0.893 | 0.8910] 0.889 | 0.887 
+04 | 0.815 | 0.812 | 0.8090} 0.806 | 0.803 
+05 | 0.715 | 0.711 | 0.7071 | 0.703 | 0.699 
+0.6 | 0.597 | 0.592 | 0.5878] 0.583 | 0.578 
+0.7 | 0.464 | 0.459 | 0.4540} 0.449 | 0.443 
+08 | 0.318 | 0.313 | 0.3090} 0.304 | 0.300 
+09 | 0.162 | 0.159 | 0.1564] 0.153 | 0.150 
| 
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of C, by the relation 

C:= — (2b, 37+ 54-7) eC, (94) 
to a quite sufficient degree of accuracy, we find 
that 


ik 4le8a,/9—ib\m, 
A = ——leés,— 
3mc (4/e*a, 9)? + ( bim,)? 
2b; | 
ui() Yi(m) — ——e'u3(£) V3(m) ¢ (1 — &) te f* 
| 32.58. (95) 


At a distance from the nearest point on the 
antenna not less than the semi-interfocal distance 
f the field components are specified to better 


than one percent by using the functions (21) 
and (22) for Yi(m) and Y3(n), respectively, 


although very close to the antenna it is necessary 


to use (27) and (33). 
The components of the radiation field are 


k M i]es, 1 
wat (eye! 
3reX\kK f (f-#) 


1 
{ (4le%a,/9)? + (bym)*} ? 


2b,é4 | 
x | m(8)- : u3(E) pe'er-e'trie-y, (96) 
32.547 | 
k les, 1 
IT,=—- 
3mc f (n?—1)! 
1 
x a ; i 
{ (48a, /9)?+ (bym)*} 
2b,é4 
xt ef) -—— unt) eiten—ottr/2—y) (97) 
54 
bim, 
where tan y= —, (98) 
4/le*a, 7) 


TABLE VIII. Radiation resistance at resonance. 


l (x u)tRp ohms) 
0.000 73.1 
0.020 73.0 
0.050 72.8 
0.080 72.3 
0.100 71.8 
0.125 71.1 
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e find 
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For small / it follows that y is nearly zero at 
resonance, nearly 7 2 below resonance, and 
nearly — 7 2 above resonance. 

As the Poynting flux is —cE;:H, in the 
direction of increasing 7 and the element of area 
of a spheroidal surface is 2 xf? | (n? — &)(m?—1)} 'dé, 
the mean rate of radiation of energy is 


dR 16 K\? a*f* 
( ) = ( ) cEe—a\' 
dt/» 9 im “ 


(Jes \c;")? 


af 
x | [ui(é) Pdé, (99) 


(4/e%a,//9)? + (bym,)*7 _) 
provided we make use of the relation a,c,;=1 
and neglect, as heretofore, the square of the 
third harmonic. From (14) we find 


i 
| [u1(t) Fade 
e l 


4 2 8 4 
| ef+---b=—, (100) 


S4-7 = (358-7 35) 


TABLE IX. Radiation resistance (x/y)*R, in ohms for 
frequencies near resonance. 


¢/¢r 
| | 
l | 0.850 0.900 0.950 | 1.000 1.050 1.100 1.150 
0.050 | 55.3 | 61.0 | 66.9 | 72.8 | 78.7 | 84.6 | 90.6 
0.100 | 54.5 | 60.2 | 66.0 | 71.8 | 77.7 | 83.6 | 89.4 


| 
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Fic. 7. Radiation resistance Ry, in ohms plotted against 
ratio €/e, of frequency to frequency at resonance. 


where s; is the series defined by (74). Therefore 


dR 64 K\? a*f* a, 
( ) 7 r° c Ey 
dt Ay 27 MK te Sj 


(lesyc;*)* 
x (101) 


(4]e%a 9)*+ (b,m1,)* 


in agreement with (90). 
In a future communication we hope to discuss 
the field in more detail. 


Note added in proof. Since submitting this article for 
publication there has come to the attention of the authors 
a recent paper by L. V. King [Trans. Roy. Soc. London, 
236, 381 (1937) ] in which the radiation from a cylindrical 
antenna is investigated by a very different method from 
that pursued here. 











MAY 15, 1938 


PHYSICAL 


REVIEW VOLUME 53 


Statistical Theory of a Scaling Circuit 
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The over-all efficiency of a counter circuit embodying a 
scaling circuit to reduce the rate of recording is investi- 
gated. The probability of detecting counts (G(n,r)) of a 
scaling circuit and recorder fed with random impulses is 
shown to be: 1—J(ur, n), where 7 is the scale of the circuit, 
7 the resolving time of the recorder, u the statistical average 
of incoming counts, and where J(ur, m) is the ratio of the 
incomplete gamma-function with upper limit wr of n 
to the complete gamma-function of n. Efficiency curves 
(100 G(n, r)) are plotted showing that the distribution in 
time of the impulses leaving the scaling circuit is no longer 


random, but is a distribution in which the extreme short 
and long intervals between impulses are averaged out. The 
average rate of recording counts (u,) is (u/n)G(n, r). The 
analysis is extended to include the effect of the finite 
resolving time of the amplifier (¢), vielding when 7 2no: 
un = (u/n)[1/(uo +1) ][1 — 7 (u(r — no), n) ],and when r Sno: 
un = (u/n)[1/(uo+1)], the latter being the response of the 
amplifier alone. Thus, one need only make n>+r/o and the 
scaling circuit and recorder will operate with 100 percent 
efficiency. 





I. INTRODUCTION 


COUNTER-amplifier-recorder circuit is 

limited in the number of impulses per 
second it can detect and record mainly by the 
mechanical recorder which usually has a re- 
covery time much greater than the recovery 
time of any other part of the circuit. If an impulse 
to the recorder is preceded by an interval of 
time ¢ from the last count which is less than the 
recovery time 7 of the recorder, then the latter 
fails to advance one count. 

In order to record counts at a faster rate a 
scaling circuit is usually introduced which sends 
only every mth impulse into the recorder. The 
fact that a scale-of-m counting circuit has an 
advantage other than simply dividing the count- 
ing rate by » has been noted by Stevenson and 
Getting.! They stated that the distribution in 
regard to time intervals of the impulses leaving 
the scale-of-” circuit was no longer random, but 
had a standard deviation less than random 
counts coming in at the same average rate. 
This means that there are fewer short intervals, 
and therefore the mechanical recorder will not 
miss so many impulses that come to it. It is 
desired here to derive an expression for the over- 
all response of a counter-amplifier and scale-of-n 
recorder set in order that limits of usefulness of 
the scaling circuit may be evaluated. 


1 E. C. Stevenson and I. A. Getting, Rev. Sci. Inst. 8, 414 
(1937). 


II. RESPONSE OF PERFECT SCALING CIRCUIT AND 
RECORDER FED WITH RANDOM COUNTS 


It is of interest to investigate first the response 
of the scaling circuit and recorder alone. There- 
fore we assume that the counter, amplifier, and 
scaling circuits are perfect, and that only the 
recorder has a resolving time 7 not zero. Let 
Q(t)dt be the probability that the interval of 
time between a count and its predecessor be 
greater than ¢ and less than ¢+dt. Let L(n, #), 
G(n, t) be the probability that the time elapsed 
when m successive counts follow after an arbi- 
trary count be less (greater) than ¢. Then always 
L(n, t)+G(n, t)=1. 

The probability that »—1 counts follow after a 
first count in time less than ‘—x is L(n—1, t—x); 
the probability that the mth count follow after 
the last of these in time x to x+dx is Q(x)dx. 
Hence the probability that the » counts come 
together in this fashion is 


L(n—1, t—x)Q(x)dx. (1) 


The probability that the combined n-interval be 
less than ¢ is therefore :* 


? The functions L(n, ¢) satisfy a simple recursion relation: 
: sd " 
L(n+m, t)= I, L’(n, x)L(m, t—x)dx, (2c) 


where ” and m are integers. This relation can be proved by 
an induction on +m, assuming the formula holds for L(m, t) 
and L(m+n-—1, t) the integral on the-right of Eq. (2) can 
be reduced to that of the definition of L(n+m, t). From 
this it follows that the output of a scale-of-2* made up of 
k stages of scales-of-two in series is that of a scale-of-2* in 
one stage. The definition of L(2, ¢) shows that this is true for 
a single stage of two. Suppose it is true for the output of 
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Fic. 1. The efficiency (one hundred times the probability 
of recording count = 100 G(m, r)) of a scaling circuit and 
mechanical recorder as a function of the statistically 
averaged input frequency (u) of random counts (solid lines) 
multiplied by the resolving time (7) of the recorder. The 
numbers on the solid lines refer to the scaling ratios. In an 
actual case the scaling circuit is not fed with random 
counts from the amplifier, but with a distribution in which 
the short intervals are removed, increasing the efficiency 
of the scaling circuit alone. This is shown by the dotted 
curve (8A, for a scale-of-eight), for which 7 and o are taken 
to be 0.0133 sec. and 0.0005 sec., respectively. If the scale 
of the circuit were made greater than r/o (=26.6), the 
efficiency curve would move off to the right to infinity, and 
the circuit would detect all counts coming out of the 
amplifier. 


at 


L(n,t)= | QO(x)L(n—1,t—x)dx, (2a) 
Jo 
when L(1,#)= | O(x)dx. (2b) 


If now we take Q(t) to be that for random 
counting,® Q(t)=ye“', where uw is the average 
rate of receiving impulses, by an induction on n, 


a 
a | ue**(ux)""'dx (3a) 
0 


(n—1)!« 


L(n, t)= 


“ak (ux) "te-*d(ux). (3b) 
ian 


This is the same expression as that obtained by 
using a Poisson distribution of intervals. By 
definition of the I'-function, 


L(n, t)=T (ut, n)/T(n) =I (ut, n), (3c) 
k stages. Then the probability that an impulse in the out- 
put of the &th stage follow its predecessor in time ‘— x or less 
is L(2*,—x). The probability that the next impulse come 
in time x to x+dx is L’(2*, x)dx. Hence the probability that 
these two impulses come as specified is 


L'(2*, x) L(2*, t—x)dx. (2d) 

The probability that the two impulses come together in 
total time less than ¢ is 

fre, x) L(2*, t—x)dx =L(2*4, t) (2e) 


by Eq. (2c), and the induction is completed. 
"3 For proof of this see Marsden and Barratt, Proc. Phys. 
Soc. 23, 367 (1911) and 24, 50 (1911). 
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Fic. 2. The same functions (one hundred times proba- 
bility of recording counts) that are plotted in Fig. 1 are 
plotted here as functions of ur/n, i.e., as functions of 
equal average rates of recorder input. For a scale of one the 
efficiency is that for random counts fed into the recorder; 
for higher scales the distribution is not random. In these 
cases, the short intervals are averaged out if w<m/r, in- 
creasing the efficiency; for u>n/r the long intervals are 
averaged out, decreasing the efficiency in comparison to 
that for random counts. The curves show that this effect 
becomes more pronounced as the scale of the circuit 
increases. 


where I(ut, 2) is the incomplete I'-function with 
upper limit uf of , (7) the complete I'-function, 
and J(ut, m) this ratio. J(ul, m) is given in tables 
of the incomplete [I-function.4 The function 
G(n, t) has the form 


G(n, t)=1—L(n, t)=1—TI (ul, n). (3d) 


On carrying out the indicated integration in 
Eq. (3a) one can obtain by means of a recursion 
formula :° 


(ut)* 
L(n, t)=1 —e(1 +nh+—— 


2! 
(ut)"— 
ess Qemeen - ) (4a) 
(n—1)! 
or 
(ut)? (ut)"=' 
G(n, t)=e m(t+urt ——++-++-+——. -). (4b) 
2! (n—1)! 


G(n, t) is equal to e~“‘ times the first » terms of 
the expansion of e+‘ as a power series in /. 

If we now consider the impulses being fed into 
a recorder which will not respond unless the 


* Karl Pearson, editor, Tables of the Incomplete T- Function 
(H. M. Stationery Office, London, 1922). 

5B. O. Peirce, A Short Table of Integrals, No. 403 (Ginn 
and Company, Boston, 1929). 
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interval between the count and its predecessor is 
greater than 7, the fraction of the impulses 
recorded will be G(n, r). The average number, 
Mn, recorded per second will then be 


n=G(n, 7)/tn, (5a) 


where ?, is the average interval between counts 
in the output of the scale: 


in= | tL’(n, t)dt 


“0 
1 ie 2) 
- | tue-*"(ut)""'dt=n/p. (Sb) 
(n—1)!. 0 


Mn=(u/n)(1—L(ur, n)). (Sc) 


Graphs of the efficiency of the scaling circuit 
(100 G(n, r)) against wr as the argument are 
shown in Fig. 1 for values of » from one to 
thirty-two for scales-of-two in series. In order to 
compare the efficiencies of the various scaling 
circuits with equal numbers of counts per second 
going in to the mechanical recorder, 100 G(n, r) 
is also plotted against ur/n in Fig. 2. The dis- 
tribution of impulses according to the time 
interval between counts is the negative deriva- 
tive of G(n, t) with respect to time, and one can 
see from an inspection of these graphs that a 
scaling circuit, in averaging the time intervals n 
at a time, makes fewer extreme short and long 
intervals to be passed on to the recorder. Because 
there are fewer short intervals for a given rate 
of recorder input when a scaling circuit is used, 
it is better to increase the efficiency of recording 
counts by employing a scaling circuit than to 
attempt to build a mechanical recorder with a 
smaller time-constant and operate it directly out 
of the amplifier. 

The higher the scale of the circuit, the more 
the short intervals are averaged out. In Eq. (4b) 
it can be seen that as m increases, more and more 
of the terms of the expansion of e**" are included 
so that un, approaches the value (for each ur) 


(5d) 





Hence 


bn =(u/n)e*tet4#* = p/n, n>>urT. 


That is, the scaling circuit and recorder will 
register all of the incoming counts divided by n. 
In the next section we shall see that when the 
effect of the losses in the amplifier is considered, 
n will not have to increase without bound, but 


nN. M. SMITH, FR. 
for a finite and practical value of m all the im- 
pulses (divided by m) received from the amplifier 


will be registered. 


Ill. THe Errect DUE TO THE FINITE RESOLVING 
TIME OF THE AMPLIFIER AND COUNTER 


In an actual experimental set-up, the scaling 
circuit is not fed with random impulses, but with 
a distribution in which the short intervals have 
been cut out by the counter and amplifier. We 
assume that the counter and amplifier act to- 
gether with an effective resolving time, ¢.° Under 


such conditions Ruark and Brammer’ have 
shown that 
mi=pe*T-) /(ua+1), 720 (6a) 
Mi=pu/(uo+1), tTKo (6b) 
when no scaling circuit is used. 
With a scale-of-m introduced between the 


amplifier and recorder, we proceed to calculate 
the probability that an impulse in the output of 
the scaling circuit shall come within a time ¢ 
of its predecessor. The probability that an 
impulse in the amplifier output shall come 
within time ¢ and ¢+dt of its predecessor is 


Q*(t)dt, where 
Q*(t)=pe#(, 
=(), 


(7a) 
(7b) 


t2oa 
to 


Eq. (7b) is obvious. The amplifier cuts out all 
these counts. The probability that a count should 
come within a time o or more of its predecessor 
is e-**, This is the fraction of the total counts 
which makes up the output of the amplifier. 
The fraction of counts between ¢ and ¢+dt for 
t2o is e~“tdt. Hence the probability that the 
output impulse of the amplifier be between ¢ and 
t+dt of its predecessor is e~*'-%dt, the ratio of 
the two fractions. 

We also assume that if the scaling circuit is 
one of a single stage, that its recovery time does 


6 This time a is not necessarily the recovery time of the 
amplifier. The scaling circuit reacts if a variation of a cer- 
tain size occurs in the output of the amplifier. If impulses 
come too closely, this variation is not sufficiently large; it 
depends on the recovery time of the amplifier, the distortion 
caused by the amplifier, and the sensitivity (bias control) of 
the scaling circuit. 

7 These conditions and the following formulae are given 
by A. E. Ruark and F. E. Brammer, Phys. Rev. 52, 322 


(1937). 
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not exceed that of the amplifier; if the circuit is 
a series of scales-of-two, that the recovery time 
of the kth stage does not exceed 2*~'c. 

Now let L*(n, t), (G*(n, t)) be the probability 
that the time elapsed when m successive amplifier 
impulses follow an arbitrary one be less (greater) 
than ¢. 

Then, as before, 


t 


L*(2, i)= | O*(x)L*(n—1,t—x)dx, (8a) 


“0 
a 
where L*(1, ) =| O*(t)dt. (Sb) 
0 


These functions can also be expressed by means 
of the incomplete I’-functions: 


L*(n, t)=I(u(t—no), n), t2ne, (9a) 
=(, tSno. (9b) 


This can be shown by induction. For n=1, 
Eq. (9a) and (9b) hold. Let then Eq. (9a) and 
(9b) hold for m. Then: 


at 
L*(n+1, t)= | O*(x)L*(n, t—x)dx. 
“0 
By hypothesis, L*(m, t—x) is zero when (t— x) 
no; i.e. if x2 (t—no). Also Q*(x)=0 if x<o. 
Hence L*(n+1,#)=0 if (t—no)<o, ie. if 
t<(n+1)c. 
Let 2 (n+1)c. Then: 


at—no 
L*(n+1, t)= | pe #(2-9) 
XI (u(t—no—x), n)dx. (9c) 
Substitute x—o0=y, then 


t—(n+l)o 


L*(n+1, )= { pe~#y 
xX I(u(t—(n+1)o—y), n)dy (9d) 
= I(u(t—(n+1)e), n+1) (9e) 


and the induction is completed.* 


8 As before it can be shown that 
vt 
L*(n+m, t) = J, L*’(n, x) L*(m, t—x)dx. (9f) 


If the transformation y=x—no be made, this is reduced to 
the same form as that of Eq. (2c) in footnote 2. From this 
it can be shown that for a scale-of-2* made up of a series of 
k scales-of-two, that the recovery time of the /th stage can 
be made as great as 2'¢, and have L*’(2*, t) give the dis- 
tribution of intervals in the output of the scale-of-2*. 
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The function G*(n, ¢) has the form: 
G*(n, t)=1—J(u(t—no),n), t2no0 (10a) 
= 1, tSne. (10b) 


Plots of this function can be obtained by re- 
placing the abscissa of Fig. 1 by u(t—ne). It is 
seen that this probability is increased by the 
removal of the short intervals by the counter 
and amplifier. 

The dotted line of Fig. 1 has been drawn for 
a scale-of-eight for which values of r= 0.0133 sec. 
and «=0.0005 sec. are assumed. The abscissae 
are as indicated in the figure. 

As before, the rate of recording is 


pn* =G*(n, r)/t,*, (11a) 
where 
i,* = -| 1G*!(n, t)dt (11b) 
v0 
1 | 
—— | tue") u(t—na) dt (1c) 
(n—1)! ny 
=(n/u)(1+y0). (11d) 


Hence® 


Mn* =(u/n)[1—/(u(r—n0), n) |/(1+p0), 


t>mno_ (12d) 
ae 
SR ae eceeen, tone. (12b) 
ni+ypo 


For n=1 these are the formulae of Ruark and 
Brammer, for the function G*(n, 7) when not 
unity is e-“‘"-"? multiplied by the first ” terms 
of the expansion of et#(t~", 


IV. CONCLUSIONS 


For a_ scale-of-2*, the probability that an 
impulse in the output of the /th stage comes 
before a time ¢ from its predecessors is L*(2!, 2). 
This vanishes if t< 2'¢. Hence the resolving time 
of the (/+1)st stage can be made as great as 2'c. 
If for any stage the resolving time exceeds this 
figure, the formulae given do not hold, and some 
rather complicated expressions will enter into the 
formula for L*(2*, 2). 

® Note added in proof. For amplifiers not of the Neher- 


Harper type, un=(u/n)P(c)G*(n, r), where P(c) is the 
amplifier’s probability of counting. 
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We do not attempt to draw response curves for 
various values of +r and o, for these are best 
drawn for a particular apparatus. Instead we 
wish to point out a practical limitation to the 
scale of the circuit. In Eqs. (12) the term ex- 
pressing the losses of the scaling circuit and 
recorder alone drops out if t&no. In this case 
the recorder counts the full output of the ampli- 
fier; i.e., if the scale n2 r/o. Thus for r=0.0133 
sec. and ¢=0.0005 sec., 7/o=26.6, and a scale- 
of-32 would be sufficient; a higher scale would 
not increase the efficiency, for the losses due to 
the amplifier are independent of the scale. 

This analysis shows that a determination of 
the amplifier’s resolving time’ would be necessary 
no matter what the scale; for small values of yuo 





the error is equal to wo, for large uo it becomes 
quite appreciable, being about 1 
« =0.0005 sec. and u = 20 impulses per second and 
when »=200 impulses per 


percent for 
about 9 percent 
second. 
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Electromagnetic and Gravitational Radiation 
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Guided by the analogy of the electromagnetic field we consider here the influence of gravi- 
tational radiation on the gravitational field and on the equations of motion in the case of a two- 
body problem. This leads to a generalization of the previously developed methods of solving 








the gravitational equations and finding the equations of motion. 


HE “new approximation method” applied to 

the gravitational equations! leads to the 
formulation of the equations of motion of mass 
points represented by singularities of the gravi- 
tational field. In I it was claimed that this new 
approximation method in which no direction of 
time flow is distinguished differs from the old one, 
based on the linearized equations and on the use 
of retarded potentials. There seems to be no 
place for radiation in the new approximation 
method as developed in I. It is just this point 
which could arouse suspicion as to whether the 
new method of solving the field equations and 
formulating the equations of motion leads to 
results corresponding to reality. The purpose of 
this note is to discuss the point and to investigate 
1 Einstein, Infeld and Hoffmann. ‘‘The Gravitational 


Equations and the Problem of Motion,” Ann. Math. 39, 
65 (1938) quoted here as I. 


the connection between the new and the old 


approximation methods. 


§1. THe New AppROXIMATION METHOD AND 


THE EQUATIONS OF ELECTRODYNAMICS 
The equations of electrodynamics are well 
suited to show the essential features of the new 
approximation method. We write them in the 
following form: 


Yo\ss — Yo\00 = — P, (1.1a) 

Y mies — Ym\o0 = pe”, (1.1b) 

Y m| m= Yolo- (1.1c) 

Here yo is the electric, ym (m=1, 2, 3) the 


magnetic potential, p the charge density and & 
the velocity vector. (Latin indices run from 1 to 3 
and their repetition implies summation.) Finally: 
Yi m=O7y/Ox™: ¥i0=Oy/OX*; (x°=t; c=1). 
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The advanced and retarded potentials are 
solutions of (1.1). The usual procedure is to 
regard the advanced potential as physically 
meaningless and to distinguish the direction of 
the flow of time by taking solutions corresponding 
only to the retarded potential. We shall show 
that in the new approximation method we are 
not faced with the choice between ‘‘past and 
future,’ not, at any rate, in the old form in 
which the choice between advanced and retarded 
potential had to be made. 

Let us assume that the velocity of a particle 
having the charge e is small compared with 
“1” that is, with the velocity of light. Then, if yo 
is of the order ‘‘2’’ (by convention), ym is of the 
order |yo&™"|, that is, small compared with yo; 
we shall regard y,, as of the order ‘‘3.’’ Therefore: 


yo] ~d?; [ym] ~A* (1.2) 


We can devise a method of solving the Eqs. 
(1.1) a method particularly well suited to the 
quasi static case. The first step is to express 
explicitly that the time derivative of yo, ym is 
small compared with yo, ym and with its spatial 
derivatives. To do this we introduce an auxiliary 
time coordinate 

T=)x? (1.3) 


and assume that every field quantity is a function 
of (7, x', x?, x*) rather than (x°, x!, x’, x*). If y’s 
are such quantities then y, y;m, 0y/dr are of the 
same order of magnitude, so that yjo is of the 
order of Ay. Denoting now 


Tia 7 a™ dy/dx™; 


(1.4) 
Vio = Oy /dx° = dy /dT = Ayo, 
we may write (1.1): 
Y0.ss = — PHA*Yo0,00, (1.5a) 
Y mas = pE™ +727 m,005 (1.5b) 
Y m,m=AY0,0- (1.5c) 


Furthermore we shall assume that charges are 
represented not by given charge densities but by 
singularities with spherical symmetry, so that 
the set of equations written above takes, outside 
of the singularities, the form: 


Y0.ss =A*¥0.00, (1.6a) 
Y n.0c = A*¥m,00; (1.6b) 
Y m.m= AY0.0- (1.6c) 


AND GRAVITATIONAL 
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We now expand the y's in power series of \ 
taking the lowest powers of the expansion to be 
of the order determined by (1.2). From the fact 


that only second derivatives of the y's with 
respect to the time enter the equations (1.6a, b) 
it follows that the powers of \ in successive terms 
of the expansion of the y's may differ by “two.” 


Thus we are led to the simple assumption that: 


Yo=NMovotMiyotrA eyot: ++, — (1.7a) 
Y¥m=Ny¥mtMsy¥mt:*:. (1.7b) 

and the equations (1.6) split into: 
21'Y0,ss = 21-2Yo,0 (1.8a) 
214.1Y mss = 21-1Y m,005 (1.8b) 
£141¥m.m=217¥00; /=1, 2:: (1.8c) 


The numbers before the y's denote their order. 
It should be noted that ‘‘zeros’’ after a comma 
change the order; thus, for example, eyo, 9 is of 
the order ‘‘3,”’ evo, oo is of the order ‘‘4”’ etc. 

Let us consider the case of one point singularity 
with its motion represented by £(r). The 
equations (1.8) in the lowest order (J=1) take 
the form: 


2¥0,ss=0; 3¥m.e2=0; 3Y m,m = 2Y0,0- (1.9) 


The solution of these equations corresponding 
for £"=0 to the static solution is: 


oyo=e/r; s¥m= —e/ré™; 


(1.10) 


Indeed, oyo, sym are harmonic functions and 


e¥0,0= —e(0/dx*)(1/r)E*=s7m.m, (1.11) 


so that the equations (1.9) are satisfied. 

The nature of the solution is determined by 
this initial choice of the harmonic functions for 
eyo and sym if we agree to add no arbitrary 
harmonic functions in the further approximation 
steps. We obtain: 


e d2!-2 
21Yo= on ——(r?!-*), 


) — (1.12a) 
(21—2)! dr?! 


9 


e d2!-2 
. . 
21419 m = — ——(r?! se). 


oe (1.12b) 
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and the power series, 


es 1 d?!-? 
vye™meD WO ———— ——_(r**), (1.13a) 
ist (21—2) ! dr? 
x 1 d2!-2 
Yn = —€>d*+1____ —__(p*-8§),  (1.13b) 


l=1 (21—2)! dr? 


if they are convergent, satisfy the electrodynamic 
equations. These power series for yo and y,, are 
well known.2 They correspond to a standing 
wave, that is, to the sum 


} advanced +43 retarded potential, 


whenever the series are convergent. 

The possibility of this solution of the electro- 
dynamic equations was frequently discussed at 
the time of Bohr’s theory. It was intended to 
show that there is no contradiction between 
Bohr’s theory of stationary electron paths 
without radiation and classical electrodynamics. 
Here we have tried only to show that the 
familiar solution can be gained by applying our 
approximation method which, in its simplest 
form, leads to the standing and not to the 
outgoing wave. 

Is it possible to obtain only the outgoing or 
only the incoming wave by the use of our 
approximation method? The existence of such a 
possibility can be shown although, from the 
formal point of view, it seems to be more 
artificial and less simple than the above solution. 

Instead of the power series (1.7) we can take a 
more general expansion, 


Yo=NoVotA sV0+MVOHAM sYO+: °°, 
Ym = A*37 m +447 m +r 5¥m+ laa 


and apply the approximation method as before. 
The odd and even powers of \ do not mix since 
the electrodynamic equations are linear and of 
the second order. The final solution in this more 
general case is characterized not only by the 
choice of 270, xym but also by the choice of the 


(1.14a) 
(1.14b) 


harmonic functions 370, «Ym. We now have 
besides (1.8), the equations: 
214170,ss = 21-17Y0,00, (1.15a) 


2 Nordstrom, Proc. Acad. Amsterdam 22, 145 (1920); 
Page, Phys. Rev. 24, 296 (1924); comp. also Wentzel, 
Zeits. f. Physik 86, 479, 635; 87, 726 (1933). 


INFELD 


= 21Y m,00» (1.15b) 


214+2Y mas 


J=1,2---. (1.15c) 


214.2Y m,m = 214+170,0; 
Let us take for 370, «ym the simple harmonic 
functions, satisfying (1.15) for 7=1: 


syo=0; 4¥m=e(d7é dr’). (1.16) 


Then in the further approximations we find: 


e d?! 1 
wre? ——— —(r?!-?), (1.17a) 
(21—1) ! dr?! 
e d?! 1 
2142Y m = pti 2 (r?! 75). (1.17b) 


(21—1)! dr?!“ 

The general power series in (1.14) if (1.12) and 
(1.17) are inserted, correspond exactly whenever 
they are convergent, to Lenard-Weichert’s 
solution, that is, to the retarded potential. By 
taking only even powers in the expansion for Yo 
(and odd in yn) we obtain the ($ advanced +3 
retarded) potential. Similarly, taking only odd 
powers in the expansion for yo (and even in Ym) 
and starting with (1.16) we obtain a solution 
representing the (3 retarded — 4 advanced) 
potential. By taking both the solutions, that is, 
by inserting (1.12) and (1.17) in (1.14) we obtain 
the retarded potential developed into a power 
series. 

It follows from our argument that for two 
charges e; and é2 moving with velocities 7’, ge we 
should have for eyo, 3¥m; 3¥0, 4¥m 


2Yo>= e:/r1 +¢é2/r2 ’ 


vm — (1.18) 
3¥n= —e rin” — es reo”; 
d?n™ d?—™ 
syvo=O0; 4¥m=e1— r os 
dr dr>— (1.19) 


7,2 = (x*—n*)(x*— 9"); 72? =(x*—f")(x*— $F"). 


§2. GENERALIZATION OF THE GRAVITATIONAL 


CASE 


The method of approximation used in I is 
essentially that presented here in §1. The differ- 
ence consists, however, not only in the fact that 
the calculations are, in the gravitational case, 
incomparably more complicated owing to the 
nonlinear character of the gravitational equa- 
tions. There is another, much more fundamental, 











|.15b) 


l.15c) 


nonic 


(1.16) 


17a) 


.17b) 


) and 
never 
1ert’s 
l. By 
‘Or Yo 
1 +3 
- odd 
1 Ym) 
ution 
iced) 
at is, 
btain 
ower 


two 
r? we 


1.18) 


19) 


INAL 


ffer- 
that 
case, 

the 
qua- 
ntal, 














ELECTROMAGNETIC AND GRAVITATIONAL RADIATION 839 


difference. The motion of the singularities is 
quite arbitrary in the case of the electrodynamic 
field and determined on appropriate choice of 
coordinate conditions in the case of a gravi- 
tational field. We shall quote here, very briefly, 
the results obtained in I with respect to the 
problem of motion of the field singularities. 

In the gravitational case we had 10 gravi- 
tational*® potentials y,,. (The Greek indices run 
from 0 to 3.) The power development used for 
Yu» in I was: 


Yoo= N20 +A *aV00-+A *e¥00 + sti (2.1a) 
Yom = A33Vom tA som + Pe. (2.1b) 
Ymn= Ay mn tr e¥ mnt ° oe, (2.1c) 


In I only odd powers in X for yom and only 
even for Yoo, Ymn were taken. The whole treat- 
ment corresponds, therefore, to that in the 
electrodynamic case for which the solution is a 
standing wave, that is, the sum of the advanced 
and retarded potentials. 

The following set of split equations was in- 
vestigated in I: 


21'Y00,ss = 2arAoo, (2.2a) 
2141Y0m,ss = 22141A0m; (2.2b) 
2141 0m ,m = 21700, (2.2c) 
2142 mn,ss = 22142A mn; (2.3a) 
ei42Amnn =O; J=1,2---. (2.3b) 


The similarity between the equations (2.2) and 
the electrodynamic equations (1.8) is evident. 
The only difference in the left-hand side of these 
two sets of equations is the additional index ‘‘0”’ 
in (2.2). In both cases the right-hand side is 
known and determined by the previous ap- 


. proximation steps. The A’s in (2.2) and (2.3) are 


functions of y’s calculated in the previous 
approximations.‘ (The definition of A’s given 
fully in I is repeated in the footnote.) The main 


* The yyy are not the usual gravitational potentials but 
their linear combinations: yyy = hy» — }nuvn’’hop, where 
huy = Ouv and : Iman = Shan + Ino= 0 : et: S 

* The linear terms are only present in 22;,2A,,- They are: 

—2141°Yom, On — 241Yon, Om +5 nn2t'Voo, 00+ 21Y%mn, 00- 
The nonlinear terms are the result of the development of 
2Apv = 2L uv — nuvn’?’ Lop, 
where Lap is equal to the nonlinear part of Rap. 


idea in the solution of the gravitational equations 
is, as in the electrodynamic case, to bring them, 
by the use of our approximation method and by 
application of the coordinate conditions into the 
form 


where y are the unknown and A the known 
functions. The equations (2.2) (2.3) have been 
derived from the gravitational equations for 
empty space (Ruvy=0) together with the coordi- 
nate conditions and their relation to the gravi- 
tational equations has been investigated. The 
equations (2.2) (2.3) are more general than the 
gravitational equations. They enable us to 
determine the field in every approximation step, 
provided the motion of the singularities is given. 
In (2.2) (2.3) the motion of the singularities can 
be prescribed in an arbitrary way and is not 
restricted by these equations. It turned out to be 
convenient to investigate the equations (2.2) 
(2.3) which represent a more general case than 
that of the gravitational equations, and then to 
ask: what additional conditions must be imposed 
upon the equations (2.2) (2.3) so that their 
solution becomes a solution of the gravitational 
equations for empty space? 

Let us assume, for the sake of definiteness, that 
only two singularities are present, their motion 
being described by the functions n°(r), ¢*(7). Let 
us form, in each approximation step, the integrals 


1 1 
—{ 2ers2A mn cos (n: N)dS =2i42Cm} 
dr 
(2.4) 
1 2 


—| 2o:49A mn cos (n- N) dS =9)40d m. 
4r 


Here the integrals are taken over an arbitrary 
surface enclosing only the first and the second 
singularity respectively, as indicated by the 
number on the top of the integral. Furthermore 
(n-N) denotes the ‘‘angle’’ between the direction 
of x" and the “‘normal”’ to S, and the summation 
convention applies to the m. These integrals do 
not depend on the shape of the surface enclosing 
the singularity, because of (2.3b). Both e:,2¢m 
and 942d, will be functions of r only. We can 
now state the additional condition which must be 
fulfilled in order that the solution of (2.2) (2.3) 
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should also be a solution of the gravitational 
equations for empty space. It is: 


A4GCmaA EC mA sCm+ oo+ =Q), (2.5a) 
Midmt+A ed tA sdmt+- ++ =0. (2.5b) 


There are 6 differential equations for yn’, ¢* and 
their solution gives the possible motions of the 
particles represented by the singularities of the 
field. in this essential point the gravitational 
equations differ from the electrodynamic one, 
because they together with the coordinate con- 
dition determine the motion of the field singu- 
larities. 

By carrying the calculation of A,» up to, say, 
the sixth approximation (as in I), approximate 
equations of motion can be obtained. They are: 


ACmt+A 6m =O, (2.6a) 


AMidmt+A*ed»=0. (2.6b) 


In the first approximation step the equations 
for e700 and sYom are 


2¥00,ss = 0; 3YV0m,ss =(); 


3VO0m,m = 2Y00,0) 


which are exactly the equations (1.9). The 
solution from which we started in I corresponding 
to the existence of ‘wo point singularities with 
central symmetry was 


eyo=229; 2= —2m,,71—2M2/ Po, 


3Vom = 4m, /719" —4m2/ roe", (2.8) 


ry? =(x*—n")(x*— nn); re? =(x*—F*) (x*— F*). 


The constant and positive coefficients:m, and 
mz are the masses of the point singularities. 
Their sign is determined by the fact that two 
masses attract each other. Again we see a 
similarity between (2.8) and (1.18a) (1.18b), the 
only difference being in the constant coefficients. 
The next step was to calculate ;A,, and the 
integrals 4Cm, «dm. As was to be expected, the 
equations of motion up to this approximation, 
that is the equations 


fa=0; daa=0 (2.9) 


proved to be the Newtonian equations of motion. 
We may note for future use that in the same 


INFELD 


approximation we had 


d*n™ d*t* 
Cmt+ad m SS —-+ Me ; = 40; 
dr? TT 


(2.10) 


therefore, m,(d?n/dr*) +mz.(d*¢ dr?) is, at least, of 
the sixth order. 

So far we have tried to repeat briefly those 
results obtained in I which we shall use later. 

Our problem now is to investigate for the 
gravitational case a more general solution in 
which the powers of \ in the development for y's 
do not change by ‘‘2”’ but only by “1.” As in the 
electromagnetic case, we shall seek here a solu- 
tion corresponding not to a standing wave, but 
to an outgoing one. It is indeed possible to 
generalize our procedure by introducing for y’s 
instead of (2.1) the more general development 


Yoo = VN eV00 FA*sV00 +A ay otr svot+ oe es (2.1 la) 
Yom = 8 3Vom FA Vom tA sVom + 7 (2.1 1 b) 
Ymn= NY mn tr sY mat ° (2.11c) 


We shall add the new arbitrary expressions in 
such a way that the arbitrary harmonic functions 
by which we start our calculations are the 
beginning of series representing the retarded 
gravitational potentials and not, as in I, the 
standing wave. It is of some interest to investi- 
gate the influence of such a more general pro- 
cedure upon the equations of motion, or in other 
words, the part which can be played by the 
radiation in the gravitational case, brought in 
through this generalization. In this more general 
case the equations (2.2) (2.3) would have to be 
supplemented by the equations 


214+1Y00,s0 = 22141Ao0, (2.12a) 

214.2Y0m es = 22142A0m: (2.12b) 

214+2Y0m,m = 2141 00,05 (2.12c) 

2143Y mn,se = 22143A mny (2.13a) 
2143Amnn=0; 1=1,2---, (2.13b) 

and instead of (2.5) we should have: 

NMC mA SCmEA 6C mA 7Cmbe* ++ =O, (2.14a) 
Aid mtr sdmtr edmtA Md m+: ++ =0, (2.14b) 
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where 

1 1 

214+3Cm = f Zora ms cos (n- N)dS; 
4r 

(2.15 

1 2 

914.34 » = —{ 2erssA mn cos (n- N)d5. 
4r 


For /=1 the equations (2.12) take the form 
3Y00,e8— 0; 4Y0m,s3 = 0; 
4Y0m,m = 3Y00,0- (2. 1 6) 
For their solution we take, as in (1.19) 
d*n™ d*—™ 
syo=0; som = —4m,—— —4m2.——.._ (2.17) 
dr° dr* 
But in this case we have 
5A mn =O, (2.18) 
and therefore 
sCm=0; sd..=0. (2.19) 
This shows that the radiation terms do not 
influence the equations of motion up to the sixth 
order and that the right hand side of gyom in 
(2.17) is of the sixth order at least. We can, 
however, go one step further and calculate the 


equations of motion up to the seventh order. For 
=2 the equations (2.12) again take the form 


5Y00,ss=0; 6Yom,s2=0; 


(2.20) 
6YOm,m = 5Y00,0- 
We take as a solution of these equations: 
4 d 
Yoo =— ——(mir?-+mer2’), (2.21a) 
6dr 


d 
6Yom = —4—(m yy" +m") 
dr 
4 da 
—- (mri? + morert™). 2.21b) 
6d7 


This solution can be deduced from (1.17) if we 
take into account the different constant coeth- 
cients for the gravitational problem and the fact 
that we are now dealing with a two-particle 
problem. Both functions in (2.21) are harmonic 
and of the desired order because of the equations 
of motion. We can now calculate in turn 27A,. 
and 7m, 7dm. The result is: 


27A i 6YOm — 6Y0n m0 TOm n(s'Y00,00) 
_ (2¢) (sYo0, mn) sae (s 10) (2%, mn) = (3 )(sYo0,m) (2¢,n) 
dis (3) (sY00,n) (2¢, m) + 36 m n(sYoo.1) (29,2) (2.22) 


The calculation of the surface integrals gives: 


4 d*n™ é*t* 
m= — ~My my + ms » (2.23a) 
3 5 


dr dr* 
4 d*n™ de” 

d m= — mam + moe ), (2.23b) 
3 dr* dr* 


But again, because of the equations of motion, 
the right-hand side is not of the seventh, but at 
least of the ninth order. From Robertson’s 
calculation it follows® that the time mean value 
of (2.23) is at least of the eleventh order. The 
result shows the astonishingly small role played 
by the gravitational radiation in the motion of 
double stars. 

I wish to thank Professors Einstein and 
Robertson for their interest in this problem and 
for much stimulating discussion. 


5 Robertson, Ann. Math, 39, 101 (1938). 
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Magnetic Moment of Li’ in the a-Particle Model 


The new method of Rabi! yields values of the nuclear 
magnetic moments of unprecedented accuracy. For Li’, a 
moment of 3.265+0.016 nuclear magnetons has been 
found! while preliminary measurements? of the proton 
moment gave 2.78 with an accuracy of about one percent. 
This gives a difference between the Li? and H' moments of 
0.48, with a probable error of about ten percent. In con- 
trast to this, existing theories based on the one-particle 
(Hartree) model: 4 give only 0.28 to 0.32 for this difference, 
a discrepancy of about 50 percent. It seems therefore worth 
while to re-examine the theoretical situation. 

It seems almost certain that the total angular momentum 
of the Li’ nucleus (3) is composed of a spin } and an 
orbital momentum 1, and that the spin is due to the odd 
proton in the nucleus. These conclusions follow from the 
fact that pairs of like particles have a strong tendency to 
have resultant spin zero. There is considerable direct 
evidence for this tendency in that all nuclei with even 
numbers of neutrons and protons have zero spin. The 
tendency is also understandable theoretically because a 
strong repulsive force between like particles with parallel 
spins is required in order to give saturation of the nuclear 
forces, #.e., to balance the attraction of like particles with 
opposite spin. Direct evidence for the repulsion of like 
particles in the triplet state may be found in the high mass 
of light nuclei of high isotopic number® such as Li’ and 
B!, In the particular case of Li’, the four neutrons and 
two of the protons will give resultant spin zero, leaving 
only the spin of the third proton. Inglist has shown by 
explicit calculation that there is only an extremely small 
probability for a state of affairs in which the two neutrons 
outside the alpha-particle have their spins parallel to the 
total spin and the proton opposite. Therefore, in first 
approximation, the magnetic moment of Li’ should be 
equal to the proton moment plus the moment due to the 
orbital motion. 

The calculation of the orbital magnetic moment requires 
a much more extensive use of a model than does the mag- 
netic moment due to the spin. In the single particle model, 
it is assumed essentially that Li’ consists of an a-particle 
core around which the three other particles are moving in 
a symmetric fashion. Thus the orbital momentum 1 will 
be shared about equally between the three particles, giving 
a momentum } for the proton. Since the neutron orbital 
momentum contributes nothing to the magnetic moment, 
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the orbital magnetic moment will be }, as pointed out by 
Inglis.4 Only very small corrections have to be added to 
this figure which reduce ‘t to about 0.30. 

It may seem more satisfactory to consider the Li’ nucleus 
as made up of an a-particle and a triton (H*). This model 
seems more appropriate in our case than in many others 
because the binding energy of Li’, as compared with 
He*+ Hs, is only 2.76 mMU (milli-mass-units) whereas the 
internal binding energies of H® and Het (as compared with 
free neutrons and protons) are 8.94 and 30.23 mMU, 
respectively. According to Teller,® it is plausible that the 
system of a-particle and triton must have an orbital 
momentum unity around its center of mass, in agreement 
with the experimental spin of Li’ (see above). 

With this model, the magnetic moment of Li’ becomes 
equal to the intrinsic moment of the triton plus the orbital 
moment. The triton moment is certainly very close to the | 
proton moment because the two reutrons in H? will have 
zero resultant spin (see above). The orbital angular mo- 
mentum is shared between triton and a-particle in the 
inverse mass ratio. Since H’ has a specific charge of }, Het 
one of 24 times the proton specific charge, we find for the 
contribution of the orbital motion to the magnetic moment 


of Li? 


po= 47-1343¢-2, - 1755 =0.40;. 


This is in much better agreement with the observed value 
of 0.48 than the result of the single particle model, the 
difference being just within the experimental error. The 
reason for the improvement is that now the a-particle re- 
ceives a share of the orbital momentum, and since its 
specific charge is greater than that of the triton, it gives 
a larger contribution to the magnetic moment. This may 
be considered as some evidence in favor of the a-particle 
model for lightly bound nuclei. For Li® which also falls in 
this class the model would predict exactly the deuteron 
moment. 
H. A. BeETHE 
Cornell University, 


Ithaca, New York, 
April 27, 1938. 


' Rabi, Zacharias, Millman and Kusch; Phys. Rev. 53, 318, 495 
(1938). 

? Private communication from Professor Rabi. 

Rose and Bethe, Phys. Rev. 51, 205, 993 (1937) 

‘Inglis, Phys. Rev., to appear shortly. 

® Feenberg and Wigner, Phys. Rev. 51, 95 (1937). 

Teller and Hafstad, Phys. Rev., to appear shortly. 
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Radial Lines in Laue Spot Photographs 


When taking Laue spot photographs with a simple type 
of Laue spot camera which is described elsewhere,! radial 
lines were noticed apparently emanating from the edge of 
the shadow cast by the lead button which was used to 
prevent halation around the central spot. As the size of 
the button was reduced these radial lines became more 
definite and with a button of 5 mm diameter the result 
obtained with a 7-hour exposure is shown in Fig. 1 which 
was taken at approximately perpendicular incidence to a 
100 plane of a natural sodium chloride crystal. The film 
was 5 cm from the crystal and the x-ray tube a Coolidge 
water-cooled type with tungsten target operated at about 
40 kv peak and 15 ma. The same pattern was obtained 
with a natural potassium chloride crystal but not with 
Iceland spar, mica, zinc blend, fluorite, or gypsum, each 
perpendicular to a natural cleavage face. If the crystal is 
not as nearly perpendicular to a 100 plane for sodium 
chloride’ the same 8 radial lines show up but they begin 
on a rectangle rather than a square. 

The result of a much longer exposure, 15 hours, is shown 
in Fig. 2 for the same crystal as for Fig. 1 but after it had 
been re-mounted some weeks later. The longer exposure 
was undertaken after the negatives had been shown to 
several_men at the New York meeting in February and 
one of them had suggested that it might be due to a sort 
of crossed lattice effect, much as a similar pattern is ob 
tained in the visible region with crossed gratings. Another 
observation on this plate which does not show up in the 


print, is the peculiar “structure’’ of those spots which 
appear exactly round. Instead of being of uniform intensity 
they have the appearance shown in Fig. 3, the rectangular 
portion being much darker and having its long side per- 


pendicular toa radius from the center. Crystals from several 





Fic. 1. Laue spot photograph taken at approximately perpendicular 
incidence to a 100 plane of a natural sodium chloride crystal. 
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Fic. 3. 


sources far removed geographically were supplied by our 
geology department but no change in pattern was observed, 
i.e., no radial lines were exhibited except for sodium chloride 
and potassium chloride crystals, but invariably for these 
In order to obtain comparison photographs sodium chlo 
ride and potassium chloride crystals were grown from a 
saturated solution of the respective salts. Chemically pure 
salts were first allowed to crystalize and these crystals then 
used to form a second saturated solution from which the 
crystals finally used were obtained. These crystals were 
not clear although apparently well formed. Neverthele 
the potassium chloride crystals used showed a pattern 
exactly like that shown in Fig. 1. To date the author has 
had no explanation of these lines nor is he prepared to offe: 
one at present. If they are due to strain why should they 
appear only with the two crystals mentioned and not with 


the others examined. 


A. P. R. WApLUND 


Jarvis Physical Laboratory, 
Trinity College, 
Hartford, Connecticut, 
April 2, 1938. 


‘A. & B. Wadlund, “A Portable Laue Spot Camer “te Am. P 
Teacher, 6, 103 (1938). 
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Comments on the Article by A. P. R. Wadlund: Radial 
Lines in Laue Spot Photographs 


It is well known that the diffraction spots in a Laue 
photograph of a strained crystal are drawn out into radial 
streaks and that new streaks appear as well. The Laue 
photographs reproduced in the article by Mr. Wadlund 
show well-defined spots, indicating that the bulk of the 
crystal is subjected to little or no strain. It seems unreason- 
able therefore that the streaks are due to strain. 

Indeed, it proves possible to interpret the radial streaks 
as due to two-dimensional lattice effects. The directions of 
the diffraction maxima from a two-dimensional lattice are 
given by: u—uo=AL//ibit+-/hb.]+pk. Here up and u 
designate the directions of incident and scattered beams, 
a; and @ are the lattice vectors, b; and be their reciprocals, 
while k is a unit vector normal to the lattice plane. \ is 
the wave-length, JJ; and J/, two integers. The quantity Pp 
has to be fixed so that the right-hand side of the equation 
equals |u—wo!. 

If a surface layer of a sodium chloride cube is normal to 


the incident beam, the equation becomes: 


2 7} 
(H+ He | k 


a 


‘ , 
u=~--[/fait+Ji.a,]+] 1— 
a? 

; The two-dimensional 
structure factor is: Fu,u,=2Lfna+(—1)""'fel] with 1, +/h 
even. Thus the intensity of the diffration effects will be 

very small unless //; and //, are both even integers. 

It is readily seen from the second equation that the eight 
reflections with J/,//. assuming the values (+2,0) (0,+2) 
(+2,+2) will give rise to eight spots (A constant) on the 
photographic plate. The eight spots will lie at the corners 
of a square (J)//,= +2,+2) and at the midpoints of the 
edges of the square (//,//,= +2,0 and 0,+2).! As the wave- 
length increases the eight spots move radially away from 


with a;-a2=0 and a;=a.=5.63A. 


the central spot and thus produce the group of eight inner 
streaks. (See Fig. 1 of Wadlund’s article.) The edge of the 
square formed by the lower wave-length limits of the eight 
streaks is 1.1 cm. This gives a value of 0.31A for the 
minimum wave-length present in the incident beam. This 
value agrees remarkably well with Ajjn=0.30A obtained 
by interpretation of the regular Laue pattern. (Mr. Wad- 
lund states that a peak voltage of 40 kv was used. This 
corresponds to Amin =9.31A, which is slightly too high.) 

The much weaker streaks lying in the region of the photo- 
graphic plate where the regular Laue spots occur, are not 
sufficiently well defined in the reproductions to permit 
quantitative measurements. It is clear, however, that the 
six surface layers of a NaCl cube will give rise to a dis- 
tribution of radial streaks of the type observed. 

It may be stated with confidence that an exact treatment 
of the x-ray diffraction from a finite crystal must give rise 
to boundary solutions corresponding to diffraction from 
the boundary layer. Usually these effects are too weak to 
be observed. The reason for their occurrence in the case of 
NaCl and KCI may be sought in the pronounced mosaic 
character of these crystals. It is well known that these 
crystals may be considered as consisting of a huge number 
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of independently scattering crystal blocks in nearly parallel 
orientation. The great number of discontinuities naturally 
means a corresponding enhancement of the surface effects. 
University of Chicago, W. H. ZACHARIASEN 

Chicago, Illinois, 

April 18, 1938. 

1 The edges of the square are actually hyperbol i¢, as a simple con- 
sideration shows. 





Gamma-Radiation from Proton Bombardment of Carbon 


The resonance capture of protons by carbon has recently 
been investigated by Dee, Curran and Petrzilka.' Of the 
480 
associated with the capture of protons by C” 
subsequent emission of 2.6 Mev y-radiation and the other 


two resonances observed, the one at Kev is to be 


with the 
at about 560 Kev to the capture by the less abundant 
isotope C'* with the formation of N“. The assignment of 
the first resonance can be made unambiguously since the 
positron emission from the residual nucleus N'™ can be 
observed.2?, Moreover the y-ray energy agrees fairly well 
with the expected energy release (binding energy of proton 
in N'=1.9 Mev 
was carried out by the absorption method of Bothe and 
Becker.’ A similar analysis of the 2.6 Mev ¥ 
from Th C” showed practical coincidence between this 


The measurement of the y-ray energy 
radiation 


absorption curve and that from the radiation which is to 
be attributed to C®+H'—N'%+/y. This radiation is very 
probably monochromatic or composed of two lines of 
almost the same energy (cf. reference 2). 

The assignment of the 560 Kev resonance to the proton 
capture by C® is based essentially on the fact that the 
y-ray energy in this case is much greater than the energy 
release (including proton energy) in the C® capture. From 
the fact that with 11 mm of Al the absorption curve 
merges with the continuous background, Dee and _ his 
collaborators deduce a y-ray energy of 7.4 Mev, whereas 
the energy release in the reaction 


C84H! 


+560 Kev—N'"'+ hy 


is 8.2 Mev. However the use of the end of the absorption 
curve for the determination of the y-ray energy is decidedly 
questionable since the absorption should approach zero 
intensity asymptotically. Moreover the half-value thick 
ness (thickness required to reduce the intensity to 50 per- 
cent) is only 2 mm Al and from an analysis‘ of the trans 
mission of the electrons through the absorber an absorption 
curve with this half-value thickness should correspond to 
a y-ray energy of only 4.6 Mev, cf. Fig. 1, assuming that 
the radiation is monochromatic. 

It is almost certain, however, that the radiation is not 
monochromatic. The absorption curves for monochromatic 
radiation consist of an almost straight portion arising 
from energy loss with negligible scattering for small ab- 
sorber thickness, and then a “‘tailing off’? due to diffusion 
of the electrons at larger thicknesses. With increasing 
y-ray energy the straight portion of the absorption curve 
will be relatively longer. An examination of the observed 
absorption curve shows that it is composed essentially 
the presence of at 


of two straight portions indicating 


least two y-rays. 
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Fic. 1. Halt-value thickness in aluminum for the absoprtion of 
Compton and pair electrons generated by monochromatic y-rays of 
energy Av. The inset gives the curve for small energy in greater detail 
The experimental data for Th C”’ (2.6 Mev) and H? (2.2 Mev) were 
obtained by Fleischmann (reference 3) and for Ra C (1.8 Mev) by 
Mitchell and Langer (reference 7). 


Further support of this interpretation is obtained from 
the observation of y-rays in the bombardment of carbon 
by deuterons which is to be ascribed to the reaction® 


e C8+H*+N'4+4n!, 


This furnishes evidence for the existence of an excitation 
level of N'* at about 4 Mev. We may then expect three 
y-rays, one of 8.2 Mev corresponding to the transition to 
the ground state and two y-rays of about 4 Mev corre- 
sponding to the cascade process involving the 4 Mev 
excitation level.° We shall make the plausible assumption 
that no other y-rays are present with appreciable intensity. 

The observed absorption curve should then be a super- 
position of the three individual absorption curves weighted 
according to the intensities of the y-rays. Inasmuch as the 
energies Of the y-rays differ sufficiently so that the indi- 
vidual absorption curves do not interfere strongly, the 
relative intensities of the components of the radiation may 
be found by extrapolating to zero thickness the straight 
line portion corresponding to the hard component. In this 
way we find the relative intensity 2:1 of soft to hard 
component. 

In conclusion it may he pointed out that while the 
half-value thickness versus y-ray energy curve given in 
Fig. 1 is limited in its application to monochromatic 
radiation, in some cases it may be used as above in con 
junction with disintegration data as a rough indication of 
the number of lines and relative intensities involved in the 
y-ray spectrum. 

Cornell University, M. E. Ros 
Ithaca, New York, 

April 26, 1938. 


1P. I. Dee, S. C. Curran and V. Petrzilka, Nature 141, 642 (1938). 

2 Such activity associated with a resonance peak at 480 Kev has also 
been observed by Hafstad and Tuve, Phys. Rev. 48, 306 (1935) who 
observed in addition a resonance peak at 400 Kev proton energy 

+ Bothe and Becker, Zeits. {f. Physik 76, 421 (1932); see also Fleisch- 
mann, Zeits. tf. Physik 103, 113 (1936) 

4H. A. Bethe, M. E. Rose and L. P. Smith, Proc. Am. Phil. Soc, 78, 
573 (1938). 

*Tuve and Hafstad, Phys. Rev. 48, 106 (1935); Bonner and Bru- 
baker, Phys. Rev. 50, 308 (1936); cf. also Livingston and Bethe, Rev 
Mod. Phys. 9, 245 (1937), §101, B and C. 

® We assume that selection rules do not forbid these transitions. 

7 Mitchell and Langer, Phys. Rev. 52, 137 (1937) 
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A Low Temperature Thermal Source of Li Ions 


It was found by Blewett and Jones! that the 8-eucryptite 
(LixO-Al,O 3-2SiOs) they had prepared had an emissivity 
of Li ions about twice as large as that of spodumene, the 
well-known Li* source. Rumbaugh and others? * working 
with Kunsman catalysts were also able to get a much 
larger Li* current than with spodumene. Another type of 
complex compound, which occurs in natural minerals, 
known under the name amblvgonite (LiAI(F, OH)PO, 
Was tried in a mass spectrometer equipped with an elec 
trometer tube amplifier. It was finely powdered and, with 
water as the binder, was mounted on a thin platinum cup 
welded to a tantalum strip filament. Li ions were emitted 
at a temperature as low as about 700°C with a total current 
of about 10~* ampere per cm*. The Li’ current obtained at 
this temperature was more than 10 times larger than that 
obtained from spodumene at about 1000°C. It is believed 
that this natural compound will be suitable as a low tem 
perature Li* source for large scale separation of the 
isotopes. 

Since the two isotopes were completely resolved, it was 
possible to observe simultaneously the time change of the 
emission of each individual ion. In a continuous 12-hour 
run, it was observed that while the total current remained 
fairly steady, the Li emission (peak value) of either kind 
increased slowly with time for the first four hours and then 
more rapidly to a maximum. After that the emission de 
creased gradually with time and apparently attained a 
final constant value. The curves in log-scale are shown in 
Fig. 1. One can only speculate on what causes the initial 
rising emission. The latter decay part of the curve is prob 
ably due to the gradual withdrawal of moisture from the 
interior of the salt. In another run without water binding 
the highest emission was reduced by about } and reached 
the final value earlier. That moisture enhances emission 
has been found by Birkenberg.* 

rhe variation with time of the abundance ratio of the 
two isotopes of Li is of interest. These time change curves 
show that the ratio Li’/Li® increased rapidly with time in 


10°] | sa ] 
f | 
4 he 
y lb A Rin 
i >>} + yo 
co y ; f a 
2 Py / x 
Ww 
: f ‘ 
« ‘ A Lig * 
= é A —_— 
Fic. 1 10,14 F 
6 ] x 
| f 
| f LI ION CURRENTS 
if FROM 
10.44 AMBLYGONITE 
} AT ~ 700°C 
| 
LI 4 i 4 4 ‘ 4 rl 4 i } 
0 ! 2 3.6UC«4 5 6 7 8 9 0 tt t2 
TIME IN HOURS 











846 LETTERS TO 


the first hour, but then attained a fairly constant value for 
several hours, and then gradually increased again. The 
latter part is, therefore, in agreement with that observed 
by Brewer® from an impregnated filament. The ratio of 
the areas under the time change curves is 12.29 for Li*/Li®. 


Horr Lt 


Department of Physics, 
University of Minnesota, 
Minneapolis, Minnesota, 
April 26, 1938. 


' Blewett and Jones, Phys. Rev. 50, 464 (1936). 
* Rumbaugh, Phys. Rev. 49, 882 (1936) 
Smythe, Rumbaugh, and West, Phys. Rev. 45, 724 (1934). 
‘ Birkenberg, Ann. d. Physik 1, 157 (1929). 
* Brewer, J. Chem. Phys. 4, 350 (1936). 





Note on the Range of Occurrence of Stable Isotopes 


If we arrange the known stable (nonradioactive) isotopes 
of the elements in an A—2Z, Z (Harkins) diagram, we 
notice that for each value of the isotopic number A —2Z, 
the isotopes extend in almost unbroken sequence from the 
lowest to the highest atomic number Z for which stable 
isotopes have been found. 

For odd values of A —2Z the sequence of atomic numbers 
is continuous, for even values of A —2Z only, even values 
of Z are found for elements Z>7, in accord with the 
known rule that for elements beyond nitrogen no isotope 
of odd atomic number has even atomic weight. There are 
in all only 19 gaps in the sequences, corresponding to 
isotopes which have not yet been observed. In view of the 
large number (about 275) of stable isotopes already known, 
we may predict with considerable certainty that if iso- 
topes corresponding to the 23 gaps are discovered, they 
would also be found to be nonradioactive, defining the 
to correspond to half-lives less 


limit of “radioactivity” 
than about 10" years. It is evident that nonradioactivity 
does not necessarily imply abundance, although the con- 
verse (rarity of short lived elements) is necessarily true. 
On the above basis the following elements should be stable: 
13A37, 1sA%, o3V%, ogCr®!, oo Ti, ogCr5®, osMn'3, o¢Fe5, 
ogFe®™, os Ni*, 339r?9, aoZr3, 43Ma%, aaMa®, s2le!29, 56Ba!4?, 
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stLa'3’, 5,Ce!, ssCe!4!, Ce, s9Pr, 6:11, 6 11'9?. The iso- 
tope i9K* is known to be radioactive as might have been 
expected since it has odd Z and even A. For 3;Rb* the odd 
even rule is obeyed, but because of its radioactivity we 
need not expect stable isotopes between it and 42Mo*” in 
the A —2Z=13 sequence. 

HARRISON BROWN 


Department of Chemistry, 
University of California, 
Berkeley, California, 
March 29, 1938. 





A Relation Between Atomic Number and Isotopic Number 
for Stable Isotopes 


In Fig. 1 the isotopic numbers A —2Z of all stable iso- 
topes are plotted against the atomic number Z. Known 
isotopes are represented by dashes, and the gaps for which 
stable isotopes are expected in the preceding letter by 
circles. In Fig. 2 the value of A—2Z at the midpoint of 
each sequence is plotted against the corresponding atomic 
number Z, for all even values, and in Fig. 3 for all odd 
values of A—2Z. Between cusps which lie at values of 
A—2Z fourteen units apart, with an additional cusp at 
A—2Z=8, the points in Fig. 2 lie on remarkably smooth 
curves. In Fig. 3 the spacing of 14 units between cusps is 
also evident, and in addition the ranges of 14 units are 
broken up into 8 and 6, respectively. In no case do the 
points deviate from the smooth curve more than could be 
accounted for by the lack of a single isotope at the ex 
tremity of a sequence, with the exception of the two 
points marked by crosses. The point at A—2Z=20 in 
Fig. 2 falls on the curve if we exclude ¢:Sm'**. The single 
cross in Fig. 3 corresponds to 3;Rb*’, which has already 
been proved to be radioactive. Possibly 6.5m" is a-radio- 
active, if not it is the only one besides ,Be’ which appreci- 
ably disturbs the smoothness of the curves. . 

G. E. GiBson 
HARRISON BROWN 


Department of Chemistry, 
University of California, 
Berkeley, California, 
March 29, 1938. 
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Long-Lived Radio Cobalt Isotopes 


We have already reported! confirmation of a long-lived 
radio-cobalt produced by irradiating cobalt with neutrons. 
rhis effect was first observed by Sampson, Ridenour and 
Bleaknev,? who gave the half-life as over a vear. Recently 
Risser’ has quoted the figure as 2.0+0.5 years, after ob- 
serving the decay for six weeks. 

We have been following chemically separated radio- 
cobalt samples, obtained from neutron and from deuteron 
bombardment of C.P. cobalt oxide, for about a year and 
find the decay curves not to be simple exponentials. For 
the first 100 days the rate of decay does indeed suggest a 
2-year period, but for the remaining interval the half-life 
is much longer, appearing to be about 10 vears or more. 
When this is well established, subtraction will show the 
initial decay to be due to a much weaker activity of the 
order of several months half-life. A somewhat older deu- 
teron-activated sample of cobalt oxide, but without 
chemical analysis, bears out these observations, as may be 
seen in Fig. 1. 

Absorption curves in Al and in Pb of the 10-vear cobalt 
(chemically separated from a Co+I1) bombardment) are 
shown in Fig. 2. No appreciable difference is found in such 
curves taken a year apart, indicating that the shorter 
period is in fact very weak. The gamma-ray is reduced to 
half-value by 12 grams per cm? Pb, indicating an energy 
of about 1.3 Mev. The range of the electrons is not well 
defined and may indicate two electron groups with ranges 
of 0.03 to 0.06 and perhaps 0.3 gram per cm?® Al, respec- 
tively. These absorption measurements are in substantial 
agreement with those of Risser.* 

It is not yet certain that the shorter-lived and low in 
tensity radio-cobalt is actually formed from cobalt; slight 
traces of iron or nickel in the cobalt oxide could yield radio- 
cobalt with a period of this magnitude, as we have shown 
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from bombardments of iron with deuterons! and of nickel 
with neutrons and deuterons.‘ 

Perrier, Santangelo and Segré® have reported radio 
cobalts of half-lives about 1 month and about 215 days 
obtained from copper filings from the deflector plate of the 
Berkeley cyclotron, the decay having been followed for 
about a vear. Since cobalt should not be formed from 
copper by deuterons, they suggest that stray neutrons on 
copper have produced Co" and Co® by the (”, a) reaction. 
However, they were unable to confirm this, after irradiating 
copper for several weeks with the neutrons from 500 mC 
of Rn+Be. We have looked for such an activity by examin 
ing borings taken from beneath the surface of the copper 
utility target holder of the cyclotron after it had received 
the neutrons from at least 1000 microampere hours of 4 to 
6 Mev deuterons on Be, and must also report negative 
results. We would like to suggest that iron sputters onto 
the deflecting plate from the iron lid and base plate of the 
cyclotron vacuum tank (marks of this are always evident) 
and hence that Perrier et al. are really observing the two 
cobalt periods that are known! to be produced from the 
deuteron bombardment of iron. In fact, their curve is 
identical with ours (for Co from Fe+D) for about a year; 
beyond that, our half-life has lengthened somewhat and 
there is no assurance that it will not continue to increase 
even more, with consequent lengthening of the period of 
the shorter lived component. 

The bombardment of iron for 1 hour with 100 micro- 
amperes of 3.2 Mev protons yields a long-lived activity 
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(as well as the 18-hour Co**) and a slow decay is also found 
after activating manganese for 1 hour with 0.1 micro- 
ampere of alpha-particles at 16 Mev. The present apparent 
half-lives of these samples are about 180 days and 80 days, 
respectively. But since these activities, as well as those ot 
the above-mentioned radio-cobalts from the activated 
nickel’ and iron, apparently have not yet settled down to 
a single rate of decay, it seems premature to attempt to 
analyze them into component periods which might be 
identified with each other or to speculate on their isotopic 
origin. 

We gratefully acknowledge the support of the Chemical 
Foundation, the Research Corporation and the Josiah 
Macy, Jr. Foundation. 

J. J. LivVINGOOD 
G. T. SEABORG 
Radiation Laboratory, Dept. of Physics (J.J. 
Department ot Chemistry (G.T.S.) 
University ot Calitornia, 


Berkeley, Calitornia, 
April 29, 1938 


' Livingood, Seaborg and Fairbrother, Phys. Rev. 52, 135 (193; 
?Sampson, Ridenour and Bleakney, Phys. Rev. 50, 382 (1936 

' Risser, Phys. Rev. 52, 768 (1937 

‘ Livingood and Seaborg, Phys. Rev. 53, 765 (1938) 

’ Perrier, Santangelo and Segré, Phys. Rev. 53, 104 (1938) 





Erratum: Ferromagnetic Anisotropy of Nickel-Iron 
Crystals at Various Temperatures 


I wish to take this opportunity to point out and correct 
a mistake in the discussion of ‘‘Results’’ in my paper in 
the Physical Review." 

The fourth sentence in this section of the paper which 
now reads ‘“The most unexpected result is that the direc- 
tion for easiest magnetization changes not from <111> 
to <100> as iron increases beyond a critical amount, but 
from <111> to <110>,” should read ‘The most un- 
expected result is that the direction for most difficult 
magnetization changes not from <100> to <111> as 
iron increases beyond a critical amount, but from <100> 
to <110.>" 


J. D. KLEIs 
Research Laboratory, 
Fansteel Metallurgical Corporation, 
North Chicago, Illinois, 
April 20, 1938. 


! Phys. Rev. 50, 1178-1181 (1936) 





On the Penetrating Component of the Cosmic Radiation 


In a previous work! we have pointed out, that at a great 
depth (about 800 meters water equivalent) the ionizing 
part of the cosmic radiation consisted only of shower 
particles produced by a nonionizing radiation. We empha- 
sized also our view, that these penetrating nonionizing 
rays are neutrinos, created somewhere in the atmosphere 
or in the earth, in multiple shower processes (as described 
C. Wilson® comes 


in consequence of his measurements to the same conclu- 


by Heisenberg”). In a recent article V. 
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sion, but he even believes that at great depth successively 
two kinds of radiation are perceivable ‘‘heavy electrons’’ 
and “‘neutrinos.”’ In Fig. 1 we can see the intensity varia- 
tion of cosmic-ray showers from 4000 m above sea level 
down to 1000 m water equivalent depth. The points repre 
* We may 


notice bends on the curve 2: at about 10-20 m and a 


sent experimental data of different authors. 


second somewhere near 250 m water equivalent, but the 
place of the second bend cannot be located with great 
accuracy. The first bend may perhaps signify that from 
this depth onward the majority of the observed showers 
are produced by heavy electrons and from 250 m on by 


neutrinos. This interpretation is also in best agreement 


with the results arrived at by V. C. Wilson. 
]. BARNOTHY 


Institute for Experimental Physics, M. Forro 
University of Budapest, = si 
Budapest, Hlungary, 
\pril 13, 1938 


J. Barnéthy and M. Forr6, Zeit Physik 104, 74! (1937) 

2 W. Heisenberg, Zeits Physik 101, 533 (1936 

V.C. Wilson, Phys. Rev. 53, 337 (1938 

'R. HL. Woodward, Phys. Rev. 49, 711 (1936); S. Kikuchi and 
J. Watasa, Proc. Phys. Math. Soc. Japan 18, 210 (1936); P. Auger and 
G. Meyer, Comptes rendus 204, 572 (1937); A. Ehmert, Zeits Physik 
106, 751 (1937 
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1. Radioactivity. Kart K. Darrow, Bell Telephone Laboratories. 

2. Thermodynamical Viewpoints in the Study of Nuclear Reactions. L. W. Norpuerim, 
Duke University. 

3. Twenty-Five Years of Wilson Cloud Chamber Research; an Illustrated Survey of the 
More Important Discoveries Made with the Cloud Chamber. Artuur E. RvuARk, L’niversity 





of North Carolina. 


ABSTRACTS OF CONTRIBUTED PAPERS 


1. Measurement at High Frequencies of the Diamag- 


netic Faraday Effect in Liquids. G. C. Comstock, The 


Citadel.—Direct measurements of the Verdet constants of 
four diamagnetic liquids: carbon disulfide, glycerine, iso- 
propyl and amyl alcohol’, for optical frequencies far 
removed from their absorption regions, have been made in 
the comparatively high radiofrequency range from 1 to 
1.710 c.p.s. by the modification of a method due to 
Bretscher,' in which a photo-cell coupled to a tuned radio- 
frequency amplifier was used as the detector of the rotations 
of the plane of polarization. The variations of the plate 
current in the output circuit were calibrated in terms of 
the peak Faraday rotations. The radiofrequency values of 
the Verdet constants were compared with those obtained 
by measurement at 60 cycles and with direct current using 
a.f. and d.c. photoelectric amplifiers. The Verdet constants 
were found to be independent of the frequency of the field 
up to 1.710% c.p.s. within the experimental accuracy of 
about 5 percent. These measurements are in agreement 
with the theoretical predictions of Persico* of the constancy 
of the Faraday rotation with variable magnetic field 
frequency for optical frequencies far removed from the 


absorption frequencies of the liquids. No anomalous effects 


849 


were found at radiofrequencies for which these liquids show 
strong absorption of the radio waves. 


Ek. Bretscher, Helv. Phys. Acta 9, 42 (1936). 
I. Persico, Rend. Linc. 3, 561, 603 (1926) 


2. Optical Activity of Nickel Sulfate, a-Hexahydrate. 
EDWARD B. NELSON, Vanderbilt University.— Nickel sulfate 
crystals of the a-hexahydrate form have been examined to 
determine the relative number of laevo- and dextro- 
rotatory crvstals, the activity of their solutions, and the 
rotatory dispersion of the crystals. The relative abundance 
of the laevo- and dextrorotatory crystals is a one to one 
ratio. The crvstals are found to be mirror-image crystals, 
since sections of equal thickness give numerically equal 
rotations for any particular wave-length. Aqueous solutions 
of only laevo- and of only dextrorotatory crystals were 
separately examined throughout the visible spectrum, but 
no optical activity was observed. Visual observations with 
a Lippich half-shade polarimeter and photographic ob- 
servations with an adaptation of the method of Fizeau and 
Foucault were made of the rotatory dispersion of the 
crystals which was found to be of an anomalous nature. 
The rotation is zero for light of wave-length 5030A. Above 
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this point the rotation increases gradually in magnitude to 
about 6200A where it begins to increase rapidly. Below 
5030A the rotation increases gradually in magnitude to 
about 4600A where it begins to increase rapidly. The sign 
of the rotation of each crystal changes on passing the zero 
point. The rotation for the D line of sodium is approxi- 


mately 1.55° per mm. 


3. The Effect of Humidity and Irradiation on the 2.0 cm 
Sphere Gap. ArTHUR B. Lewis, University of Mississippi. 

Quantitative data are presented showing the effects of 
humidity and irradiation on the 2.0 cm sphere-gap volt- 
meter in the lower voltage ranges. Preliminary work has 
been confined to a single gap setting, 4.0 mm (about 10 kv 
peak). Irradiating this gap with a 500-watt cored-carbon 
arc lowers the sparking voltage by 2.0 percent from the 
value obtained when the gap is completely dark, and 
decreases the probable error of a single observation from 
about +0.7 percent with the gap dark to about +0.1—0.2 
percent with the gap irradiated. A sixfold increase in the 
intensity of this irradiation causes no further significant 
decrease in either the sparking voltage or the scattering of 
individual observations. The presence of water vapor in 
the gap causes a roughly linear increase in the sparking 
voltage of +0.1 percent in voltage per mm (of Hg) increase 
in vapor pressure. The control of these two variables 
permits the reproduction of the sparking voltage of this 
particular gap to +0.3 percent. This increased accuracy 
of indication, however, reveals the presence of other 
variables which affect the indications of the sphere-gap 
voltmeter and which must be isolated and controlled if 
the full accuracy of the instrument is to be realized. These 


are being investigated. 


4. A Sensitive Differential Barometer. Wittiam Hurst, 
Duke University. (Introduced by W. M. Nielsen.)—A sensi- 
tive differential barometer has been designed and con- 
structed for application in clinical medicine. Its sensitivity, 
small cost, and other characteristics suggest that it may 
have a wider range of usefulness. In place of the more 
conventional metal bellows, a drum with light weight 
rubber stretched over it is used. Motion of the rubber due 
to differential pressure is applied through a train of gears 
from a small jeweled wrist watch to an indicating arm 
attached to the last gear. The device is sensitive to pressure 
changes equivalent to 0.002 cm of mercury and is capable 
of indicating pressure changes as rapidly as 10 pulses per 
second. Because of the design and character of the com- 
ponent parts, the instrument is quite stable. It should be of 
particular usefulness in the measurement of pressure 
variations in systems of small volume. 


5. Separation of Solutions by the Tubular Centrifuge.* 
J. W. Beams, University of Virginia.—The tubular vacuum 
type air-driven centrifuge! has been modified and applied 
to the separation of solutions. The centrifuge ‘‘bowl’’ 
consists of an alloy steel tube which is spun in a vacuum 
chamber. The solution enters the spinning tube at the top 
continuously through a stainless steel tube, gauge 14. The 
ighter and heavier fractions are collected at the axis and 
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periphery respectively by channels in the bottom of the 
spinning tube. These fractions are led to the outside through 
coaxial tubes, gauges 15 and 12, which support the rotor. 
Vacuum tight oil glands seal the vacuum chamber and 
separate the chambers in which the two fractions are col- 
lected. The ratio of the quantities of light to heavy frac- 
tions collected can be varied over wide ranges. The oil 
glands or bearings are maintained at approximately the 
same temperature as the rotor by circulating water or oil. 
Possible stirring where the solution enters the centrifuge 
and where the fractions are drawn off is prevented by 
baffles. Experiments also show that the remixing caused 
by comparatively wide variations in rotational speed could 
be prevented by a coil of screen wire which filled the centri- 
fuge. The observed separation was in rough accord with 
the calculations of Archibald.? 


* Supported by a grant from the Division of Natural Sciences of the 


Rocketeller Foundation 
Beams, Linke and Skarstrom, Science 86, 293 (1937). 


2 In press. 


6. Practical Applications of Electronic Control Principles 
to the Operation of a Wilson Cloud Chamber. I. I. Hop- 
KINs, R. T. DICKERSON AND W. M. NIELSEN, Duke Uni- 
versity.—A circuit for the complete electronic control of a 
Wilson cloud chamber has been designed by one of us 
(R.T.D.), constructed, and put into operation. The general 
idea is based on the maximum possible use of gaseous dis- 
charge valves and typical resistance capacitance time 
circuits. Tracks are illuminated by two capillary helium- 
mercury discharge lamps, each operating from a five micro- 
farad condenser charged to about five thousand volts. 


7. Probable Errors in Counts from Geiger-Miiller 
Counters. Eric RopGers, University of Alabama.—When 
a large number N of random events are counted with a 
Geiger-Miiller counter, the result is usually written as 
N+(N)}!. This is based on the assumption that the events 
being counted follow a Poisson distribution, in which case 
the standard deviation equals the square root of the aver- 
age. If x is assumed to represent the ‘‘true’’ average count 
for a given time interval, the probability that a single trial 
does not deviate from x by nfore than (x)! lies between 
(2/mwe)t and (2/m)! if x is sufficiently large. Thus (x)! is 
perhaps not a bad approximation to the probable error, 
which is that deviation from x for which the chances are 
equal that a single trial yield a result within or without 
the range. However, when a single observation is made, 
the result may differ widely from the “true’’ average for 
that interval. It therefore seems that not much faith can 
be put in the practice of taking the square root of the 
number of counts obtained in an interval and calling it the 
probable error even if the interval is of such length that 


many counts are obtained. 


8. Construction and Use of Lissajous Figures. F. I 
MitTCcHELL, University of Alabama.—Beginning with space 
models of Lissajous figures, their equations may be deduced 
and the component sine waves analyzed without their being 
in phase. The results are checked by the actual figures as 


produced by the oscillograph. 
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9. Hysteresis Testing with Cathode Rays. C. T. Razor, 
The Citadel—Application of the conventional type of 
cathode-ray oscillograph with electrostatic deflection plates 
to draw hysteresis loops of magnetic specimens. No modifi- 
cation of the instrument is needed. Hysteresis tests may 
also be made on the cores of completed machines without 
damage or dismantling. The method is applicable either 
to quantitative measurements or to quick and easy demon- 
stration. The necessary accessories are found in most 
laboratories. 


10. A Device for Controlling the Temperature of Ab- 
sorption Cells. DUDLEY WILLIAMS AND RICHARD TASCHEK, 
University of Florida.—A device for controlling the tem- 
perature of small absorption cells is described. The ar- 
rangement involves a thermocouple, a galvanometer, a 
photronic cell, and a miniature relay. A light beam re- 
flected from the galvanometer mirror falls upon the pho- 
tronic cell when the temperature of the absorption cell 
increases to a certain value. The photronic cell operates a 
miniature relay which serves to increase the resistance in 
series with the heating coil, thus reducing the heater cur- 
rent. By proper adjustment of resistances the temperature 
of the absorption cell can be kept within a certain range. In 
the arrangement used at present it is possible to control 
the temperature of the absorption cell to 0.1 degree C 
at any point between 25 degrees C and 160 degrees C. By 
the use of a more sensitive galvanometer it should be possi- 
ble to reduce the range of variation. This arrangement may 
prove useful in other experimental work where the geom- 
etry of the apparatus precludes the possibility of using 
ordinary forms of thermostatic controls. 


11. A Spectroscopic Study of Deuterium Bonds. WALTER 
Gorpy, Mary Hardin-Baylor College-—Mixtures of deu- 
terium oxide with dioxane and with acetone have been 
studied in the region 3 uw to 12 uw. The effect of dioxane and 
of acetone on the spectrum of the deuterium oxide was to 
shift the intense 4 « band to the shorter wave-lengths by 
about 0.2 un. In the most dilute solutions, however, the 
position of the band appeared at wave-lengths appreciably 
longer than that of the corresponding band of deuterium 
oxide vapor. In the change from the liquid to the solution 
state, the shift to the shorter wave-lengths probably indi- 
cates the breaking of deuterium bonds between the deu- 
terium oxide molecules. The appearance of the band at 
longer wave-lengths in the solution than in the vapor state 
may be a result of the formation of deuterium bonds 
between the solute and the solvent; variation of certain 
bands of the solvent provide additional evidence for this 
interpretation. 


12. Calibration Wave-Lengths for Infra-Red Spectrom- 
eters. P. FE. SHEARIN AND E. K. PLYLER, University of 
North Carolina.—About 20 wave-lengths which have been 
measured on grating spectrometers have been given for the 


calibration of prism spectrometers in region from about 2 4 
to 24u. A number of lines due to HO and CO, can be 
normally observed due to atmospheric absorption. Since 


not many sharp bands have been measured on grating 
spectrometers from 74 to 134, several sharp bands in 
ethyl alcohol and amyl alcohol in the liquid state have been 
measured in this region. The values found in the present 
work check well with those previously found by Wenigetr 
except those in the 9 u and 10 u region. Here some of the 
bands differ by about 0.2 « from those found in this work 
This difference is probably due to the fact that Weniger did 
not have the more accurate values of the index of refraction 
of rocksalt due to Paschen at his disposal. 


13. A Portable Eagle Mounting for a 3 m Grating.* 
HERTHA SPONER AND WILLIAM Hurst, Duke University 
Confronted with the necessity of building a grating spectro 
graph which would take little room and could be easily 
moved, we have tried the following construction of an 
Eagle mounting of relatively light weight. The instrument 
consists of two independent parts. One is the actual spectro 
graph having an | beam chassis, supporting an optical 
bench for movement of the grating, and the plate holder. 
This mechanism is surrounded by a double wall insulating 
box. The second part is a steel carrier on which the instru 
ment rests on three points, at a convenient height above 
the floor. The carrier has three rubber tired wheels which 
may be removed and easily replaced for moving the instru 
ment around the building. Mr. F. B. Key adjusted the 
instrument. Results have proven so far to be satisfactory 
with no broadening of the lines even when the instrument 
was subjected to light blows in short exposure experiments 
It should be added that the floor of the room where the 
instrument is now set up is floated. Another important 
feature of this apparatus is that the construction was 


relatively inexpensive. 





* Read by title. 


14. X-Ray Crystal Structure Study of Tetramethyl Am- 
monium Tri-Iodide. KR. C. L. Moonry, Newcomb College, 
Tulane University.—Tetramethyl ammonium tri-iodide was 
prepared by a method analogous to that used for ammo- 
nium tri-iodide; crystals of suitable size and perfection were 
obtained from slow cooling of an alcoholic solution. They 
were definitely identified as the tri-iodide by quantitative 
chemical analysis. X-ray diffraction data were obtained by 
means of the single crystal oscillation method. The complete 
set consisted of six overlapping fifteen degree oscillations 
about each axis. The Patterson modification of the Fourier 
series method was used to locate atomic positions. The 
cell size of the crystal is 10.43A, 8.60A, 14.62A along the 
a, b, and c directions, respectively. The class is orthorhombic 
bypyramidal, with the space group symmetry of Q,” 
(Pnmn). All of the iodines are in special positions, lying 
in the reflection plane perpendicular to the y axis. The 
six parameters necessary to locate them in the plane were 
found by consideration of observed diffraction intensities. 
Some assumptions of reasonable interatomic distances 
were required to decide the positions of the cations. Three 
iodines are associated to form a monovalent complex anion 
group. The group is linear, with an iodine to iodine distance 
of 2.75A. A comparison of the details of this structure 
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with those of ammonium tri-iodide previously examined 
demonstrates the effects of the relative sizes of anion and 


cation. 


15. Infra-Red Absorption Spectrum of a Liquid Crystal. 


RICHARD TASCHEK AND DUDLEY WILLIAMS, University of 


The infra-red spectrum of the anisotropic liquid 
been 


Florida. 


(peri-ethoxy-benzalamino) cinnamate has 


peri- 
studied at temperatures from 20°C to 160°C, in the region 
between 2u and 12. Intense characteristic absorption is 
found at 3.3u, 5.85, 6.254, 7.74, 8.0u, 8.6u, 9.74, and 11.6u. 
All but the 8.0u and 8.64 bands can be assigned to known 
inner vibrations of the groups occurring in the molecule. 
In the anisotropic liquid phase the 3.34 and 7.74 bands are 
found to be shifted to slightly longer wave-lengths than in 
the solid or isotropic liquid. At constant wave-length, the 
transmission as a function of temperature remains constant 
while the compound is solid. Between the melting point 
and clearing point the transmission increases, the rate and 
continuity of increase depending upon the wave-length. 
The transmission remains constant beyond the clearing 
point and shows a marked hysteresis on cooling the com- 
pound. For most wave-lengths one or more “‘arrests’’ are 
temperature- 
with 


temperatures in the 
transmission curves. The increase in 
temperature is attributed to lessened scattering from the 


observed at definite 


transmission 


constituent ‘‘swarms’”’ of the anisotropic phase as it passes 


to the clear liquid. 


16. The Infra-Red Spectra of Nitrate Solutions. DUDLEY 
WILLIAMS AND LORAINE DECHERD, University of Florida.- 
The absorption of several nitrate solutions has been meas- 
ured in the region between 2u and 16u. The results are com- 
pared with earlier measurements of nitrate absorption. In 
the solutions the nitrate bands near 7.5u and 15.5 are 
intense, while the bands near 3u, 4.8u, and 12u, are com- 
paratively weak. It is found that the characteristic fre- 
quencies of the NOs; group in solution are slightly lower 
than those of the NO; group in crystals, a result in accord 
with Pfund’s early reflection measurements. 


17. Design and Calibration of an Infra-Red Spectrom- 
eter with a Fluorite Prism. E.S. BARR, Tulane University. 
—An infra-red spectrometer of the Wadsworth type, with 
an aperture of 0.218, has been built for use with a 60° 
fluorite prism (sides 45 by 75 mm.). Arrangements have 
been made for observing prism settings by means of a 
calibrated dial, or directly to 10” from a circle. There is 
provision for interchanging prisms to vary the range of the 
instrument. Using the index of refraction values of fluorite 
given by Coblentz,! tables have been prepared for circle 
readings of angles of incidence and for effective slit widths 
as functions of the wave-length and are shown in the form 
of curves. By the use of a secondary prism, the fluorite 
prism is set at minimum deviation for sodium light. The 
Wadsworth mirror was set visually for sodium light from 
a sodium arc lamp. For finer adjustment of the mirror and 
for calibration, atmospheric absorption bands due to CO, 
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and H:O were used. Other well-established points from 
spectra measured by grating instruments serve as a check 
on the correctness of the calibration. 


W. W. Coblentz, J. Opt. Soc. Am. 4, 432 (1910 


18. The Refractive Index of Air in the Visible and 
Photographic Infra-Red. Davip BENDER, California Insti- 
tute of Technology.—The refractive index of dry CO.-free 
air has been measured from 5300A to 11,200A by means of a 
Fabry-Perot interferometer with invar separating posts. 
The following dispersion curve was obtained for standard 


conditions (0°C and 760 mm 


15.77 0.421 
(to—1)107 = 2885.59 + = 
d? r4 

where \ is expressed in uw. A comparison of the present re- 
suit with those of other observers shows that this curve lies 
above the others by about 10 parts in 10’, and that all 
curves are nearly parallel. This suggests that the dis- 
crepancies are due to svstematic errors, providing, as seems 
to be the case, that the constitution of air remains suffi- 
ciently constant. In the present experiments systematic 
errors due to phase change at reflection, methods of purifi- 
cation of the air, and bulk modulus of the invar posts were 
considered. A calculation concerning the oxygen absorption 
band at 7600A indicates that this band has no measurable 
effect on the dispersion of air. 


* Now at Louisiana State University. 

19. The Polar Moment of Cellosolves. W. H. Byers 
AND A. A. BLEss, University of Florida.—The polar moment 
of 2-methoxy ethanol and of 2-ethoxy ethanol was deter- 
mined using the heterodyne beat method for finding the 
dielectric constant. With benzene as a solvent the moment 
is 2.05 X107'8 e.s.u. for the first compound and 2.13 X107'* 
e.s.u. for the second. The relation between these values and 
the values of the moment of the related alcohols is dis- 
cussed. Work is now in progress on the other members 


of the series. 


20. A Sensitive Method of Detecting Magnetic Im- 
W. Constant AND J. M. Formwatt, Duke 


purities. F. 
The presence of ferromagnetic impurities in 


University.- 
any small solid specimen may be detected and magnetic 
moments per cc as small as 1077 (compared with 10° for 
pure iron), measured by the following method. The speci- 
men is first placed in a strong magnetic field which serves to 
magnetize all the magnetic impurities in the same direction. 
Their residual magnetism is then measured by suspending 
the specimen by a delicate quartz fiber at right angles toa 
weak field supplied by Helmholtz coils. By noting the 
deflection and calibrating the fiber the residual magnetic 


9 


moment is determined. It was found to range from 7 X10 
per cc for brass, to 31077 per cc for Pyrex, all materials 
tested showing some ferromagnetic impurity. It is of 
interest that such small impurities are still ferromagnetic; 
further study of their magnetic properties is being made. 
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21. A Displacement Polarimeter. Newron UNDERWOOD, 


Vanderbilt University —The device is based on a deforma- 
tion of the fringe system seen when a uniaxial crystal cut 
with its faces perpendicular to its optic axis is observed in 
convergent elliptically polarized light. The polarimeter is 
constructed using Polaroid, quartz and Cellophane, ar- 
ranged in the following order; the polarizer, the specimen, 
the wave plate which covers one-half of the field, the con 
verging lens, the quartz, and finally the analyzer Polaroid. 
The field of the instrument consists of the rings due to the 
quartz and is bisected by the edge of the Cellophane wave 
plate. The rings form continuous smooth curves when the 
axis of the wave plate is parallel to either the analyzer or 
polarizer. A rotation of the plane of polarization produced 
by inserting a specimen immediately in front of the polar- 
izer causes a break to appear in the curves at the edge of the 


x 
ean 
ow 


wave plate. Rotation of the polarizer or analyzer to compen- 
sate for the rotation introduced by the specimen causes the 
curves to line up again and thus affords a measurement of 
the rotation due to the specimen The measurements are 


accurate to 0.1 

22. Collision Relations. Akrutir | RUARK AND 
CreiGuron Jones, University of North Carolina.— Braun 
bek! has published formulas dealing with energy and 
momentum relations in the collision of two unequal rela 
tivistic particl His results cover both radiative and non 
radiative collisions. We have worked out more formulas, 
with a view to cloud-chamber applications, and will present 


a systematic summary of those which appear most usetul 


1 Zeits. ft. Physik 96, 600 and 104, 619 





Organization Meeting of the 
New York State Section 


of the 


American Physical Society 


UNION COLLEGE, SCHENECTADY 
April 2, 1938 


INAL steps were taken in the organization of the New York State Section 

of the American Physical Society at a meeting held at Union College, 
Schenectady on Saturday, April 2. One hundred seventy-five persons attended, 
representing institutions in all parts of the state (exclusive of the metropolitan 
area) and including a large group of High School physics teachers as well as 
physicists in college and industrial work. After the constitution and by-laws 
were adopted 83 persons enrolled as charter members and paid their annual 


dues. 


The following offcers were elected for a two year term: 


Chairman, 


P I. Wop, Union College 





Vice Chairman, W. B. Rayton, Bausch & Lomb Optical Co. 


> 


Secretary, P. R. GLEASON, Colgate Universit, 


Treasurer, G. H. CAMERON, Hamilton College 


Members of Executive Committee: 


Terms expiring in 1942: 


Mrs. ANNA W. PEARSALL, Hamilton High School 


R. E. BURROUGHS, Eastman Kodak Co. 


Terms expiring in 1940: 


R. C. Gispps, Cornell University 


C. A. Moose, Milne High School, Albany 


At the dinner meeting Professor Gibbs brought official greetings to the new 
organization as President of the Optical Society of America; Dr. Irving 


Langmuir brought similar greetings from the American Chemical Society. 
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The program of invited papers was as follows: 
Morning Session 


Address of Welcome. D). R. Fox, President, Union College 
Television Transmission. E. H. B. BARTELINK, General Electric Co. 
Biological Effects of Cathode Ray Irradiations. C. P. Haskins, L’nion College. 





Afternoon Session 


Recent Developments in Glass Manufacture. H. P. GaGr, Corning Glass Works. 
Physics in Aeronautical Science. P. E. HEMKE, Rensselaer Polytechnic Institute. 
Electron Optics. ne JOHNSON, General Flectric Ce. 
PauL R. GLEASON 
Secretary 




















